Effect of the nature of the solid substrate on spatially heterogeneous activated dynamics in glass forming supported films.

We extend the force-level elastically collective nonlinear Langevin equation theory to treat the spatial gradients of the alpha relaxation time and glass transition temperature, and the corresponding film-averaged quantities, to the geometrically asymmetric case of finite thickness supported films with variable fluid-substrate coupling. The latter typically nonuniversally slows down motion near the solid-liquid interface as modeled via modification of the surface dynamic free energy caging constraints that are spatially transferred into the film and which compete with the accelerated relaxation gradient induced by the vapor interface. Quantitative applications to the foundational hard sphere fluid and a polymer melt are presented. The strength of the effective fluid-substrate coupling has very large consequences for the dynamical gradients and film-averaged quantities in a film thickness and thermodynamic state dependent manner. The interference of the dynamical gradients of opposite nature emanating from the vapor and solid interfaces is determined, including the conditions for the disappearance of a bulk-like region in the film center. The relative importance of surface-induced modification of local caging vs the generic truncation of the long range collective elastic component of the activation barrier is studied. The conditions for the accuracy and failure of a simple superposition approximation for dynamical gradients in thin films are also determined. The emergence of near substrate dead layers, large gradient effects on film-averaged response functions, and a weak non-monotonic evolution of dynamic gradients in thick and cold films are briefly discussed. The connection of our theoretical results to simulations and experiments is briefly discussed, as is the extension to treat more complex glass-forming systems under nanoconfinement.

Nature of dynamic gradients, glass formation, and collective effects in ultrathin freestanding films.

Molecular, polymeric, colloidal, and other classes of liquids can exhibit very large, spatially heterogeneous alterations of their dynamics and glass transition temperature when confined to nanoscale domains. Considerable progress has been made in understanding the related problem of near-interface relaxation and diffusion in thick films. However, the origin of "nanoconfinement effects" on the glassy dynamics of thin films, where gradients from different interfaces interact and genuine collective finite size effects may emerge, remains a longstanding open question. Here, we combine molecular dynamics simulations, probing 5 decades of relaxation, and the Elastically Cooperative Nonlinear Langevin Equation (ECNLE) theory, addressing 14 decades in timescale, to establish a microscopic and mechanistic understanding of the key features of altered dynamics in freestanding films spanning the full range from ultrathin to thick films. Simulations and theory are in qualitative and near-quantitative agreement without use of any adjustable parameters. For films of intermediate thickness, the dynamical behavior is well predicted to leading order using a simple linear superposition of thick-film exponential barrier gradients, including a remarkable suppression and flattening of various dynamical gradients in thin films. However, in sufficiently thin films the superposition approximation breaks down due to the emergence of genuine finite size confinement effects. ECNLE theory extended to treat thin films captures the phenomenology found in simulation, without invocation of any critical-like phenomena, on the basis of interface-nucleated gradients of local caging constraints, combined with interfacial and finite size-induced alterations of the collective elastic component of the structural relaxation process.

Toward better understanding of the high-pressure structural transformation in beryllium by the statistical moment method.

Beryllium is a vital alkaline-earth metal for plasma physics, space science, and nuclear technology. Unfortunately, its accurate phase diagram is clouded by many controversial results, even though solid beryllium can only exist with hcp or bcc crystalline structures. Herein, we offer a simple quantum-statistical solution to the above problem. Our core idea is to develop the moment expansion technique to determine the Helmholtz free energy under extreme conditions. This strategy helps elucidate the underlying correlation among symmetric characteristics, vibrational excitations, and physical stabilities. In particular, our analyses reveal that the appearance of anharmonic effects forcefully straightens up the hcp-bcc boundary. This phenomenon explains why it has been difficult to detect bcc signatures via diamond-anvil-cell measurements. Besides, we modify the work-heat equivalence principle to quickly obtain the high-pressure melting profile from the room-temperature equation of state. The hcp-bcc-liquid triple point of beryllium is found at 165 GPa and 4559 K. Our theoretical findings agree excellently with cutting-edge ab initio simulations adopting the phonon quasiparticle method and the thermodynamic integration. Finally, we consider the principal Hugoniot curve and its secondary branches to explore the behaviors of beryllium under shock compression. Our predictions would be advantageous for designing inertial-confinement-fusion experiments.

