3D Printed Cobalt-Chromium-Molybdenum Porous Superalloy with Superior Antiviral Activity

COVID-19 pandemic and associated supply-chain disruptions emphasise the requirement for antimicrobial materials for on-demand manufacturing. Besides aerosol transmission, SARS-CoV-2 is also propagated through contact with virus-contaminated surfaces. As such, the development of effective biofunctional materials that can inactivate SARS-CoV-2 is critical for pandemic preparedness. Such materials will enable the rational development of antiviral devices with prolonged serviceability, reducing the environmental burden of disposable alternatives. This research reveals the novel use of Laser Powder Bed Fusion (LPBF) to 3D print porous Cobalt-Chromium-Molybdenum (Co-Cr-Mo) superalloy with potent antiviral activity (100% viral inactivation in 30 min). The porous material was rationally conceived using a multi-objective surrogate model featuring track thickness (tt) and pore diameter (ϕd) as responses. The regression analysis found the most significant parameters for Co-Cr-Mo track formation to be the interaction effects of scanning rate (Vs) and laser power (Pl) in the order PlVs>Vs>Pl. Contrastively, the pore diameter was found to be primarily driven by the hatch spacing (Sh). The study is the first to demonstrate the superior antiviral properties of 3D printed Co-Cr-Mo superalloy against an enveloped virus used as biosafe viral model of SARS-CoV-2. The material significantly outperforms the viral inactivation time of other broadly used antiviral metals such as copper and silver, as the material’s viral inactivation time was from 5 h to 30 min. As such, the study goes beyond the current state-of-the-art in antiviral alloys to provide extra protection to combat the SARS-CoV-2 viral spread. The evolving nature of the COVID-19 pandemic brings new and unpredictable challenges where on-demand 3D printing of antiviral materials can achieve rapid solutions while reducing the environmental impact of disposable devices.


