Engineering biomaterials by inkjet printing of hydrogels with functional particulates.

Hydrogels with particulates, including proteins, drugs, nanoparticles, and cells, enable the development of new and innovative biomaterials. Precise control of the spatial distribution of these particulates is crucial to produce advanced biomaterials. Thus, there is a high demand for manufacturing methods for particle-laden hydrogels. In this context, 3D printing of hydrogels is emerging as a promising method to create numerous innovative biomaterials. Among the 3D printing methods, inkjet printing, so-called drop-on-demand (DOD) printing, stands out for its ability to construct biomaterials with superior spatial resolutions. However, its printing processes are still designed by trial and error due to a limited understanding of the ink behavior during the printing processes. This review discusses the current understanding of transport processes and hydrogel behaviors during inkjet printing for particulate-laden hydrogels. Specifically, we review the transport processes of water and particulates within hydrogel during ink formulation, jetting, and curing. Additionally, we examine current inkjet printing applications in fabricating engineered tissues, drug delivery devices, and advanced bioelectronics components. Finally, the challenges and opportunities for next-generation inkjet printing are also discussed.

On the efficacy of facial masks to suppress the spreading of pathogen-carrying saliva particles during human respiratory events: Insights gained via high-fidelity numerical modeling.

Respiratory fluid dynamics is integral to comprehending the transmission of infectious diseases and the effectiveness of interventions such as face masks and social distancing. In this research, we present our recent studies that investigate respiratory particle transport via high-fidelity large eddy simulation coupled with the Lagrangian particle tracking method. Based on our numerical simulation results for human respiratory events with and without face masks, we demonstrate that facial masks could significantly suppress particle spreading. The studied respiratory events include coughing and normal breathing through mouth and nose. Using the Lagrangian particle tracking simulation results, we elucidated the transport pathways of saliva particles during inhalation and exhalation of breathing cycles, contributing to our understanding of respiratory physiology and potential disease transmission routes. Our findings underscore the importance of respiratory fluid dynamics research in informing public health strategies to reduce the spread of respiratory infections. Combining advanced mathematical modeling techniques with experimental data will help future research on airborne disease transmission dynamics and the effectiveness of preventive measures such as face masks.

Bark and biochar in horizontal flow filters effectively remove microplastics from stormwater.

Organic materials such as bark and biochar can be effective filter materials to treat stormwater. However, the efficiency of such filters in retaining microplastics (MPs) - an emerging stormwater pollutant - has not been sufficiently studied. This study investigated the removal and transport of a mixture of MPs commonly associated with stormwater. Different MP types (polyamide, polyethylene, polypropylene, and polystyrene) were mixed into the initial 2 cm material of horizontal bark and biochar filters of 25, 50, and 100 cm lengths. The MP types consisted of spherical and fragmented shapes in size ranges of 25-900 µm. The filters were subjected to a water flow of 5 mL/min for one week, and the total effluents were analyzed for MPs by µFTIR imaging. To gain a deeper insight, one 100 cm bark filter replica was split into 10 cm segments, and MPs in each segment were extracted and counted. The results showed that MPs were retained effectively, >97%, in all biochar and bark filters. However, MPs were detected in all effluents regardless of filter length. Effluent concentrations of 5-750 MP/L and 35-355 MP/L were measured in bark and biochar effluents, respectively, with >91% of the MP counts consisting of small-sized (25 µm) polyamide spherical particles. Combining all data, a decrease in average MP concentration was noticed with longer filters, likely attributed to channeling in a 25 and 50-cm filter. The analyses of MPs in the bark media revealed that most MPs were retained in the 0-10 cm segment but that some MPs were transported further, with 19% of polyamide retained in the 80-90 cm segment. Overall, this study shows promising results for bark and biochar filters to retain MPs, while highlighting the importance of systematic packing of filters to reduce MP emissions to the environment from polluted stormwater.

Stray Magnetic Field Variations and Micromagnetic Simulations: Models for Ni0.8Fe0.2 Disks Used for Microparticle Trapping.