Theoretical insights into non-Arrhenius behaviors of thermal vacancies in anharmonic crystals.

Vacancies are prevalent point defects in crystals, but their thermal responses are elusive. Herein, we formulate a simple theoretical model to shed light on the vacancy evolution during heating. Vibrational excitations are thoroughly investigated via moment recurrence techniques in quantum statistical mechanics. On that basis, we carry out numerical analyses for Ag, Cu, and Ni with the Sutton-Chen many-body potential. Our results reveal that the well-known Arrhenius law is insufficient to describe the proliferation of vacancies. Specifically, anharmonic effects lead to a strong nonlinearity in the Gibbs energy of vacancy formation. Our physical picture is well supported by previous simulations and experiments.

Theory of Pressure-Induced Rejuvenation and Strain Hardening in Metallic Glasses.

We theoretically investigate high-pressure effects on the atomic dynamics of metallic glasses. The theory predicts compression-induced rejuvenation and the resulting strain hardening that have been recently observed in metallic glasses. Structural relaxation under pressure is mainly governed by local cage dynamics. The external pressure restricts the dynamical constraints and slows down the atomic mobility. In addition, the compression induces a rejuvenated metastable state (local minimum) at a higher energy in the free-energy landscape. Thus, compressed metallic glasses can rejuvenate and the corresponding relaxation is reversible. This behavior leads to strain hardening in mechanical deformation experiments. Theoretical predictions agree well with experiments.

Coupling between structural relaxation and diffusion in glass-forming liquids under pressure variation.

We theoretically investigate structural relaxation and activated diffusion of glass-forming liquids at different pressures using both Elastically Collective Nonlinear Langevin Equation (ECNLE) theory and molecular dynamics (MD) simulations. An external pressure restricts local motions of a single molecule within its cage and triggers slowing down of cooperative mobility. While the ECNLE theory and simulations generally predict a monotonic increase of the glass transition temperature and dynamic fragility with pressure, the simulations indicate a decrease of fragility as pressures above 1000 bar. The structural relaxation time is found to be linearly coupled with the inverse diffusion constant. Remarkably, this coupling is independent of compression. The theoretical calculations quantitatively agree well with the simulations and are also consistent with prior work.

Theoretical Model for the Structural Relaxation Time in Coamorphous Drugs.

We propose a simple approach to investigate the structural relaxation time and glass transition of amorphous drugs. Amorphous materials are modeled as a set of equal sized hard spheres. The structural relaxation time over many decades in hard-sphere fluids is theoretically calculated using the elastically collective nonlinear Langevin equation theory associated with Kramer's theory. Then, new thermal mapping from a real material to an effective hard-sphere fluid provides temperature-dependent relaxation time, which can be compared to experiments. Numerical results quantitatively agree with previous experiments for pharmaceutical binary mixtures having different weight ratios. We carry out experiments to test our calculations for an ezetimibe-simvastatin-Kollidon VA64 mixture. Our approach would provide a simple but comprehensive description of glassy dynamics in amorphous composites.

Confinement effects on the solar thermal heating process of TiN nanoparticle solutions.

We propose a theoretical approach to describe quantitatively the heating process in aqueous solutions of dispersed TiN nanoparticles under solar illumination. The temperature gradients of the solution with different concentrations of the TiN nanoparticles are calculated when confinement effects of the container on the solar absorption are taken into account. We find that the average penetration of solar radiation into the solution is significantly reduced upon increasing the nanoparticle concentration. At high concentrations, our numerical results show that photons are localized near the surface of the solution. Moreover, the heat energy balance equation at the vapor-liquid interface is used to describe the solar steam generation. The theoretical time dependence of temperature rise and vaporization weight losses is consistent with experiments. Our calculations give strong evidence that the substantially localized heating near the vapor-liquid interface is the main reason for the more efficient steam generation process by floating plasmonic membranes when compared to randomly dispersed nanoparticles. The validated theoretical model suggests that our approach can be applied towards new predictions and other experimental data descriptions.

Theory of the spatial transfer of interface-nucleated changes of dynamical constraints and its consequences in glass-forming films.