Introduction
Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is the third human coronavirus [1][2][3][4] that is much more contagious than SARS-CoV and MERS-CoV [5][6][7][8][9][10][11][12]. Its rapid transmission rate has provoked the current coronavirus disease (COVID-19) pandemic. SARS-CoV-2 is a highly pathogenic enveloped positive-sense single-stranded RNA virus [13][14][15] that belongs to the Baltimore group IV [16]. This global life-threatening situation needs the development of new antimicrobial approaches that could treat or prevent COVID-19 infections [17][18][19][20][21][22]. In this regard, non-woven fabrics are currently used for the fabrication of infection prevention clothing, such as face masks, caps, scrubs, shirts, trousers, disposable gowns, overalls, hoods, aprons, and shoe covers. These infection prevention tools are needed, especially in hospitals during surgical operations, in microbiological and biomedical biosafety laboratories, and face masks should be worn by most citizens as a demonstrated prevention tool in the current COVID-19 pandemic [23]. Nevertheless, these items of clothing for infection prevention are produced with materials that do not possess antimicrobial properties, and some progress has been achieved thus far in the development of antimicrobial prevention fabrics [24].
Many disinfectants such as household bleach, hand soap solution, ethanol, povidoneiodine, chloroxylenol, chlorhexidine and benzalkonium chloride have shown potent antiviral activity against SARS-CoV-2 so far [25]. Thus, non-woven fabrics for the fabrication of face masks and face shields have been treated with benzalkonium chloride to produce antimicrobial infection prevention tools [26,27]. Very recently, non-woven fabrics have been coated with solidified hand soap to produce antimicrobial face masks capable of inactivating SARS-CoV-2 in one minute of contact [28]. Other authors have produced next generation infection prevention materials using other antimicrobial agents such as antimicrobial polymers, salts, carbon nanomaterials, metals or metal oxides [24,[29][30][31]. However, very few studies about the use of natural and biodegradable compounds such as cranberry extracts have been reported [24].
Cranberry extracts have shown antiviral activity against other enveloped viruses such as the herpes simplex virus type 1 (HSV-1) and type 2 (HSV-2) due to the presence of antimicrobial A-type proanthocyanidins (PACs) that provoke alterations of their envelope glycoproteins. However, although HSV-1 and HSV-2 belong to a different Baltimore group I [16] than SARS-CoV-2 because they are double-stranded DNA viruses, they are also enveloped viruses like SARS-CoV-2. A cranberry extract has also exhibit antiviral activity against influenza virus (IFV) [32]. IFV is a negative-sense single-stranded RNA virus that belongs to the Baltimore group V [16]. However, it is also enveloped like HSV-1 and HSV-2. Therefore, since it seems that the PACs present in cranberry extracts effectively interact with the envelope glycoproteins achieving viral inhibition, we hypothesize here that a commercial non-woven fabric treated with two different commercial extracts produced with different cranberries will show antiviral activity against the enveloped SARS-CoV-2 and phi 6 viruses. The phage phi 6 is also an enveloped double-stranded RNA virus (group III of the Baltimore classification [16]) that can be used as surrogate of SARS-CoV-2 and other enveloped viruses such as influenza due to biosafety reasons [26]. Furthermore, atypical viral pneumonia is associated with SARS-CoV-2 infection [2,33] that can increase its risk by co-infection with Gram-positive bacteria [34][35][36][37], including clinically relevant antibioticresistant strains. Additionally, bacterial resistance to pneumonia treatments is increasing at an alarming rate [38,39]. Since the PACs present in cranberry extracts are well-known for their antibacterial properties against Gram-negative Escherichia coli [40] and antifungal activity against Candida albicans [41], we hypothesize here also that the two non-woven fabrics dip-coated with cranberry extracts will show also antibacterial activity against two Gram-positive multidrug-resistant bacteria, the methichillin-resistant Staphylococcus aureus (MRSA) and the methichillin-resistant Staphylococcus epidermidis (MRSE).
The existence and evolution of SARS-CoV-2 coronavirus variants are spreading more effectively than earlier ones [42,43] due to its evolutionary advantage over SARS-CoV primarily in enhanced host receptor binding [44]. According to Sun and Ostrikov [45], this gives the virus a longer-lasting ability to retain activity on diverse surfaces. As such the lessons so far dealing with SARS-CoV-2 transmission suggests that the new developments of antiviral materials are critical in establishing pandemic preparedness [45,46].
The development of effective materials with antiviral capacity is critical in the development of reusable protective devices such as mask filters, high-efficiency particulate air filters (HEPA) and other antiviral devices [47]. Antiviral portable HEPA filters, air purifiers and aerosol decontaminants are critical also for hospital isolation wards and temporary anterooms [48]. Currently, contaminated devices such as masks and filters are disposable under strict protocol, which is often ignored and possess a serious risk of secondary contamination [49]. A recent study by Maclntyre et al. [50] confirmed that even decontamination of reusable masks through washing still offer the probability of infection; further highlighting the need for antiviral materials. Careless disposal of masks is another issue that is resulting in a potential source of microplastic pollution [51] and environmental damage threatening aquatic and animal habitats [52][53][54].
The development of antiviral face mask that contain porous filters or antiviral face shields made of antiviral materials capable of inactivating SARS-CoV-2 is an effective strategy to provide extra protection against the COVID-19 transmission [26][27][28]55] and produce reusable devices to reduce environmental impacts [54,56]. The COVID-19 pandemic has also highlighted the potential for supply chain disruptions causing shortages of essential supplies including PPE, swabs, and ventilators [57][58][59]. As such, there is a requirement for on-demand and onsite manufacturing of antiviral materials suitable for a range of applications. Digital manufacturing techniques such as additive manufacturing (3D printing) demonstrated in this study offers significant potential making on-demand and onsite fabrication of antiviral devices accessible [58,60]. Additive manufacturing (AM) is transforming medical supplies by allowing personalisation and onsite fabrication enhancing resilience against supply chain disruption [61][62][63][64]. This subsequently results in rapid development and deployment of potential solutions which is critical when it comes to pandemic preparedness [65][66][67]. Nevertheless, achieving this requires the development of antimicrobial materials [68][69][70] and processes that can be additively manufactured without the requirement for complex pre/post-processing as demonstrated in this study.
According to Doremalen et al. [71], aerosol and surface transmission of SARS-CoV-2 facilitates infection as the virus can persist for extended periods in a range of common materials as listed in Table 1. Kumar et al. [72] highlight that not many studies have examined the effectiveness of metallic materials against SARS-CoV-2 leading to data scarcity for decision making. Generally, the long duration required for inactivation of coronaviruses by metals such as silver (Ag) and copper (Cu) indicates that they might be ineffective when rapid disinfection is required [73,74]. Furthermore, the high cost of these materials is also prohibitive when it comes to their mass adoption as an effective antiviral material. When it comes to 3D printed metallic materials, only Cu coating deposited using cold spray have reported 99.2% virus inactivation at 5 h [82]. Antiviral properties of Cu against SARS-CoV-2 are increasingly being documented at around 4-5 h [73,84]. Other than Ag and Cu based compositions, studies on the antiviral properties of alternative alloys against SARS-CoV-2 are yet to be reported. As such, the current state-of-art inactivation times for metallic alloys is around 4-5 h against SARS-CoV-2 which demands composition with superior antiviral efficacy where this study contributes.
According to Hatamie et al. [85], cobalt (Co) is an essential trace element that is significantly less toxic in comparison to non-essential metals. Studies on the antimicrobial properties of Co complexes [86] have shown effectiveness against seven microbial strains namely Bacillus subtilis, Escherichia coli, Klebsiella pneumoniae, Salmonella enterica, Shigella flexneri, Proteus vulgaris and Staphylococcus aureus. Recently, Kota et al. [87] showed Co-Cr to be ef-fective in suppressing S. aureus and P. acnes proliferation both in vitro and in vivo infection models. Overall, Co-based alloys seem to offer a broad-spectrum antimicrobial possibility making them worthy candidates to explore for potential antiviral activity. Currently, literature on the antiviral properties of Co-based alloys is scarce let alone their performance against SARS-CoV-2.
This research, therefore, is the first step towards drastically improving the state-ofthe-art antiviral alloys against SARS-CoV-2 at the interface of 3D printing and surrogate modelling. The research reveals the first Co-Cr-Mo porous material with superior antiviral activity that can be 3D printed on-demand where the innovation pipeline is kept open. The study is directed towards a process-structure-property relationship where both the material and its processing parameters collectively inform an optimum functional architecture. The influence of the LPBF 3D printing process on the structure-property relationship at the sub-micrometre is also analysed to make the digital and on-demand fabrication of the antiviral material accessible and easily scalable.
Overall, the development, analysis, and optimisation Co-Cr-Mo LPBF porous superalloy that feature high antiviral activity against an enveloped virus such as SARS-CoV-2 are demonstrated for the first time. The effect of LPBF process parameters on the characteristics of the porous architecture such as the thickness of the laser melted track (t t ) and pore diameter (φ d ) are also carried out. This was done with the help of a surrogate model that features laser power (P l ), hatch spacing (S h ) and scanning rate (V s ) as LPBF process parameters. The surrogate model was validated and subsequently used for parametric analysis which characterised the order of influence and interaction effects between the process parameters and the resulting printed Co-Cr-Mo architecture.