Patterned micro-scale thin-film magnetic structures, in conjunction with weak (~few tens of Oe) applied magnetic fields, can create energy landscapes capable of trapping and transporting fluid-borne magnetic microparticles. These energy landscapes arise from magnetic field magnitude variations that arise in the vicinity of the magnetic structures. In this study, we examine means of calculating magnetic fields in the local vicinity of permalloy (Ni0.8Fe0.2) microdisks in weak (~tens of Oe) external magnetic fields. To do this, we employ micromagnetic simulations and the resulting calculations of fields. Because field calculation from micromagnetic simulations is computationally time-intensive, we discuss a method for fitting simulated results to improve calculation speed. Resulting stray fields vary dramatically based on variations in micromagnetic simulations-vortex vs. non-vortex micromagnetic results-which can each appear despite identical simulation final conditions, resulting in field strengths that differ by about a factor of two.

Coalescence-Induced Self-Propelled Particle Transport with Asymmetry Arrangement.

We experimentally investigated the coalescence-induced droplet-particle jumping phenomenon on a submillimeter scale in symmetric and asymmetric particle arrangements with poly(methyl methacrylate) (PMMA) particles and stainless steel (SS) particles. Coalescence-induced droplet-particle jumping exhibited excellent capability and interesting behavior for both droplet jumping enhancement and particle transport. The particle increased the normalized droplet jumping velocity from 0.250 for no particle case to 0.315 and 0.320 for symmetric and asymmetric particle cases. Compared with similar-sized macrostructures fixed between droplets, better jumping performance with particles may be attributed to avoiding the work of adhesion during droplet-macrostructure separation. Besides, all particles always sunk at the bottom in the symmetric cases, while the stick mode for PMMA particles and sink, wander, and jet modes for SS particles appeared in the asymmetry cases. We revealed that the asymmetric particle arrangement induces an unbalanced surface tension force, which may provide a driving force in the vertical direction. Additionally, a small enough resistive force caused by hydrophobic particles is another necessary condition for the wonder and jet mode. Finally, we realized a significantly superior particle transport in the asymmetric SS particle cases with maximum particle height reaching ∼2.1 mm, ∼12.4 times the particle radius, the most significant vertical self-propelled transport distance currently.

Parameter characteristics in intranasal drug delivery: A key to targeting medications to the olfactory airspace.

BACKGROUND: The nose is a viable pathway for topical drug delivery to the olfactory cleft for treatment of obstructive smell loss and nose-to-brain drug delivery. This study investigates how variations in nasal vestibule morphology influence intranasal spray drug transport to the olfactory cleft and olfactory roof/bulb regions. METHODS: The unilateral nasal vestibule morphology in three healthy subjects with healthy normal nasal anatomy was classified as Elongated (Subject DN001), Notched (Subject DN002), and Standard (Subject DN003). Computational fluid and particle dynamics modelling were used to simulate nasal airflow and drug particle transport to the olfactory cleft and olfactory roof/bulb regions in each subject-specific nasal cavity. To evaluate highest drug depositions in these regions, the particle transport simulations involved extensive parameter combination analyses: 6 inspiratory flow rates mimicking resting to sniffing (10-50 L/min); 5 spray release locations (Top, Bottom, Central, Lateral, and Medial); 5 head positions (Upright, Tilted Forward, Tilted Back, Supine, and Mygind); 3 particle velocities (1, 5, and 10 m/s); 350,000 µm-particles (1-100 µm) and 346,500 nanoparticles (10-990 nm). FINDINGS: Particle size groups with highest depositions in olfactory cleft: DN001 left = 28.4% at 11-20 µm, right = 75.3% at 6-10 µm; DN002 left = 16.8% at 1-5 µm, right = 45.3% at 30-40 nm; DN003 left = 29.1% at 21-30 µm, right = 15.9% at 6-10 µm. Highest depositions in olfactory roof/bulb: DN001 left = 6.5% at 11-20 µm, right = 26.4% at 11-20 µm; DN002 left = 3.6% at 1-5 µm, right = 2.6% at 1-5 µm; DN003 left = 2.8% at 21-30 µm, right = 1.7% at 31-40 µm. INTERPRETATION: DN001 (Elongated nasal vestibule) had the most deposition in the olfactory regions. Micron-particles size groups generally had better deposition in the olfactory regions.

Analysis of Lanthanum Oxide Based Double-Gate SOI MOSFET using Monte-Carlo Process.