We formulate a new theory for how caging constraints in glass-forming liquids at a surface or interface are modified and then spatially transferred, in a layer-by-layer bootstrapped manner, into the film interior in the context of the dynamic free energy concept of the Nonlinear Langevin Equation (NLE) theory approach. The dynamic free energy at any mean location (cage center) involves contributions from two adjacent layers where confining forces are not the same. At the most fundamental level of the theory, the caging component of the dynamic free energy varies essentially exponentially with distance from the interface, saturating deep enough into the film with a correlation length of modest size and weak sensitivity to the thermodynamic state. This imparts a roughly exponential spatial variation of all the key features of the dynamic free energy required to compute gradients of dynamical quantities including the localization length, jump distance, cage barrier, collective elastic barrier, and alpha relaxation time. The spatial gradients are entirely of dynamical, not structural or thermodynamic, origin. The theory is implemented for the hard sphere fluid and diverse interfaces which can be a vapor, a rough pinned particle solid, a vibrating (softened) pinned particle solid, or a smooth hard wall. Their basic description at the level of the spatially heterogeneous dynamic free energy is identical, with the crucial difference arising from the first layer where dynamical constraints can be weakened, softened, or hardly changed depending on the specific interface. Numerical calculations establish the spatial dependence and fluid volume fraction sensitivity of the key dynamical property gradients for five different model interfaces. A comparison of the theoretical predictions for the dynamic localization length and glassy modulus with simulations and experiments for systems with a vapor interface reveals good agreement. The present advance sets the stage for using the Elastically Collective NLE theory to make quantitative predictions for the alpha relaxation time gradient, decoupling phenomena, Tg gradient, and many film-averaged properties of both model and experimental (colloids, molecules, and polymers) systems with diverse interfaces and chemical makeup.

Theory of activated glassy dynamics in randomly pinned fluids.

We generalize the force-level, microscopic, Nonlinear Langevin Equation (NLE) theory and its elastically collective generalization [elastically collective nonlinear Langevin equation (ECNLE) theory] of activated dynamics in bulk spherical particle liquids to address the influence of random particle pinning on structural relaxation. The simplest neutral confinement model is analyzed for hard spheres where there is no change of the equilibrium pair structure upon particle pinning. As the pinned fraction grows, cage scale dynamical constraints are intensified in a manner that increases with density. This results in the mobile particles becoming more transiently localized, with increases of the jump distance, cage scale barrier, and NLE theory mean hopping time; subtle changes of the dynamic shear modulus are predicted. The results are contrasted with recent simulations. Similarities in relaxation behavior are identified in the dynamic precursor regime, including a roughly exponential, or weakly supra-exponential, growth of the alpha time with pinning fraction and a reduction of dynamic fragility. However, the increase of the alpha time with pinning predicted by the local NLE theory is too small and severely so at very high volume fractions. The strong deviations are argued to be due to the longer range collective elasticity aspect of the problem which is expected to be modified by random pinning in a complex manner. A qualitative physical scenario is offered for how the three distinct aspects that quantify the elastic barrier may change with pinning. ECNLE theory calculations of the alpha time are then presented based on the simplest effective-medium-like treatment for how random pinning modifies the elastic barrier. The results appear to be consistent with most, but not all, trends seen in recent simulations. Key open problems are discussed with regard to both theory and simulation.

Repulsive interactions of a lipid membrane with graphene in composite materials.

The van der Waals interaction between a lipid membrane and a substrate covered by a graphene sheet is investigated using the Lifshitz theory. The reflection coefficients are obtained for a layered planar system submerged in water. The dielectric response properties of the involved materials are also specified and discussed. Our calculations show that a graphene covered substrate can repel the biological membrane in water. This is attributed to the significant changes in the response properties of the system due to the monolayer graphene. It is also found that the van der Waals interaction is mostly dominated by the presence of graphene, while the role of the particular substrate is secondary.

Electric potential profile of a spherical soft particle with a charged core.

The electrostatic potential profile of a spherical soft particle is derived by solving the Poisson-Boltzmann equations on a spherical system both numerically and analytically. The soft particle is assumed to consist of an ion-permeable charged outer layer and a non-permeable charged core with constant charged density. The contribution of the core to the potential profile is calculated for different charges and dielectric constants. Our results show that the charged core heavily influences the local potential within the soft particle. By contrast, the potential distribution outside the particle in the salt solution is found to be weakly dependent on the core features. These findings are consistent with previous experiments showing the minor impact of the core of the MS2 virus on its overall electrical properties. Our studies also indicate that while a change in temperature from 290 K to 310 K only slightly varies the potential, the ionic strength in the range of 1-600 mM has a significant effect on the potential profile. Our studies would provide good understanding for experimental research in the field of biophysics and nanomedicine.