Morphology LPBF Samples
Although Co-Cr-Mo-based superalloys are suitable for LPBF [88][89][90], no studies have demonstrated a rationally conceived process informed by porosity at a track thicknesses below 300 µm. Before identifying the influence of the process parameters on the fabricated samples, the working limits of the LPBF process variables were assessed. This was a critical step as the porosity of Co-Cr-Mo is dictated by the process parameters as opposed to geometry.
The Co-Cr-Mo porous test specimens were printed for all LPBF parametric combinations informed by the rationally conceived surrogate model. The samples were removed from the build plate and analysed under SEM to characterise their porosity. The resulting morphologies of LPBF Co-Cr-Mo porous architecture informed by randomised parametric combinations are shown in Figure 1. The parametric combinations of laser power (P l ), scan speed (V s ), and hatch spacing (S h ) are sufficient for raising the required energy density at the powder bed for printing Co-Cr-Mo at varying porosity. The entire parametric range selected for the training matrix was found to be suitable to fully laser melt Co-Cr-Mo showing stable and continuous track formation.
As the pore size reduced significantly as shown in Figure 1j,k, some inconsistency in porosity and spatter contamination can be observed. This is due to the high energy density at the centre of the laser spot causing a recoil pressure at the melt pool while other parts of the melt pool are solidifying, which result in the expulsion of the molten material. The spattered metal subsequently cool down forming particles of varying sizes depending on the duration of the condensation process. These spatter formations are not unique to Co-Cr-Mo and are widely observed in a range of materials processed using powder bed fusion as summarised by Young et al. [91]. Although some spatters are observed, these are not extensive and can be seen to not obstruct the overall porous architecture being generated. Figure 1. SEM data of additively manufactured Co-Cr-Mo informed by randomised process parameter (S h is the hatch spacing, V s is the scanning rate and P l is the laser power) combinations that were subsequently used to train the surrogate model used to generate optimum Co-Cr-Mo microporous superalloy showing (a) S h = 0.60 mm, 60 mm, V s = 1025 mm/s, P l = 110 W, (j) S h = 0.14 mm, V s = 1025 mm/s, P l = 110 W, (k) S h = 0.14 mm, V s = 775 mm/s, P l = 110 W, (l) S h = 0.37 mm, V s = 775 mm/s, P l = 135 W.
Although the printed samples established the suitability of conceiving process-induced porosity, the ideal combination of the process parameters that will result in an optimum porous construct to maximise surface contact is required. However, the parameters should ensure that the generated φ d is not so low making the powder removal impossible. To establish such an optimum parametric combination, the order of influence of the process parameters on the resulting track thickness and pore size is required. This gives rise to a multi-objective optimisation problem which was established using a surrogate model.

Training Matrix and Regression Analysis
Controlling the three parameters for the formation of targeted track width, and pore size requires establishing a process-property relationship. The attempt is to fabricate Co-Cr-Mo architecture that features the thinnest track width and smallest possible distance between adjacent tracks while preserving porosity. Achieving this requires characterising both the interaction effects and order of influence of the LPBF process variables on the printed Co-Cr-Mo porous material.
To achieve the process-property relationship, a randomised BBD training matrix is generated as shown in Table 2 informed by the factors identified in Table 8. Keeping the layer thickness constant, three remaining primary LPBF process parameters that can vary the energy density at the powder bed were used as the variable factors. Consideration was also given when selecting the maximum and minimum limits of these process parameters to make sure sufficient energy density was available for Co-Cr-Mo melt pool generation.
Test prints were carried out in Co-Cr-Mo for all parametric combinations dictated by the matrix. Their characteristic results measured using SEM for track thickness and pore diameter are listed in Table 2. The results of the analysis were used to identify LPBF process parameters that had the most and least significance on the responses of interest. While the ideal combination of LPBF process parameters for fully dense Co-Cr-Mo is being increasingly documented [92][93][94], the literature is rather scarce when it comes to processing porosity and the fabrication of tracks with thickness below 300 µm. Table 2. Co-Cr-Mo surrogate model training matrix showing randomised parameters and the measured responses. S h is the hatch spacing, V s is the scanning rate, P l is the laser power, t t is the track thickness and φ d is the pore diameter.