INTRODUCTION: This work proposes a Double-Gate (DG) MOSFET with a Single Material made of Silicon On-Insulator (SOI). The Lanthanum Oxide material with a high k-dielectric constant has been used as an interface between two gates and the channel. The Monte Carlo analysis has been used to determine the Conduction Band Energy (Ec) profiles and electron sheet carrier densities (ns) for a Silicon channel thickness (tsi) of 10 nm at 0.5 V gate drain-source voltages. The transverse electric fields are weak at the midchannel of DG SOI MOSFETs, where quantum effects are encountered. The Monte Carlo simulation has been confirmed to be effective for high-energy transport. A particle description reproduces the granularity property of the transport for nanoscale modeling. METHODS: This work utilizes a Monte Carlo (MC) Simulation for the proposed Double Gate Single Material Silicon On Insulator MOSFET with (La2O3=2 nm) as dielectric oxide on upper and lower gate material. The electrical properties of the DG SOI MOSFETs with Lanthanum Oxide were analyzed using Monte Carlo simulation, including the conduction band energy, electric field, potential distribution, particle movement, and average velocity. RESULTS: The peak electric field (E) simulation results and an average drift velocity (υavg) of 6Í105 V/cm and 1.6Í107 cm/s were obtained, respectively. The conduction band energy for the operating region of the source has been observed to be 4 % to the drain side, which obtained a value of -0.04 eV at the terminal end. CONCLUSION: This proposed patent design, such as double-gate SOI-based devices, is the best suggestion for significant scalability challenges. Emerging technologies reach the typical DG SOI MOSFET's threshold performance when their geometrical dimensions are in the nanometer region. This device based on nanomaterial compounds has been more submissive than conventional devices. The nanomaterials usage in the design is more suitable for downscaling and reducing packaging density.

Combined Funnel, Concentrator, and Particle Valve Functional Element for Magnetophoretic Bead Transport Based on Engineered Magnetic Domain Patterns.

Controlled actuation of superparamagnetic beads (SPBs) within a microfluidic environment using tailored dynamic magnetic field landscapes (MFLs) is a potent approach for the realization of point-of-care diagnostics within Lab-on-a-chip (LOC) systems. Making use of an engineered magnetic domain pattern as the MFL source, a functional LOC-element with combined magnetophoretic "funnel", concentrator, and "valve" functions for micron-sized SPBs is presented. A parallel-stripe domain pattern design with periodically decreasing/increasing stripe lengths is fabricated in a topographically flat continuous exchange biased (EB) thin film system by ion bombardment induced magnetic patterning (IBMP). It is demonstrated that, upon application of external magnetic field pulses, a fully reversible concentration of SPBs at the domain pattern's focal point occurs. In addition, it is shown that this functionality may be used as an SPB "funnel", allowing only a maximum number of particles to pass through the focal point. Adjusting the pulse time length, the focal point can be clogged up for incoming SPBs, resembling an on/off switchable particle "valve". The observations are supported by quantitative theoretical force considerations.

Microplastic retention in marine vegetation canopies under breaking irregular waves.

The present study provides indications and underlying drivers of wave-induced transport and retention potential of microplastic particles (MP) in marine vegetation canopies having different densities. The anthropogenic occurrence of MP in coastal waters is well documented in the recent literature. It is acknowledged that coastal vegetation can serve as a sink for MP due to its energy dissipating features, which can mimic a novel ecosystem service. While the transport behavior of MP in vegetation has previously been investigated to some extent for stationary flow conditions, fundamental investigations for unsteady surf zone flow conditions under irregular waves are still lacking. Herein, we demonstrate by means of hydraulic model tests that a vegetation's retention potential of MP in waves increases with the vegetation shoot density, the MP settling velocity and decreasing wave energy. It is found that particles migrating by traction (predominantly in contact with the bed) are trapped in the wake regions around a canopy, whereas suspended particles are able to pass vegetated areas more easily. Very dense canopies can also promote the passage of MP with diameters larger than the plant spacing, as the canopies then show characteristics of a solid sill and avoid particle penetration. The particle migration ability through a marine vegetation canopy is quantified, and the key drivers are described by an empirical expression based on the particle settling velocity, the canopy length and density. The findings of this study may contribute to improved prediction and assessment of MP accumulation hotspots in vegetated coastal areas and, thus, may help in tracing MP sinks. Such knowledge can be considered a prerequisite to develope methods or new technologies to recover plastic pollutants and rehabilitate valuable coastal environments.