Screening and collective effects in randomly pinned fluids: a new theoretical framework.

We propose a theoretical framework for the dynamics of bulk isotropic hard-sphere systems in the presence of randomly pinned particles and apply this theory to supercooled water to validate it. Structural relaxation is mainly governed by local and non-local activated process. As the pinned fraction grows, a local caging constraint becomes stronger and the long range collective aspect of relaxation is screened by immobile obstacles. Different responses of the local and cooperative motions results in subtle predictions for how the alpha relaxation time varies with pinning and density. Our theoretical analysis for the relaxation time of water with pinned molecules quantitatively well describe previous simulations. In addition, the thermal dependence of relaxation for unpinned bulk water is also consistent with prior computational and experimental data.

Confinement Effects on the Spatially Inhomogeneous Dynamics in Metallic Glass Films.

This work develops the elastically collective nonlinear Langevin equation theory to investigate, for the first time, the glassy dynamics in capped metallic glass thin films. Finite-size effects on the spatial gradient of structural relaxation time and glass transition temperature (Tg) are calculated at different temperatures and vitrification criteria. Molecular dynamics is significantly slowed down near rough solid surfaces, and the dynamics at location far from the interfaces is sped up. In thick films, the mobility gradient normalized by the bulk value obeys the double-exponential form since interference effects between two surfaces are weak. Reducing the film thickness induces a strong dynamic coupling between two surfaces and flattens the relaxation gradient. The normalized gradient of the glass transition temperature is independent of vitrification time scale criterion and can be fitted by a superposition function as the films are not ultrathin. The local fragility is found to remain unchanged with location. This finding suggests that one can use Angell plots of bulk relaxation time and the Tg spatial gradient to characterize glassy dynamics in metallic glass films. Our computational results agree well with experimental data and simulation.

Approaching Polymer Dynamics Combining Artificial Neural Networks and Elastically Collective Nonlinear Langevin Equation.

The analysis of structural relaxation dynamics of polymers gives an insight into their mechanical properties, whose characterization is used to qualify a given material for its practical scope. The dynamics are usually expressed in terms of the temperature dependence of the relaxation time, which is only available through time-consuming experimental processes following polymer synthesis. However, it would be advantageous to estimate their dynamics before synthesizing them when designing new materials. In this work, we propose a combined approach of artificial neural networks and the elastically collective nonlinear Langevin equation (ECNLE) to estimate the temperature dependence of the main structural relaxation time of polymers based only on the knowledge of the chemical structure of the corresponding monomer.

A comparative study of receptor interactions between SARS-CoV and SARS-CoV-2 from molecular modeling.

The pandemic of COVID-19 severe acute respiratory syndrome, which was fatal for millions of people worldwide, triggered the race to understand in detail the molecular mechanisms of this disease. In this work, the differences of interactions between the SARS-CoV/SARS-CoV-2 Receptor binding domain (RBD) and the human Angiotensin Converting Enzyme 2 (ACE2) receptor were studied using in silico tools. Our results show that SARS-CoV-2 RBD is more stable and forms more interactions with ACE2 than SARS-CoV. At its interface, three stable binding patterns are observed and named red-K31, green-K353 and blue-M82 according to the central ACE2 binding residue. In SARS-CoV instead, only the first two binding patches are persistently formed during the MD simulation. Our MM/GBSA calculations indicate the binding free energy difference of about 2.5 kcal/mol between SARS-CoV-2 and SARS-CoV which is compatible with the experiments. The binding free energy decomposition points out that SARS-CoV-2 RBD-ACE2 interactions of the red-K31 ([Formula: see text]) and blue-M82 ([Formula: see text]) patterns contribute more to the binding affinity than in SARS-CoV ([Formula: see text] for red-K31), while the contribution of the green-K353 pattern is very similar in the two strains ([Formula: see text] and [Formula: see text] for SARS-CoV-2 and SARS-CoV, respectively). Five groups of mutations draw our attention at the RBD-ACE2 binding interface, among them, the mutation -PPA469-471/GVEG482-485 has the most important and favorable impact on SARS-CoV-2 binding to the ACE2 receptor. These results, highlighting the molecular differences in the binding between the two viruses, contribute to the common knowledge about the new corona virus and to the development of appropriate antiviral treatments, addressing the necessity of ongoing pandemics.