Variable Factors
Responses (µm) The regression analysis of the training data from Table 2 based on best-fit indicators revealed that the thickness of Co-Cr-Mo track has a quadratic relationship with the LPBF process parameters as listed in Equation (1). The pore size or the linear distance between two adjacent tracks however was found to follow a linear relationship as listed in Equation (2). A quadratic dependency when it comes to track thickness indicates strong interaction effects between the LPBF process parameters considered. This means that specifying each of the LPBF process parameters require a critical understanding of the linking parameters and their interaction effects to accurately control the resulting t t .

Model Accuracy
The accuracy of the surrogate models was evaluated using the analysis of variance (ANOVA) technique where the significant model terms are as summarised in Table 3. The relevant accuracy indicators include the probability (p-value), coefficient of determination R 2 , Adjusted R 2 , and Adequate precision. It can be seen that all models exhibit high F-values and very low p-values confirming that the models are significant. In statistical terms, surrogate models with a p-value less than 0.05 and an adequate precision ratio greater than 4 signifies an accurate model [95]. A closer to unity R 2 and Adj-R 2 also indicates that the surrogate model is accurate for all the responses analysed. The relationship between actual responses and those based on the surrogate model are shown in Figure 2. As can be seen, the residuals are close to the predicted results which validate the accuracy of the surrogate model. Overall, the ANOVA demonstrates that all the models developed in this study are suitable for making valid predictions. This means that Equations (1) and (2) adequately characterise the relationship between the laser power, scanning rate, and hatch spacing and to that of the resulting properties of porous Co-Cr-Mo. This means that the surrogate model can be used to analyse the interaction effects between the process parameters and to identify their optimum combinations to print an optimum porous architecture.

Laser Power
When the laser scans over the powder layer, the absorbed energy from the laser beam heat the Co-Cr-Mo particles creating a melt pool. Heating, time-evolution of the melt pool, and the solidification process depend on the powder and process characteristics [96]. The process parameters influence the phases, recoil pressure, surface tension, Marangoni effect, and hydrodynamics that affects the size, and shape of the melt pool. When the laser beam leaves the melt area, the melt pool starts to cool down and solidifies [97,98]. As such the creation of a stable melt pool with a regular shape and geometrical characteristics, the powder characteristics of Co-Cr-Mo and laser processing parameters are the significant factors. These parameters also have a direct influence on the manufacturing time and quality of the final parts. Figure 3 shows the Co-Cr-Mo track thickness and pore diameter at varying laser power at a constant hatch spacing and scanning rate. As the laser power increases, the amount of interacting powder rises increasing the track thickness as shown in Figure 3a. At higher laser power, the Co-Cr-Mo track formation is contributed both from the particles directly and adjacent to the laser spot resulting in a denudation zone. The powder denudation zone defines the volume of powder involved in the track formation and spattering process. The process is not unique to Co-Cr-Mo and is similar to most metals when processed using LPBF [99].
Although there is a consistent increase in t t with laser power, the rate of change was found to reduce gradually post 115 W. This is because the rise in laser energy, melt convection, and thermal diffusion increases the unusable energy which cannot be completely absorbed in the powder bed. Consequently, the track thickness and laser power can be observed to have a significant influence but without an absolute proportionality throughout the laser power range considered; which is consistent with the literature [94,100,101]. Overall, a variation in laser power from 85 to 135 W was found to increase the track thickness by ∼26%. When it comes to the influence of laser power on the pore size (Figure 3b), a linear relationship was observed. It was found that the pore size was highest at 274 µm when the laser power was at its lowest signifying a thinner track. As the laser power increased to 135 W, the pore diameter reduced by 5% to 260 µm. As such, without considering any interaction between the process parameters, laser power has the highest influence on the track thickness, which subsequently affects the pore diameter to a lesser extent.

Scanning Rate
The scanning rate is the speed at which the laser spot travels across the powder bed. The first consideration when choosing the parameter is to achieve a consistent, fully dense track. For any given build, the powder characteristics and layer thickness are constant parameters. As such, using a fixed laser spot size leaves laser power, scanning rate, and hatch spacing as the process variables. Generally, all three of these parameters required being carefully controlled, for example, a high scanning rate can result in insufficient energy at the powder bed. This will lead to an unmelted or partially melted track featuring 'lack of fusion' porosity. In contrast, lowering the scanning rate leads to a high energy density that can overheat the melt pool, causing deeper energy penetration leading to keyhole formation and an unstable track. Therefore, there is an optimum scanning rate window that achieves a continuously fused material track. A fully dense material track is generally the aim, while the target of laser processing in this study, however, was to identify the thinnest most stable track.
The effect of scanning rate on the track thickness and the pore shape were analysed at a constant laser power and hatch spacing of 110 W and 0.37 mm respectively. The entire scanning rate range resulted in fully melted tracks. However, at a low scan rate, the laser interacts with the powder longer resulting in a thicker track. It can be seen from Figure 4a that the thickness of the molten Co-Cr-Mo track decreased as the scanning speed increased. However, the track thickness shows that the scanning rate has a threshold character characterised by the flattening of curves at a scanning rate around 775-875 mm/s (Figure 4a).
Although there is a consistent increase in t t with laser power, the rate of change was found to reduce gradually post 115 W. This is because the rise in laser energy, melt convection, and thermal diffusion increases the unusable energy which cannot be completely absorbed in the powder bed. Consequently, the track thickness and laser power can be observed to have a significant influence but without an absolute proportionality throughout the laser power range considered; which is consistent with the literature [94,100,101]. Overall, a variation in laser power from 85 to 135 W was found to increase the track thickness by ∼26%. Overall, the highest track thickness of 121 µm was observed at a scanning rate of 825 mm/s which reduced by 22% when the scanning rate was increased to 1025 mm/s. This is because as the scanning rate increases, the residence time of the laser per unit area reduces engaging in a lesser number of powdered particles leading to smaller melt-pool width. Nevertheless, Figure 4b shows that the changing track thickness has a small influence on the pore size where the diameter of the pore increased by ∼4% when the scanning rate was increased from 775 to 1025 mm/s. Similar to the case of laser power, the scanning rate was found to have a significant effect on the track thickness which subsequently affected the pore size to a lesser extent. As such, both the laser power and scanning rate did not directly affect the spacing between the melted tracks leading to a lower influence on pore size.