Geomorphological changes in young soils with sparse vegetation: Mathematical modeling and numerical simulation.

This article presents an improved mathematical model and numerical simulation for weathering of large areas with complex topography. It uses the equations of momentum, temperature, and humidity in turbulent air and for heat and water infiltration into soils. A mathematical model is also presented to calculate the soil porosity fraction produced by physical rock weathering in areas where soil is produced from intrusive rocks (batholiths). An algorithm based on air velocity, humidity (rainfall), temperature variation, and soil topography was developed to quantify soil erosion and change of relief at each point and time step in air, at the ground surface, and within the soil. This results in a complete air-soil model based on conservation laws that have not previously been applied to large areas of the earth's surface. The mathematical model is solved using large-scale numerical simulations applied to an area of 6.6 km2 in the Sierra Nevada batholith of California, USA. The results show that the wind velocity and resulting erosion is greater in areas with steeper slopes and that moisture accumulates mainly in low and flat areas; therefore, erosion is not uniform throughout the study area. In addition, computer simulations localized calculations to discrete grid cells within the porous (saprolite) fraction of the soil produced by freezing and thawing of water in rock. Results indicate that this physical mechanism is the primary contributor to weathering of rock at the study area.

Evolution Mechanism of Sputtered Film Uniformity with the Erosion Groove Size: Integrated Simulation and Experiment.

In this work, Cu thin films were experimentally fabricated at different target-substrate distances by 2-inch and 4-inch circular planar magnetron targets. Meanwhile, the sputtering deposition of Cu thin films was investigated via an integrated multiscale simulation, where the magnetron sputtering discharge was modeled using the Monte Carlo (MC) method, and the sputtered particle transport was simulated using a coupled Monte Carlo (MC) and molecular dynamics (MD) method. Experimental results indicated that, as the target-substrate distance increased from 30 to 120 mm, the film thickness distribution of the 2-inch target sputtering changed from a bell-shaped curve to a line-shaped curve, while that of the 4-inch target sputtering varied from a saddle-shaped curve to a line-shaped curve. The simulation results were accordant with the experimental results. The simulation results revealed that, at a target-substrate distance of 30 mm, the sputtering particle flow from the 2-inch target overlapped strongly near the substrate center, leading to a bell-shaped film thickness distribution, while the increased diameter of the erosion groove on the 4-inch target reduced the superposition effect of the sputtering particle flow near the substrate center, resulting in a saddle-shaped film thickness distribution. In addition, when the target-substrate distance ranged from 30 to 120 mm, the film thickness uniformity of 4-inch target sputtering was superior to that of 2-inch target sputtering, and the underlying mechanism was discussed in detail.

Vertical connectivity of microbiome and metabolome reveals depth-dependent variations across a deep cold-seep water column.

Deciphering the vertical connectivity of oceanic microbiome and metabolome is crucial for understanding the carbon sequestration and achieving the carbon neutrality. However, we lack a systematic view of the interplay among particle transport, microbial community, and metabolic trait across depths. Through integrating the biogeochemical, microbial, and metabolic characteristics of a deep cold-seep water column (∼1989 m), we find the altered connectivity of microbial community and dissolved organic matter (DOM) across depths. Both the microbial communities (bacteria and protists) and DOM show a clear compositional connectivity from surface to the depth of 1000 m, highlighting the controls of sinking particle over microbial connectivity from the epipelagic to mesopelagic zone. However, due to the biological migration and ocean mixing, the fecal-associated bacteria and protistan consumers unexpectedly emerge and the degradation index of DOM substantially alters around 1000-1200 m. Collectively, we unveil the significance of multi-faceted particle dispersion, which supports the connectivity and variability of deep ocean microbial communities.

Hydrodynamics and effect of velocity on particle filtration due to bridging in water-saturated porous media using CFD-DEM simulation.