Effects of surface charge and environmental factors on the electrostatic interaction of fiber with virus-like particle: A case of coronavirus.

We propose a theoretical model to elucidate intermolecular electrostatic interactions between a virus and a substrate. Our model treats the virus as a homogeneous particle having surface charge and the polymer fiber of the respirator as a charged plane. Electric potentials surrounding the virus and fiber are influenced by the surface charge distribution of the virus. We use Poisson-Boltzmann equations to calculate electric potentials. Then, Derjaguin's approximation and a linear superposition of the potential function are extended to determine the electrostatic force. In this work, we apply this model for coronavirus or SARS-CoV-2 case and numerical results quantitatively agree with prior simulation. We find that the influence of fiber's potential on the surface charge of the virus is important and is considered in interaction calculations to obtain better accuracy. The electrostatic interaction significantly decays with increasing separation distance, and this curve becomes steeper when adding more salt. Although the interaction force increases with heating, one can observe the repulsive-attractive transition when the environment is acidic.

Determination of Young's Modulus of Active Pharmaceutical Ingredients by Relaxation Dynamics at Elevated Pressures.

A new approach is theoretically proposed to study the glass transition of active pharmaceutical ingredients and a glass-forming anisotropic molecular liquid at high pressures. We describe amorphous materials as a fluid of hard spheres. Effects of nearest neighbor interactions and cooperative motions of particles on glassy dynamics are quantified through a local and collective elastic barrier calculated using the elastically collective nonlinear Langevin equation theory. Inserting two barriers into Kramer's theory gives the structural relaxation time. Then, we formulate a new mapping based on the thermal expansion process under pressure to intercorrelate particle density, temperature, and pressure. This analysis allows us to determine the pressure and temperature dependence of α relaxation. From this, we estimate the effective elastic modulus of amorphous materials and capture the effects of conformation on the relaxation process. Remarkably, our theoretical results agree well with experiments.

Theory of Structural and Secondary Relaxation in Amorphous Drugs under Compression.

Compression effects on alpha and beta relaxation process of amorphous drugs are theoretically investigated by developing the elastically collective nonlinear Langevin equation theory. We describe the structural relaxation as a coupling between local and nonlocal activated process. Meanwhile, the secondary beta process is mainly governed by the nearest-neighbor interactions of a molecule. This assumption implies the beta relaxation acts as a precursor of the alpha relaxation. When external pressure is applied, a small displacement of a molecule is additionally exerted by a pressure-induced mechanical work in the dynamic free energy, which quantifies interactions between a molecule with its nearest neighbors. The local dynamics has more restriction and it induces stronger effects of collective motions on single-molecule dynamics. Thus, the alpha and beta relaxation times are significantly slowed down with increasing compression. We apply this approach to determine the temperature and pressure dependence of the alpha and beta relaxation time for curcumin, glibenclamide, and indomethacin, and compare numerical results with prior experimental studies. Both qualitative and quantitative agreement between theoretical calculations and experiments validate our assumptions and reveal their limitations. Our approach would pave the way for the development of the drug formulation process.

Cooperative nanoparticle self-assembly and photothermal heating in a flexible plasmonic metamaterial.

We theoretically investigate equilibrium behaviors and photothermal effects of a flexible plasmonic metamaterial composed of aramid nanofibers and gold nanoparticles. The fiber matrix is considered as an external field to reconfigure a nanoparticle assembly. We find that the heating process tunes particle-particle and fiber-particle interactions, which alter adsorption of nanoparticles on fiber surfaces or clustering in pore spaces. Thus, it is possible to control the nanoparticle self-assembly by laser illumination. Gold nanoparticles strongly absorb radiations and efficiently dissipate absorbed energy into heat. By solving the heat transfer equation associated with an effective medium approximation, we calculate the spatial temperature rise. Remarkably, our theoretical results quantitatively agree with prior experiments. This indicates that we can ignore plasmonic coupling effects induced by particle clustering. Effects of the laser spot size and intensity on the photothermal heating are also discussed.