Hatch Spacing
The next level of complexity is the melting of multiple tracks to fill the designated part area based on the design. The part area is build-up of adjacent tracks, where their overlap distance is dictated by the hatch spacing. It is measured from the centre of one laser spot to the adjacent. To have process-induced porosity such as the ones targeted in this study, a large hatch spacing is required. Otherwise, the laser tracks will overlap resulting in a fully dense part. Figure 5 shows the influence of hatch spacing when the laser power and scanning rate are kept constant. It can be seen from Figure 5a that lower and higher hatch spacing seem to have no significant effect on the track thickness. This was expected as the hatch spacing range was deliberately chosen to not coincide or overlap with the previous track to induce a porosity. Nevertheless, a small rise in track width can be observed around a 0.37 mm hatch spacing which can be attributed to the effect of the powder particle size and some inconsistencies in the track formation. Only the analysis on the interaction effect will reveal the exact reason which is discussed in subsequent sections. From Figure 5b, it is evident that even at the smallest hatch spacing of 0.14 mm there exist a small gap between two consecutive laser tracks. This explains why the track thickness is unaffected despite the change in hatch spacing in Figure 5a. Evaluating Figure 5b, the hatch spacing has a significant effect on the pore size. The increase in pore diameter is linearly consistent with an increase in hatch spacing. The lowest hatch spacing was found to be 50 µm at a hatch spacing of 0.14 mm, however, this increased almost 10 folds when the hatch spacing was increased to 0.60. Therefore, when it comes to pore size, hatch spacing has the most significant effect despite keeping all other parameters constant. The track thickness on the other hand was primarily affected by the laser power and scanning rate with almost no influence from the hatch spacing range considered. Although the study so far has identified the influence of individual process parameters on the characteristics of the Co-Cr-Mo microporous material, the interaction effects are not understood, which is discussed in subsequent sections.

Track Thickness
Laser power, scanning rate, and hatch spacing are the three variables considered in this study. Although these parameters can be varied independently, their interaction has the most significant effect on the melt pool characteristics. Therefore, studying the interaction effects between the process parameters and identifying their order of influence is critical in identifying the optimum parametric combination. The thickness of Co-Cr-Mo track has a significant influence on the overall pore size and the surface area achievable. Ideally, a stable but thinner track offers a higher opportunity for functional porosity as opposed to thicker tracks. As such identifying and accurately characterising the interaction effects between the LPBF process parameters is carried out to inform the optimised microporous material.
The interaction effects between the process parameters affecting the Co-Cr-Mo track thickness are shown in Figure 6. As can be seen, the parameters behaved differently depending on their combination with other parameters being used. By observing the influence of scanning rate and hatch spacing as shown in Figure 6a, it can be seen that track thickness was slightly decreased as the scanning rate increased, resulting in the thinnest track at the highest scan speed with minimal influence from hatch spacing. As such, the interaction effects between scanning rate and hatch spacing on Co-Cr-Mo track thickness were insignificant, which is consistent with the single-parameter observations. This was expected as the hatch spacing range was deliberately chosen for porosity rather than to allow melt-pool overlap. A similar trend was observed regarding the interaction effects between laser power and hatch spacing, as shown in Figure 6b. Here, the thickness of the track was primarily driven by the laser energy with negligible interaction with hatch spacing. As such, t t increased consistently with the intensity of the laser power.
Overall, both Figure 6a,b show that hatch spacing has minimal influence when it comes to track thickness. As shown in Figure 6c, the Co-Cr-Mo track thickness is primarily dictated by P l and V s with significant interaction effects between the two parameters. Therefore, achieving the thinnest track required careful control of both scanning rate and laser power. The strong interaction effects mean that the thinnest track was formed when both the laser power and scanning rate were at their lowest. As such, controlling either P l or V s without considering their interaction is unlikely to achieve the thinnest possible LPBF process-induced Co-Cr-Mo track. Deriving the order of influence, the most significant terms on t t are the interaction effects of V s and P l in the order P l V s > V s > P l with the least influence from S h . Consequently, to generate the finest porous Co-Cr-Mo architecture, a higher scan speed and lower laser power that induces sufficient energy density for the thinnest but fully melted track thickness is required.