Particle bridging owing to the confinement of the pore structure affects the transport and retention of particles in porous media. Particle motion driven by gravities were well investigated, whose filtration is mainly affected by the ratio of the particle diameter to the pore throat size of the medium. However, particles whose motions are driven by the fluid is essential to be investigated for particle separation from the carrying fluid. In this study, the motion of particles was driven by the liquid when passing through a water-saturated porous medium. The fluid-particle flow in a porous medium was modeled using computational fluid dynamics-discrete element method. The motion of particles in the slurry was traced in the porous medium, which enabled particle clogging to be directly precited by the interaction between the particles and pore surfaces by assessing the exact location of each particle. The pressure and flow field of the liquid were investigated, and the variation in flow path owing to particle clogging was predicted. The hydrodynamic study also showed that the Stokes number and particle concentration determined the particle clogging at the pore throats of the porous medium. Increasing the fluid velocity of particles such that the Stokes number was almost equal to 1 increased the separation efficiency of particles. Further increasing the fluid velocity reduced the residence time, which reduced the separation efficiency of the particles.

Bayesian inference of polymerase dynamics over the exclusion process.

Transcription is a complex phenomenon that permits the conversion of genetic information into phenotype by means of an enzyme called RNA polymerase, which erratically moves along and scans the DNA template. We perform Bayesian inference over a paradigmatic mechanistic model of non-equilibrium statistical physics, i.e. the asymmetric exclusion processes in the hydrodynamic limit, assuming a Gaussian process prior for the polymerase progression rate as a latent variable. Our framework allows us to infer the speed of polymerases during transcription given their spatial distribution, while avoiding the explicit inversion of the system's dynamics. The results, which show processing rates strongly varying with genomic position and minor role of traffic-like congestion, may have strong implications for the understanding of gene expression.

Effects of tumour heterogeneous properties on modelling the transport of radiative particles.

Dose calculation plays a critical role in radiotherapy (RT) treatment planning, and there is a growing need to develop accurate dose deposition models that incorporate heterogeneous tumour properties. Deterministic models have demonstrated their capability in this regard, making them the focus of recent treatment planning studies as they serve as a basis for simplified models in RT treatment planning. In this study, we present a simplified deterministic model for photon transport based on the Boltzmann transport equation (BTE) as a proof-of-concept to illustrate the impact of heterogeneous tumour properties on RT treatment planning. We employ the finite element method (FEM) to simulate the photon flux and dose deposition in real cases of diffuse intrinsic pontine glioma (DIPG) and neuroblastoma (NB) tumours. Importantly, in light of the availability of pipelines capable of extracting tumour properties from magnetic resonance imaging (MRI) data, we highlight the significance of such data. Specifically, we utilise cellularity data extracted from DIPG and NB MRI images to demonstrate the importance of heterogeneity in dose calculation. Our model simplifies the process of simulating a RT treatment system and can serve as a useful starting point for further research. To simulate a full RT treatment system, one would need a comprehensive model that couples the transport of electrons and photons.

Diffusiophoretic Particle Penetration into Bacterial Biofilms.

Bacterial biofilms are communities of cells adhered to surfaces. These communities represent a predominant form of bacterial life on Earth. A defining feature of a biofilm is the three-dimensional extracellular polymer matrix that protects resident cells by acting as a mechanical barrier to the penetration of chemicals, such as antimicrobials. Beyond being recalcitrant to antibiotic treatment, biofilms are notoriously difficult to remove from surfaces. A promising, but relatively underexplored, approach to biofilm control is to disrupt the extracellular polymer matrix by enabling penetration of particles to increase the susceptibility of biofilms to antimicrobials. In this work, we investigate externally imposed chemical gradients as a mechanism to transport polystyrene particles into bacterial biofilms. We show that preconditioning the biofilm with a prewash step using deionized (DI) water is essential for altering the biofilm so it takes up the micro- and nanoparticles by the application of a further chemical gradient created by an electrolyte. Using different particles and chemicals, we document the transport behavior that leads to particle motion into the biofilm and its further reversal out of the biofilm. Our results demonstrate the importance of chemical gradients in disrupting the biofilm matrix and regulating particle transport in crowded macromolecular environments, and suggest potential applications of particle transport and delivery in other physiological systems.

A numerical study on microfluidic devices to maintain the concentration and purity of dielectrophoresis-induced separated fractions of analyte.