Pore Diameter
The pore diameter is primarily dictated by the hatch spacing, as evident from Figure 7a, with no influence from scanning rate or laser power, as shown in Figure 7b,c, respectively. The dependency of pore diameter on the hatch distance is also linear with the smallest and largest φ d , consistent with hatch spacing, as shown in Figure 7a,b. This response was expected as the parameter driving the distance between two adjacent tracks is hatch spacing. As shown in Figure 7a,b, variations in V s or P l cannot introduce any changes to pore diameter, as signified by identical response surfaces. It can be also seen that there are no interaction effects between V s and P l , as well demonstrated by the flat cure in Figure 7c. As such, the most significant term dictating the pore diameter of Co-Cr-Mo is the first-order effects of hatch spacing. The influence of process parameters on both track thickness and pore diameter reveals that the interaction effects can be used to simplify the Co-Cr-Mo microporous printing process. Overall, the targeted porosity can be achieved by carefully modulating the process parameters, which allows ease of customisability, and scalability. The influence of the process parameters on the material has also been analytically represented, which allows identifying the most optimal parametric combination based on the requirement. For the study under consideration a fine (<t t and < φ d ) architecture suitable for antiviral evaluation is the aim. The analysis so far shows that this can be achieved by minimising P l and S h while maximising V s which results in thinner tracks as close as possible but not so close as to introduce a dense part.

Predicting Optimal Solution
The analysis so far has demonstrated the influence of the process parameters in the development of Co-Cr-Mo porous material. However, an accurate parametric combination that will lead to the best possible porous architecture is yet to be established. Doing this requires a multi-objective numerical description of the optimisation problem. To maximise the antiviral load, the optimum Co-Cr-Mo material should be porous in a way that maximises the surface area while featuring stable and consistent tracks. In other words, the optimised structure should allow for the smallest φ d without producing a dense structure while featuring the thinnest t t . This means that the multi-objectives algorithm should look for a desirability criterion to minimise both φ d and t t to a non-zero positive integer. In this regard, the optimisation problem can be formulated as shown in Equation (3): A multi-objective optimisation problem usually yields multiple solutions. As such a desirability function is used to combine them into a single objective (D) that meets the desirable range for each response (d i ). The desirability approach is one of the most widely used methods in parametric optimisation where multiple responses are involved. It is based on the idea that the targeted outcome has multiple desired characteristics, any solution where the outcomes fall outside the desired limits is inadequate. Using this technique, the least and most desirable outcomes can be represented between 0 and 1 respectively. When n is the number of responses, the multi-objective function is a geometric mean of all transformed responses as shown in Equation (4): Figure 8 shows the results of the multi-objective optimisation as a desirability function considering all interaction effects informed by the surrogate model. The highest desirability of 0.9 can be achieved at a scanning rate of 975-1025 mm/s, laser power of 85-105 W, and a hatch spacing of 0.14 mm. Before identifying the exact numerical values, it was important to consider particle size, which in this case ranged from 3 to 18 µm. This was set to prevent the powder particles from getting stuck in the pores leading to a dense part. As can be seen from Figure 8a, despite keeping the other parameters constant increasing the laser power decreases the desirability. Overall, the lowest desirable solution is at the lowest scan speed, hatch distance, and highest laser power as shown in Figure 8c. Table 4 shows one of the parametric combinations that offer the highest desirability that was used to fabricate the Co-Cr-Mo validation sample.

Fabrication and Validation of the Optimised Co-Cr-Mo Architecture
Co-Cr-Mo samples were additively manufactured using the predicted LPBF process parameter combinations listed in Table 4. The printed sample was analysed under SEM which confirms a fine porous architecture as shown in Figure 9a. The magnified Figure 9b confirms an even and consistent track and pore formation. Any further reduction in pore diameter runs the risk of powder particles becoming embedded in the porosity. As such the predicted optimum architecture by the surrogate model is closer to the potential porosity limit that can be achieved using the current particle and laser spot size. The results from the experimentally measured optimum Cu-W-Ag microporous material when compared to the surrogate model can be seen to be in good agreement as listed in Table 5. Overall, the surrogate model overestimated the track thickness and pore diameter by 3.9 and 3.2% respectively. This small change between the predicted and experimented data can be attributed to the influence of particle size and partially melted powder at the track boundaries. Overall, the results show that the Co-Cr-Mo porous architecture with a feature size of 75 µm and pore diameter of 61 µm is achievable using LPBF, which transcends the current state-of-the-art feature size of 300 µm. As such, any further refinement in achievable pore size requires finer feedstock and further adjustments in laser modulation.
Subsequently, the as-built Co-Cr-Mo optimised porous material was used for antiviral testing against the SARS-CoV-2 viral model to evaluate its antiviral effectiveness. Table 5. Comparison between actual and predicted values of the optimum Co-Cr-Mo microporous material on track thickness (t t ) and pore diameter (φ d ).