Determining the physical and chemical properties of biologically important particles such as cells, organelles, viruses, exosomes, complexes, nucleotides, and proteins is needed to understand their function. These properties are determined with common analytical tools (mass spectrometry, cryo-EM, NMR, various spectroscopies, nucleotide sequencing, etc.) whose function can be improved when samples are pure and concentrated. Separations science plays a central role in conditioning samples, ranging from low-resolution benchtop operations like precipitations or extractions to higher-resolution chromatography and electrophoresis. In the last two decades, gradient insulator-based dielectrophoresis (g-iDEP) has emerged as a high-resolution separation technique capable of highly selective enrichment of cells, viruses, exosomes, and proteins. Specific evidence has been shown that pure homogeneous and concentrated fractions of cells and exosomes can be generated from complex mixtures. However, recovering those fractions for analysis has not been developed, limiting the technique to an analytical rather than a preparative one. Here, a finite element analysis was undertaken to identify geometries and operational parameters to efficiently remove the enriched fraction while retaining maximum concentration and providing total mass transfer. Geometric factors (e.g., side channel width and distance from the gradient-inducing gap) were studied, along with the addition of a second inlet side channel. Two flow-generating mechanisms-electroosmosis and hydrostatic pressure-were evaluated for semi-optimized device designs, including a comparison of the one- and two-inlet designs. Simulations indicate effectively one hundred percent mass transfer and a concentration increase by an order of magnitude for several device configurations and operational parameters.

Holographic generation of arbitrary ultrasonic fields by simultaneous modulation of amplitude and phase.

Acoustic holograms have been used widely to generate desired acoustic fields. Following the rapid development of 3D printing technology, the use of holographic lenses has become an efficient method to produce acoustic fields with high resolution and low cost. In this paper, we demonstrate a technique to modulate the amplitude and phase of ultrasonic waves simultaneously using a holographic method with high transmission efficiency and high accuracy. On this basis, we generate an Airy beam with high propagation invariance. We then discuss the advantages and disadvantages of the proposed method when compared with the conventional acoustic holographic method. Finally, we design a sinusoidal curve with a phase gradient and a constant pressure amplitude and realize transport of a particle on a water surface along a curve.

Fates of Emitted Particles Depending on Mask Wearing Using an Approach Validated Across Spatial Scales.

The spread of emitted potentially virus-laden aerosol particles is known to be highly dependent on whether a mask is worn by an infected person and on the emission scenario, i.e., whether the person is coughing, speaking, or breathing. The aim of this work is to investigate in detail the fates of particles emitted by a person wearing a perfectly fitting, a naturally fitted mask with leakage, and no mask depending on the emission scenario. Therefore, a two-scale numerical workflow is proposed where parameters are carried through from a micro-scale where the fibers of the mask filter medium and the aerosol particles are resolved to a macro-scale and validated by comparison to experimental measurements of fractional filtration efficiency and pressure drop of the filter medium as well as pressure drop of the mask. It turns out that masks reduce the number of both emitted and inhaled particles significantly even with leakage. While without a mask, the person opposite of an infected person is generally at the highest risk of being infected, a mask worn by an infected person speaking or coughing will deflect the flow leading to the fact that the person behind the infected person might inhale the largest number of aerosol particles.

On the Microcrystal Structure of Sputtered Cu Films Deposited on Si(100) Surfaces: Experiment and Integrated Multiscale Simulation.

Sputtered Cu/Si thin films were experimentally prepared at different sputtering pressures and characterized using X-ray diffraction (XRD) and an atomic force microscope (AFM). Simultaneously, an application-oriented simulation approach for magnetron sputtering deposition was proposed in this work. In this integrated multiscale simulation, the sputtered atom transport was modeled using the Monte Carlo (MC) and molecular dynamics (MD) coupling method, and the deposition of sputtered atoms was simulated using the MD method. This application-oriented simulation approach was used to simulate the growth of Cu/Si(100) thin films at different sputtering pressures. The experimental results unveiled that, as the sputtering pressure decreased from 2 to 0.15 Pa, the surface roughness of Cu thin films gradually decreased; (111)-oriented grains were dominant in Cu thin films and the crystal quality of the Cu thin film was gradually improved. The simulation results were consistent with the experimental characterization results. The simulation results revealed that the transformation of the film growth mode from the Volmer-Weber growth mode to the two-dimensional layered growth mode resulted in a decrease in the surface roughness of Cu thin films; the increase in the amorphous compound CuSix and the hcp copper silicide with the decrease in the sputtering pressure was responsible for the improvement of the crystal quality of the Cu thin film. This work proposed a more realistic, integrated simulation scheme for magnetron sputtering deposition, providing theoretical guidance for the efficient preparation of high-quality sputtered films.