Antiviral Characterisation
The optimised microporous Co-Cr-Mo material was printed, and antiviral tests were performed (n = 3) against the biosafe vrial model of SARS-CoV-2. The antiviral performance is summarised in Table 6, where the data at 30 min and 5 h compared against a control sample, commercial porous fabric, and current state-of-the-art copper architecture. The corresponding visuals of the samples at 30 min and 5 h of contact with the phage are presented in Figure 10. The control sample did not exhibit any antiviral activity at 30 min or 5 h as shown in Figure 10a,d, respectively. The commercial porous fabric, which was used as a reference material showed no antiviral performance at 30 min (Figure 10b) or 5 h (Figure 10e) as expected. In comparison, the optimised Co-Cr-Mo LPBF porous architecture resulted in 100% viral inactivation at both contact times where no plaques can be observed as shown in Figure 10c,f, respectively. The results revealed that the novel Co-Cr-Mo microporous architecture developed in this study have excellent antiviral properties (100% viral inactivation in 30 min) against phage phi 6 used as enveloped RNA viral model of SARS-CoV-2. The antiviral performance of the Co-Cr-Mo porous material is also superior to the current state-of-the-art copper [82], which showed 99.2% inactivation of SARS-CoV-2 in 5 h and no antiviral performance at 30 min as compared in Table 6. Generally, the SARS-CoV-2 virus can survive for 4-5 days in a range of materials [105] such as plastics, ceramics, stainless steel and glass as listed in Table 1. As such this paper presents the first metallic material that shows high (<30 min) antiviral performance against an enveloped virus like SARS-CoV-2 and influenza virus. This means that the material is suitable to be used as permanent antiviral devices aiding response against viral pathogens. Further development in this direction can aid pandemic preparedness to protect human beings in a more effective way while reducing the environmental impacts of disposable viral control devices in critical environments.

Laser Powder Bed Fusion (LPBF)
The EOS M290 LPBF machine was used to fabricate the samples at a constant layer thickness of 30 µm. The system featured a 400 W Yb-fibre laser that is modulated above the powder bed within a 250 × 250 mm build platform. Table 7 shows the composition (wt. %) of the atomised Co-Cr-Mo alloy used as feedstock at the powder bed. The LPBF process parameters namely the laser power (P l ), hatch spacing (S h ) and scanning rate (V s ) was informed by the surrogate model which is discussed in subsequent sections. The bulk density of the material is 8.3 g/cm 3 with the morphology of the particles as shown in Figure 1. The particles featured a spherical morphology with occasional irregular shapes representative of typical feedstock suitable for LPBF. Some smaller satellite particles around 3 µm can be seen attached to larger particles around 18 µm. The overall sphericity and particle size range was found to be favourable for the powder bed fusion process and resulted in even powder spread. The laser processing was carried out in an argon inert atmosphere at <0.1% oxygen on a steel substrate of temperature 35 • C. The number of samples was informed by the required parametric combinations randomised by the surrogate model. The LPBF process prints follow a layer-by-layer method where the build platform is lowered by the constant layer thickness between subsequent layers. Post-printing the samples were heat-treated at 1050 • C for six hours in an argon atmosphere. The resulting track thickness and pore dimensions were characterised using scanning electron microscopy (SEM). EVO 50 SEM produced by Zeiss that uses an incident electron beam to interact with the printed sample to generate backscattered and secondary electrons to create an image of the porous sample is used (Figure 11) [48][49][50][51][52]. Figure 11. Scanning electron micrograph of Co-Cr-Mo feedstock showing particle morphology.

Surrogate Modelling
In LPBF, the extent of material melting and track formation is dictated by the energy density (e ld ) [106] at the powder bed which is given by Equation (5): where P l is the laser power (W), V s is the scanning rate (mm/s), S h is the hatch spacing (mm), and t l is the layer thickness. The relationship between the LPBF process parameters and the energy density at the powder bed means that the parameters can be controlled for the targeted outcome. Traditionally, this approach is used to identify the optimal energy density required to fully melt and fuse the feedstock, resulting in a stable track. Building on this relationship, this study attempted to identify the required LPBF parametric combinations to create a Co-Cr-Mo construct with controllable pore size and track thickness. To carry this out, the relationship between the LPBF process parameters is linked to an objective function of the resulting architecture using Equation (6): where x = (x 1 , x 2 , . . . , x k ) is the vector of the k LPBF process variables for Co-Cr-Mo. The maximum and minimum limits associated with each process parameter are defined by x l and x u for the objective function f (x). Printing porous Co-Cr-Mo requires characterising the order of influence of the LPBF variables that leads to specific responses. The methodology was to vary the laser energy at the powder bed by algorithmically modifying the LPBF parameters P l , S h , and V s to inform the track thickness (t t ) and pore diameter (φ d ), as shown in Figure 12.
The parametric combination of the randomised training matrix for the surrogate model was conceived based on the Box-Behnken design (BBD) principle. BBD was chosen as it allows for higher-order response surfaces using fewer required runs than a normal factorial technique in the training matrix. As such it is an alternative to the 3 k factorial that results in an efficient sampling matrix of the coded factorial levels. This approach generally results in a good fit with sufficient information to test for lack of fit. The method also allows models of increasing order to be constructed sequentially and allows for an estimation of experimental error. Note: S h is the hatch spacing, V s is the scanning rate, P l is the laser power, t t is the track thickness and φ d is the pore diameter.
The limits of the LPBF process parameters considered for the coded BBD factorial matrix are summarised in Table 8. Based on the problem description, regression analysis was used to characterise the relationship between the LPBF process variables and the resulting responses of the Co-Cr-Mo printed samples. Subsequently, best fit empirical models were derived through randomised experimental data measured for the responses t t and φ d . The resulting polynomial functions were used to predict the order of influence of the contributing process parameters on the responses of the Co-Cr-Mo printed samples. The steps that led to the development of an accurate surrogate model is summarised in Figure 13. The surrogate model was trained using randomised test data to satisfy the sampling matrix and analysis of variance criteria. The predictions of the trained model were evaluated for accuracy followed using a desirability criterion to identify the optimum parametric combination that leads to the thinnest Co-Cr-Mo melt tracks and smallest pore diameter. This optimum criterion was chosen to minimise pore diameter and track thickness to accommodate the highest number of pores as possible. The surrogate model was subsequently used to quantify the interaction effects of the process parameters and their influence on the Co-Cr-Mo porous architecture.

Antiviral Characterisation
The phage phi6 host culture was carried using Pseudomonas syringae from the DSMZ-German collection of microorganisms and cell cultures (DSM 21482) in solid Tryptic Soy Agar (TSA, Liofilchem). Post culture, the bacteria were incubated in liquid Tryptic Soy Broth (TSB, Liofilchem) at 120 rpm and at 25 • C. Pseudomonas phage phi6 titration was carried out in accordance with specifications for phage infection.
Phi6 is a double-stranded RNA virus with three-part, segmented, totalling~13.5 kb in length. This type of lytic bacteriophage belongs to Group III Baltimore classification and was chosen as SARS-CoV-2 viral model (Group IV: positive-sense single-stranded RNA viruses) as it features a lipid membrane around their nucleocapsid. According to Kitajima et al. [107] low pathogenic CoV strains (such as MHV and classical human CoVs) or the same enveloped Pseudomonas phage phi6 can be used as models of SARS-CoV-2 considering biosafety [108]. Furthermore, recent experiments performed with Phi6 and SARS-CoV-2 have validated the use of this biosafe viral model [26]. Figure 13. The surrogate modelling methodology for LPBF Co-Cr-Mo porous material. Note: S h is the hatch spacing, V s is the scanning rate, P l is the laser power, t t is the track thickness and φ d is the pore diameter.
LPBF processed Co-Cr-Mo porous samples (n = 3) of 10 mm dia. were characterised. A non-woven spunlace fabric filter from NV EVOLUTIA (commercial filters used for face masks) of 10 mm in diameter was cut with a cylindrical punch and used as reference material. All samples were dried at 60 • C for 24 h between glass plates and sterilised in ethanol/distilled water (70:30 vol.%) solution for 5 min at 25 • C and UV radiated for an hour on each side.
The antiviral activity of the samples was measured at 30 min and five hours of contact with the biosafe viral model. A 50 µL volume of phage suspension in TSB was introduced into each sample at a concentration of 1 × 10 6 plaque-forming units per mL (PFU/mL) and incubated for 30 min and 5 h, respectively. Then the samples were placed in a falcon tube with 10 mL TSB and sonicated for 5 min at 24 • C and vortexed for 1 min. Serial dilutions of each falcon were made and 100 µL of each phage dilution were placed in contact with 100 µL of the host strain at OD 600mm = 0.5. The infective capacity of the phage was measured by the double-layer method where 4 mL of top agar (TSB + 0.75% bacteriological agar, Scharlau) and 5mM CaCl 2 were introduced to the phage-bacteria mix. The mixture was poured on TSA plates that were incubated for 24-48 h in an oven at 25 • C.

Conclusions
This study demonstrates the use of Laser Powder Bed Fusion for on-demand manufacturing of a novel Co-Cr-Mo porous material that shows superior antiviral activity outperforming common antiviral metals such as silver and copper. The paper reveals a surrogate model that allows porosity personalisation without the need for complex geometry data. This is achieved through controlling the LPBF process parameters such as laser power, scanning rate and hatch spacing. The proposed methodology simplifies the data requirement and pre-processing often required for printing porous materials. The surrogate model developed in this study showed that the most significant parameters for Co-Cr-Mo track thickness (t t ) were the interaction effects of scanning rate (V s ) and laser power (P l ) in the order P l V s > V s > P l . For pore diameter (φ d ), the hatch spacing (S h ) has the most significant effect. The optimised Co-Cr-Mo microporous materials showed 100% viral inactivation against the phage phi 6 used as enveloped RNA viral model of SARS-CoV-2 in 30 min. The evolution of this and future pandemics will bring unexpected situations where the ability to print and personalise on-demand antiviral materials can achieve rapid solutions. Furthermore, the proposed methodology can be adopted to conceive functional antiviral materials that can be fabricated close to point-of-care.