Search

Hydration and, in particular, the coordination number of a metal ion is of paramount importance as it defines many of its (bio)physicochemical properties. It is not only essential for understanding its behavior in aqueous solutions but also determines the metal ion reference state and its binding energy to (bio)molecules. In this paper, for divalent metal cations Ca 2 + , Cd 2 + , Cu 2 + , Fe 2 + , Hg 2 + , Mg 2 + , Ni 2 + , Pb 2 + , and Zn 2 + , we compare two approaches for predicting hydration numbers: (1) a mixed explicit/continuum DFT-D3//COSMO-RS solvation model and (2) density functional theory based ab initio molecular dynamics. The former approach is employed to calculate the Gibbs free energy change for the sequential hydration reactions, starting from [M( H 2 O)2]2+ aqua complexes to [M( H 2 O) 9]2+, allowing explicit water molecules to bind in the first or second coordination sphere and determining the most stable [M( H 2 O) n]2+ structure. In the latter approach, the hydration number is obtained by integrating the ion–water radial distribution function. With a couple of exceptions, the metal ion hydration numbers predicted by the two approaches are in mutual agreement, as well as in agreement with the experimental data. [ABSTRACT FROM AUTHOR]

31P nuclear magnetic resonance (NMR) chemical shifts were shown to be very sensitive to the basis set used at the geometry optimization stage. Commonly used energy-optimized basis sets for a phosphorus atom containing only one polarization d-function were shown to be unable to provide correct equilibrium geometries for the calculations of phosphorus chemical shifts. The use of basis sets with at least two polarization d-functions on a phosphorus atom is strongly recommended. In this paper, an idea of creating the basis sets purposed for the geometry optimization that provide the least possible error coming from the geometry factor of accuracy in the resultant NMR shielding constants is proposed. The property-energy consisted algorithm with the target function in the form of the molecular energy gradient relative to P–P bond lengths was applied to create new geometry-oriented pecG-n (n = 1, 2) basis sets for a phosphorus atom. New basis sets have demonstrated by far superior performance as compared to the other commonly used energy-optimized basis sets in massive calculations of 31P NMR chemical shifts carried out at the gauge-including atomic orbital-coupled cluster singles and doubles/pecS-2 level of the theory by taking into account solvent, vibrational, and relativistic corrections. [ABSTRACT FROM AUTHOR]

A Markov state model is a powerful tool that can be used to track the evolution of populations of configurations in an atomistic representation of a protein. For a coarse-grained linear chain model with discontinuous interactions, the transition rates among states that appear in the Markov model when the monomer dynamics is diffusive can be determined by computing the relative entropy of states and their mean first passage times, quantities that are unchanged by the specification of the energies of the relevant states. In this paper, we verify the folding dynamics described by a diffusive linear chain model of the crambin protein in three distinct solvent systems, each differing in complexity: a hard-sphere solvent, a solvent undergoing multi-particle collision dynamics, and an implicit solvent model. The predicted transition rates among configurations agree quantitatively with those observed in explicit molecular dynamics simulations for all three solvent models. These results suggest that the local monomer–monomer interactions provide sufficient friction for the monomer dynamics to be diffusive on timescales relevant to changes in conformation. Factors such as structural ordering and dynamic hydrodynamic effects appear to have minimal influence on transition rates within the studied solvent densities. [ABSTRACT FROM AUTHOR]

Chiral thin films showing electronic and plasmonic circular dichroism (CD) are intensively explored for optoelectronic applications. The most studied chiral organic films are the composites exhibiting a helical geometry, which often causes entanglement of circular optical properties with unwanted linear optical effects (linearly polarized absorption or refraction). This entanglement limits tunability and often translates to a complex optical response. This paper describes chiral films based on dark conglomerate, sponge-like, liquid crystal films, which go beyond the usual helical type geometry, waiving the problem of linear contributions to chiroptical electronic and plasmonic properties. First, we show that purely organic films exhibit high electronic CD and circular birefringence, as studied in detail using Mueller matrix polarimetry. Analogous linear properties are two orders of magnitude lower, highlighting the benefits of using the bi-isotropic dark conglomerate liquid crystal for chiroptical purposes. Next, we show that the liquid crystal can act as a template to guide the assembly of chemically compatible gold nanoparticles into 3D spiral-like assemblies. The Mueller matrix polarimetry measurements confirm that these composites exhibit both electronic and plasmonic circular dichroisms, while nanoparticle presence is not compromising the beneficial optical properties of the matrix. [ABSTRACT FROM AUTHOR]

The optical surface of extreme ultraviolet (EUV) lithography machines is highly vulnerable to contamination by hydrocarbons, resulting in the formation of carbon deposits that significantly degrade the quality and efficiency of lithography. The dynamic gas lock (DGL) has been proven as an effective approach to alleviate carbon deposition. However, the majority of existing studies on carbon deposition neglect the influence of the DGL. This paper is dedicated to investigating the phenomena of hydrocarbon adsorption, desorption, and cleavage with considering the effects of the DGL. A comprehensive mathematical model of the carbon deposition process is established, and the impact of radiation intensity, temperature, and hydrocarbon types on the depositing rate is considered. The results suggest that the primary cause of carbon deposition is the direct cracking of hydrocarbons induced by photons with a wavelength range between 12.5 and 14.5 nm. Additionally, it has been observed that the carbon deposition rate decreases exponentially as clean gas flow increases when EUV radiation intensity exceeds 50 mW/mm2. Conversely, at low EUV radiation intensity, clean gas flow has little effect on the carbon deposition rate. An effective approach to mitigate carbon deposition is to elevate the temperature of the optical surface and employ light hydrocarbon materials in the EUV process. [ABSTRACT FROM AUTHOR]

Previous studies have shown that NaB6, KB6, and RbB6 adopting Pm 3 ̄ m are superconductors with a relatively high Tc under ambient conditions. In this paper, we conducted systematic structural and related properties research on CsB6 through a genetic evolution algorithm and total energy calculations based on density functional theory between 0 and 20 GPa. Our results reveal a cubic Pm 3 ̄ m CsB6, which is dynamically stable under the pressures we studied. We systematically calculated the formation enthalpies, electronic properties, and superconducting properties of Pm 3 ̄ m MB6 (M = Na, K, Rb, Cs). They all exhibit metallic features, and boron has high contributions to band structures, density of states, and electron–phonon coupling (EPC). The calculated results about the Helmholtz free energy difference of Pm 3 ̄ m CsB6 at 0, 10, and 20 GPa indicate that it is stable upon chemical decomposition (decomposition to simple substances Cs and B) from 0 to 400 K. The phonon density of states indicates that boron atoms occupy the high frequency area. The EPC results show that Pm 3 ̄ m CsB6 is a superconductor with Tc = 11.7 K at 0 GPa, close to NaB6 (13.1 K), KB6 (11.7 K), and RbB6 (11.3 K) at 0 GPa in our work, which indicates that boron atoms play an essential role in superconductivity: vibrations of B6 regular octagons lead to the high Tc of Pm 3 ̄ m MB6. Our work about Pm 3 ̄ m hexaborides provides a supplementary study on the borides of the group IA elements (without Fr and Li) and has an important guiding significance for the experimental synthesis of CsB6. [ABSTRACT FROM AUTHOR]

This paper presents a general second-quantized form of a permutation operator interchanging n pairs of electrons between interacting subsystems in the framework of the symmetry-adapted perturbation theory (SAPT). We detail the procedure for constructing this operator through the consecutive multiplication of single-pair permutation operators. This generalized form of the permutation operator has enabled the derivation of universal formulas for S2n approximations of the exchange energies in the first and second order of the interaction operator. We present expressions for corrections of S4 approximations and assess its efficacy on a selection of systems anticipated to exhibit a slowly converging overlap expansion. Additionally, we outline a method to sum the overlap expansion series to infinity in second-quantization, up to the second order in V. This new approach offers an alternative to the existing formalism based on density-matrix formulations. When combined with a symbolic algebra program for automated derivations, it paves the way for advancements in SAPT theory, particularly for intricate wavefunction theories. [ABSTRACT FROM AUTHOR]

Modern 4-wave mixing spectroscopies are expensive to obtain experimentally and computationally. In certain cases, the unfavorable scaling of quantum dynamics problems can be improved using a generalized quantum master equation (GQME) approach. However, the inclusion of multiple (light–matter) interactions complicates the equation of motion and leads to seemingly unavoidable cubic scaling in time. In this paper, we present a formulation that greatly simplifies and reduces the computational cost of previous work that extended the GQME framework to treat arbitrary numbers of quantum measurements. Specifically, we remove the time derivatives of quantum correlation functions from the modified Mori–Nakajima–Zwanzig framework by switching to a discrete-convolution implementation inspired by the transfer tensor approach. We then demonstrate the method's capabilities by simulating 2D electronic spectra for the excitation-energy-transfer dimer model. In our method, the resolution of data can be arbitrarily coarsened, especially along the t2 axis, which mirrors how the data are obtained experimentally. Even in a modest case, this demands O (1 0 3 ) fewer data points. We are further able to decompose the spectra into one-, two-, and three-time correlations, showing how and when the system enters a Markovian regime where further measurements are unnecessary to predict future spectra and the scaling becomes quadratic. This offers the ability to generate long-time spectra using only short-time data, enabling access to timescales previously beyond the reach of standard methodologies. [ABSTRACT FROM AUTHOR]

This paper reports the energies and charge and spin distributions of the low-lying excited states in singlet and triplet N2V defects in diamond from direct Δ-SCF calculations based on Gaussian orbitals within the B3LYP, PBE0, and HSE06 functionals. They assign the observed absorption at 2.463 eV, first reported by Davies et al. [Proc. R. Soc. London 351, 245 (1976)], to the excitation of a N(sp3) lone-pair electron in the singlet and triplet states, respectively, with estimates of ∼1.1 eV for that of the unpaired electrons, C(sp3). In both cases, the excited states are predicted to be highly local and strongly excitonic with 81% of the C(sp3) and 87% of the N(sp3) excited charges localized at the three C atoms nearest neighbor (nn) to the excitation sites. Also reported are the higher excited gap states of both the N lone pair and C unpaired electron. Calculated excitation energies of the bonding sp3 hybrids of the C atoms nn to the four inner atoms are close to that of the bulk, which indicates that the N2V defect is largely a local defect. The present results are in broad agreement with those reported by Udvarhelyi et al. [Phys. Rev. B 96, 155211 (2017)] from plane wave HSE06 calculations, notably for the N lone pair excitation energy, for which both predict an energy of ∼2.7 eV but with a difference of ∼0.5 eV for the excitation of the unpaired electron. [ABSTRACT FROM AUTHOR]

Many methods to accelerate sampling of molecular configurations are based on the idea that temperature can be used to accelerate rare transitions. These methods typically compute equilibrium properties at a target temperature using reweighting or through Monte Carlo exchanges between replicas at higher temperatures. A recent paper [G. M. Rotskoff and E. Vanden-Eijnden, Phys. Rev. Lett. 122, 150602 (2019)] demonstrated that accurate equilibrium densities of states can also be computed through a nonequilibrium "quench" process, where sampling is performed at a higher temperature to encourage rapid mixing and then quenched to lower energy states with dissipative dynamics. Here, we provide an implementation of the quench dynamics in LAMMPS and evaluate a new formulation of nonequilibrium estimators for the computation of partition functions or free energy surfaces (FESs) of molecular systems. We show that the method is exact for a minimal model of N-independent harmonic springs and use these analytical results to develop heuristics for the amount of quenching required to obtain accurate sampling. We then test the quench approach on alanine dipeptide, where we show that it gives an FES that is accurate near the most stable configurations using the quench approach but disagrees with a reference umbrella sampling calculation in high FE regions. We then show that combining quenching with umbrella sampling allows the efficient calculation of the free energy in all regions. Moreover, by using this combined scheme, we obtain the FES across a range of temperatures at no additional cost, making it much more efficient than standard umbrella sampling if this information is required. Finally, we discuss how this approach can be extended to solute tempering and demonstrate that it is highly accurate for the case of solvated alanine dipeptide without any additional modifications. [ABSTRACT FROM AUTHOR]

The network structure within polymers significantly influences their mechanical properties, including their strength, toughness, and fatigue resistance. All-atom molecular dynamics (AAMD) simulations offer a method to investigate the energy dissipation mechanism within polymers during deformation and fracture; Such an approach is, however, computationally inefficient when used to analyze polymers with complex network structures, such as the common chemically double-networked hydrogels. Alternatively, coarse-grained molecular dynamics (CGMD) models, which reduce the computational degrees of freedom by concentrating a set of adjacent atoms into a coarse-grained bead, can be employed. In CGMD simulations, a coarse-grained force field (CGFF) is a critical factor affecting the simulation accuracy. In this paper, we proposed a data-based method for predicting the CGFF parameters to improve the simulation efficiency of complex cross-linked network in polymers. Here, we utilized a typical chemically double-networked hydrogel as an example. An artificial neural network was selected, and it was trained with the tensile stress–strain data from the CGMD simulations using different CGFF parameters. The CGMD simulations using the predicted CGFF parameters show good agreement with the AAMD simulations and are almost fifty times faster. The data-inspired CGMD model presented here broadens the applicability of molecular dynamics simulations to cross-linked polymers and has the potential to provide insights that will aid the design of polymers with desirable mechanical properties. [ABSTRACT FROM AUTHOR]

Stimulated Raman Adiabatic Passage (STIRAP) and its variants, such as M-type chainwise-STIRAP, allow for efficiently transferring the populations in a multilevel system and have widely been used to prepare molecules in their rovibrational ground state. However, their transfer efficiencies are generally imperfect. The main obstacle is the presence of losses and the requirement to make the dynamics adiabatic. To this end, in the present paper, a new theoretical method is proposed for the efficient and robust creation and detection of deeply bound molecules in three-level Λ-type and five-level M-type systems via "Invariant-based shortcut-to-adiabaticity." In the regime of large detunings, we first reduce the dynamics of three- and five-level molecular systems to those of effective two- and three-level counterparts. By doing so, the major molecular losses from the excited states can be well suppressed. Consequently, the effective two-level counterpart can be directly compatible with two different "Invariant-based Inverse Engineering" protocols; the results show that both protocols give a comparable performance and have a good experimental feasibility. For the effective three-level counterpart, by considering a relation among the four incident pulses, we show that this model can be further generalized to an effective Λ-type one with the simplest resonant coupling. This generalized model permits us to borrow the "Invariant-based Inverse Engineering" protocol from a standard three-level Λ-type system to a five-level M-type system. Numerical calculations show that the weakly bound molecules can be efficiently transferred to their deeply bound states without strong laser pulses, and the stability against parameter variations is well preserved. Finally, the detection of ultracold deeply bound molecules is discussed. [ABSTRACT FROM AUTHOR]

Recent experiments by Brückner et al. [Science 380, 1357 (2023)] have observed an anomalous chain length dependence of the time of near approach of widely separated pairs of genomic elements on transcriptionally active chromosomal DNA. In this paper, I suggest that the anomaly may have its roots in internal friction between neighboring segments on the DNA backbone. The basis for this proposal is a model of chain dynamics formulated in terms of a continuum scaled Brownian walk (sBw) of polymerization index N. The sBw is an extension of the simple Brownian walk model widely used in path integral calculations of polymer properties, differing from it in containing an additional parameter H (the Hurst index) that can be tuned to produce varying degrees of correlation between adjacent monomers. A calculation using the sBw of the mean time τc for chain closure predicts—under the Wilemski–Fixman approximation for diffusion-controlled reactions—that at early times, τc varies as the 2/3 power of N, in close agreement with the findings of the Brückner et al. study. Other scaling relations of that study, including those related to the probability of loop formation and the mean square displacements of terminal monomers, are also satisfactorily accounted for by the model. [ABSTRACT FROM AUTHOR]

A major limitation of amine-based post-combustion carbon capture technology is the necessity to regenerate amines at high temperatures, which dramatically increases operating costs. This paper concludes the effect of solvent choice as a possible route to modify the thermodynamics and kinetics characterizing the involved amine regeneration reactions and discusses whether these modifications can be economically beneficial. We report experimentally benchmarked computational chemistry calculations of monoethanolamine regeneration reactions employing aqueous and non-aqueous solvents with a wide range of dielectric constants. Unlike previous studies, our improved computational chemistry framework could accurately reproduce the right experimental activation energy of zwitterion formation. From the thermodynamics and kinetics of the predicted reactions, the use of non-aqueous solvents with small dielectric constants led to reductions in regeneration Gibbs free energies, activation barriers, and enthalpy changes. This can reduce energy consumption and give an opportunity to run desorption columns at relatively lower temperatures, thus offering the possibility of relying on low-grade waste heat as an energy input. [ABSTRACT FROM AUTHOR]

This is an I easier approach to implement into computer code i , because these equations [Eq. (95) of the paper] are simpler than the corresponding EOMs for { I i SB I n i sb , I i SB I n i sb }. CLARIFICATION ON THE NUMERICAL ALGORITHM USED TO PROPAGATE DYNAMICS We want to clarify the numerical algorithm we used to propagate the EOMs and generate all numerical results presented in the paper. In the above expressions, to compute HT ht , we use Eq. (C2) of the paper. [Extracted from the article]

In this paper, we demonstrate that the weak temperature dependence of the structure factor of supercooled liquids, a defining feature of the glass transition, is a consequence of the averaging of the scattering intensity due to angular averaging. We show that the speckle at individual wavevectors, calculated from a simulated glass former, exhibits a Debye–Waller factor with a sufficiently large temperature dependence to represent a structural order parameter capable of distinguishing liquid from glass. We also extract from the speckle intensities a quantity proportional to the variance of the local restraint, i.e., a direct experimental measure of the amplitude of structural heterogeneity. [ABSTRACT FROM AUTHOR]

As an optimal one-dimensional reaction coordinate, the committor function not only describes the probability of a trajectory initiated at a phase space point first reaching the product state before reaching the reactant state but also preserves the kinetics when utilized to run a reduced dynamics model. However, calculating the committor function in high-dimensional systems poses significant challenges. In this paper, within the framework of milestoning, exact expressions for committor functions at two levels of coarse graining are given, including committor functions of phase space point to point (CFPP) and milestone to milestone (CFMM). When combined with transition kernels obtained from trajectory analysis, these expressions can be utilized to accurately and efficiently compute the committor functions. Furthermore, based on the calculated committor functions, an adaptive algorithm is developed to gradually refine the transition state region. Finally, two model examples are employed to assess the accuracy of these different formulations of committor functions. [ABSTRACT FROM AUTHOR]

Modeling excited state processes in molecular crystals is relevant for several applications. A popular approach for studying excited state molecular crystals is to use cluster models embedded in point charges. In this paper, we compare the performance of several embedding models in predicting excited states and S1–S0 optical gaps for a set of crystals from the X23 molecular crystal database. The performance of atomic charges based on ground or excited states was examined for cluster models, Ewald embedding, and self-consistent approaches. We investigated the impact of various factors, such as the level of theory, basis sets, embedding models, and the level of localization of the excitation. We consider different levels of theory, including time-dependent density functional theory and Tamm–Dancoff approximation (TDA) (DFT functionals: ωB97X-D and PBE0), CC2, complete active space self-consistent field, and CASPT2. We also explore the impact of selection of the QM region, charge leakage, and level of theory for the description of different kinds of excited states. We implemented three schemes based on distance thresholds to overcome overpolarization and charge leakage in molecular crystals. Our findings are compared against experimental data, G0W0-BSE, periodic TDA, and optimally tuned screened range-separated functionals. [ABSTRACT FROM AUTHOR]

In a recent paper, Becke et al. [J. Chem. Phys. 158, 151103 (2023)] presented a novel double hybrid density functional, "DH23," whose terms are based on good physics. Its 12 coefficients were trained on the GMTKN55 (general main-group thermochemistry, kinetics, and noncovalent interactions) chemical database of Goerigk et al. [Phys. Chem. Chem. Phys. 19, 32184 (2017)]. The lowest GMTKN55 "WTMAD2" error to date for any hybrid or double hybrid density functional was obtained (1.76 kcal/mol). Here, we make some revisions to DH23 and test its efficacy on reference data beyond GMTKN55, namely, organometallic reaction energies and barrier heights. The results confirm that DH23 is robust outside its training set. In the process, a slightly smaller GMTKN55 WTMAD2 of 1.73 kcal/mol is achieved. [ABSTRACT FROM AUTHOR]

Polymer models describing the statistics of biomolecules under confinement have applications to a wide range of single-molecule experimental techniques and give insight into biologically relevant processes in vivo. In this paper, we determine the transverse position and bending correlation functions for a wormlike chain confined within slits and cylinders (with one and two confined dimensions, respectively) using a mean-field approach that enforces rigid constraints on average. We show the theoretical predictions accurately capture the statistics of a wormlike chain from Monte Carlo simulations in both confining geometries for both weak and strong confinement. We also show that the longitudinal correlation function is accurately computed for a chain confined to a slit and leverages the accuracy of the model to suggest an experimental technique to infer the (often unobservable) transverse statistics from the (directly observable) longitudinal end-to-end distance. [ABSTRACT FROM AUTHOR]

TurboGenius is an open-source Python package designed to fully control ab initio quantum Monte Carlo (QMC) jobs using a Python script, which allows one to perform high-throughput calculations combined with TurboRVB [Nakano et al. J. Phys. Chem. 152, 204121 (2020)]. This paper provides an overview of the TurboGenius package and showcases several results obtained in a high-throughput mode. For the purpose of performing high-throughput calculations with TurboGenius, we implemented another open-source Python package, TurboWorkflows, that enables one to construct simple workflows using TurboGenius. We demonstrate its effectiveness by performing (1) validations of density functional theory (DFT) and QMC drivers as implemented in the TurboRVB package and (2) benchmarks of Diffusion Monte Carlo (DMC) calculations for several datasets. For (1), we checked inter-package consistencies between TurboRVB and other established quantum chemistry packages. By doing so, we confirmed that DFT energies obtained by PySCF are consistent with those obtained by TurboRVB within the local density approximation (LDA) and that Hartree–Fock (HF) energies obtained by PySCF and Quantum Package are consistent with variational Monte Carlo energies obtained by TurboRVB with the HF wavefunctions. These validation tests constitute a further reliability check of the TurboRVB package. For (2), we benchmarked the atomization energies of the Gaussian-2 set, the binding energies of the S22, A24, and SCAI sets, and the equilibrium lattice parameters of 12 cubic crystals using DMC calculations. We found that, for all compounds analyzed here, the DMC calculations with the LDA nodal surface give satisfactory results, i.e., consistent either with high-level computational or with experimental reference values. [ABSTRACT FROM AUTHOR]

In this paper, a global and full-dimensional potential energy surface at the 2A″ ground state for the Al + O2 → AlO + O reaction was constructed, for the first time, based on extensive electronic structure calculations using the doubly hybrid functional XYG3 and potential energy surface fittings by neural networks. Details of the reaction paths have been analyzed. It was found that both two intermediates, the cyclic-AlO2 and the linear-OAlO, were able to dissociate to the AlO + O products, and the isomerization process between these two intermediates was controlled by conical intersections between two 2A″ states. Ro-vibrational state resolved integral cross sections have also been calculated at collision energies from 1.0 to 10.0 kcal/mol. The results support the harpooning mechanism in this metal-oxidant-involved reaction. [ABSTRACT FROM AUTHOR]

Instanton-based rate theory is a powerful tool that is used to explore tunneling in many-dimensional systems. Yet, it diverges at the so-called "crossover temperature." Using the uniform semiclassical transmission probability of Kemble [Phys. Rev. 48, 549 (1935)], we showed recently that in one dimension, one might derive a uniform semiclassical instanton rate theory, which has no divergence. In this paper, we generalize this uniform theory to many-dimensional systems. The resulting theory uses the same input as in the previous instanton theory, yet does not suffer from the divergence. The application of the uniform theory to dissipative systems is considered and used to revise Wolynes' well-known analytical expression for the rate [P. G. Wolynes, Phys. Rev. Lett. 47, 968 (1981)] so that it does not diverge at the "crossover temperature." [ABSTRACT FROM AUTHOR]

Chemical and photochemical reactivity, as well as supramolecular organization and several other molecular properties, can be modified by strong interactions between light and matter. Theoretical studies of these phenomena require the separation of the Schrödinger equation into different degrees of freedom as in the Born–Oppenheimer approximation. In this paper, we analyze the electron–photon Hamiltonian within the cavity Born–Oppenheimer approximation (CBOA), where the electronic problem is solved for fixed nuclear positions and photonic parameters. In particular, we focus on intermolecular interactions in representative dimer complexes. The CBOA potential energy surfaces are compared with those obtained using a polaritonic approach, where the photonic and electronic degrees of freedom are treated at the same level. This allows us to assess the role of electron–photon correlation and the accuracy of CBOA. [ABSTRACT FROM AUTHOR]

The fewest switches surface hopping (FSSH) method proposed by Tully in 1990 [Tully, J. Chem. Phys. 93, 1061 (1990)]—along with its many later variations—forms the basis for most practical simulations of molecular dynamics with electronic transitions in realistic systems. Despite its popularity, a rigorous formal derivation of the algorithm has yet to be achieved. In this paper, we derive the energy-conserving momentum jumps employed by FSSH from the perspective of quantum trajectory surface hopping (QTSH) [Martens, J. Phys. Chem. A 123, 1110 (2019)]. In the limit of localized nonadiabatic transitions, simple mathematical and physical arguments allow the FSSH algorithm to be derived from first principles. For general processes, the quantum forces characterizing the QTSH method provide accurate results for nonadiabatic dynamics with rigorous energy conservation, at the ensemble level, within the consistency of the underlying stochastic surface hopping without resorting to the artificial momentum rescaling of FSSH. [ABSTRACT FROM AUTHOR]

In Paper I [H. Francisco, A. C. Cancio, and S. B. Trickey, J. Chem. Phys. 159, 214102 (2023)], we gave a regularization of the Tao–Mo exchange functional that removes the order-of-limits problem in the original Tao–Mo form and also eliminates the unphysical behavior introduced by an earlier regularization while essentially preserving compliance with the second-order gradient expansion. The resulting simplified, regularized (sregTM) functional delivers performance on standard molecular and solid state test sets equal to that of the earlier revised, regularized Tao–Mo functional. Here, we address de-orbitalization of that new sregTM into a pure density functional. We summarize the failures of the Mejía-Rodríguez and Trickey de-orbitalization strategy [Phys. Rev. A 96, 052512 (2017)] when used with both versions. We discuss how those failures apparently arise in the so-called z′ indicator function and in substitutes for the reduced density Laplacian in the parent functionals. Then, we show that the sregTM functional can be de-orbitalized somewhat well with a rather peculiarly parameterized version of the previously used deorbitalizer. We discuss, briefly, a de-orbitalization that works in the sense of reproducing error patterns but that apparently succeeds by cancelation of major qualitative errors associated with the de-orbitalized indicator functions α and z, hence, is not recommended. We suggest that the same issue underlies the earlier finding of comparatively mediocre performance of the de-orbitalized Tao–Perdew–Staroverov–Scuseri functional. Our work demonstrates that the intricacy of such two-indicator functionals magnifies the errors introduced by the Mejía-Rodríguez and Trickey de-orbitalization approach in ways that are extremely difficult to analyze and correct. [ABSTRACT FROM AUTHOR]

The revised, regularized Tao–Mo (rregTM) exchange-correlation density functional approximation (DFA) [A. Patra, S. Jana, and P. Samal, J. Chem. Phys. 153, 184112 (2020) and Jana et al., J. Chem. Phys. 155, 024103 (2021)] resolves the order-of-limits problem in the original TM formulation while preserving its valuable essential behaviors. Those include performance on standard thermochemistry and solid data sets that is competitive with that of the most widely explored meta-generalized-gradient-approximation DFAs (SCAN and r2SCAN) while also providing superior performance on elemental solid magnetization. Puzzlingly however, rregTM proved to be intractable for de-orbitalization via the approach of Mejía-Rodríguez and Trickey [Phys. Rev. A 96, 052512 (2017)]. We report investigation that leads to diagnosis of how the regularization in rregTM of the z indicator functions (z = the ratio of the von-Weizsäcker and Kohn–Sham kinetic energy densities) leads to non-physical behavior. We propose a simpler regularization that eliminates those oddities and that can be calibrated to reproduce the good error patterns of rregTM. We denote this version as simplified, regularized Tao–Mo, sregTM. We also show that it is unnecessary to use rregTM correlation with sregTM exchange: Perdew–Burke–Ernzerhof correlation is sufficient. The subsequent paper shows how sregTM enables some progress on de-orbitalization. [ABSTRACT FROM AUTHOR]

In Paper I [J. R. Elliott, A. J. Schultz, and D. A. Kofke, J. Chem. Phys. 143, 114110 (2015)] of this series, a methodology was presented for computing the coefficients of a power series of the Helmholtz energy in reciprocal temperature, β, through density series based on cluster integral expansions. Previously, power series in β were evaluated by thermodynamic perturbation theory (TPT) using molecular simulation of a reference fluid. The present methodology uses cluster integrals to evaluate coefficients of the density expansion at each individual order of temperature. While Paper I [J. R. Elliott, A. J. Schultz, and D. A. Kofke, J. Chem. Phys. 143, 114110 (2015)] developed this methodology for square well (SW) spheres, the present work extends the methodology to Lennard-Jones (LJ) spheres, where the reference fluid is the Weeks-Chandler-Andersen potential. Comparisons of TPT coefficients computed from cluster integrals to those from molecular simulation show good agreement through third order in β when coefficients are expressed with effective approximants. Notably, the agreement for LJ spheres is much better than for SW spheres although fewer coefficients of the density series (B2–B5) are available than for SW spheres (B2–B6). The coefficients for Bi(β) of the reference fluid are shown to follow a simple relationship to the virial coefficients of hard sphere fluids, corrected for the temperature dependency of the equivalent hard sphere diameter. This lays the foundation for a correlation of the second virial coefficient of LJ spheres B2(β) that extrapolates to infinite order in temperature. This correlation of B2(β) provides a basis for estimating the low density limit of TPT coefficients at all orders in temperature, facilitating a recursive extrapolation formula to estimate TPT coefficients of fourth order and higher over the entire density range. The applicability of the resulting equation of state is demonstrated by computing the thermodynamic properties for LJ spheres and comparing to standard simulation results. [ABSTRACT FROM AUTHOR]

The present paper is the second part of a series of papers aimed at assessing the accuracy of alchemical computational approaches based on nonequilibrium techniques for solvation free energy of organic molecules in the context of molecular dynamics simulations. In Paper I [Procacci, J. Chem. Phys. 151, 144113 (2019)], we dealt with bidirectional estimates of solvation free energies using nonequilibrium approaches. Here, we assess accuracy and precision of unidirectional estimates with the focus on the Gaussian and Jarzynski estimators. We present a very simple methodology to increase the statistics in the work distribution, hence boosting the accuracy and precision of the Jarzynski unidirectional estimates at no extra cost, exploiting the observed decorrelation between the random variables represented by the Lennard-Jones solute-solvent recoupling or decoupling work and by the electrostatic work due to the charging/discharging of the solute in the solvent. [ABSTRACT FROM AUTHOR]

The polymeric Langmuir theory, developed in Paper I [J. Dudowicz et al., J. Chem. Phys. 151, 124706 (2019)], is employed to investigate the influence of van der Waals interactions and chain rigidity on the thermodynamics of the reversible molecular binding to interfaces in one component polymer fluids (polymer melts). Both van der Waals interactions and chain stiffness are found to influence the temperature variation of the surface coverage Θ, the binding transition itself, and the cooperativity of molecular binding. Re-entrancy of the surface coverage Θ(T) is found to arise when the intermolecular interactions are sufficiently attractive to cause a liquid-vapor like phase separation in the interfacial region, a phenomenon that can occur in the binding of both small molecules and polymer chains to surfaces. [ABSTRACT FROM AUTHOR]

We report the vacuum-ultraviolet threshold-photoelectron spectrum of HC5N recorded over a wide spectral range, from 84 000 to 120 000 cm−1, with a 120 cm−1 spectral resolution, better than what was achieved in previous photoelectron studies, and with mass selectivity. The adiabatic ionization potential of cyanobutadiyne is measured at 85 366 (±40) cm−1. Assignment of the vibrational bands of the four lowest electronic states X+2Π, A+2Π, B+2Σ+, and C+2Π are performed, supported by high level ab initio calculations which are fully detailed in Paper II [B. Gans et al., J. Chem. Phys. 150, 244303 (2019)] and by Franck-Condon simulations. Only vibrational stretching modes are observed in the threshold-photoelectron spectra. The ground state of HC5N+ exhibits a vibrational progression in the ν2 stretching mode involving mainly the elongation of the C≡C triple bonds, whereas the A+ and C+ excited electronic states show a progression in the stretching mode mainly associated with the elongation of the C≡N bond, i.e., ν4 and ν3, respectively. The B+ state appears almost as a vibrationless structure in close vicinity to the A+ state. [ABSTRACT FROM AUTHOR]

In this paper, we improve the collocation method for computing vibrational spectra that was presented in the work of Avila and Carrington, Jr. [J. Chem. Phys. 143, 214108 (2015)]. Known quadrature and collocation methods using a Smolyak grid require storing intermediate vectors with more elements than points on the Smolyak grid. This is due to the fact that grid labels are constrained among themselves and basis labels are constrained among themselves. We show that by using the so-called hierarchical basis functions, one can significantly reduce the memory required. In this paper, the intermediate vectors have only as many elements as the Smolyak grid. The ideas are tested by computing energy levels of CH2NH. [ABSTRACT FROM AUTHOR]

The use of computational optimal transport is investigated as a tool for comparing two molecular spectra. Unlike other techniques for comparing molecular spectra in a pattern-recognition framework, transport distances simultaneously encode information about line positions and intensities. In addition, it is shown that transport distances are a useful alternative to Euclidean distances as Euclidean distances are based on line-by-line comparisons, while transport distances reflect broader features of molecular spectra and adequately compare spectra with different resolutions. This paper includes a tutorial on the use of optimal transport and investigates several well-chosen examples to illustrate the utility of computational optimal transport for comparing molecular spectra. [ABSTRACT FROM AUTHOR]

In this paper, we propose to compute the electronic circular dichroism (ECD) spectra of chiral molecules using a real-time propagation of the time-dependent Schrödinger equation (TDSE) in the space of electronic field-free eigenstates, by coupling TDSE with a given treatment of the electronic structure of the target. The time-dependent induced magnetic moment is used to compute the ECD spectrum from an explicit electric perturbation. The full matrix representing the transition magnetic moment in the space of electronic states is generated from that among pairs of molecular orbitals. In the present work, we show the ECD spectra of methyloxirane, of several conformers of L-alanine, and of the Λ-Co(acac)3 complex, computed from a singly excited ansatz of time-dependent density functional theory eigenstates. The time-domain ECD spectra properly reproduce the frequency-domain ones obtained in the linear-response regime and quantitatively agree with the available experimental data. Moreover, the time-domain approach to ECD allows us to naturally go beyond the ground-state rotationally averaged ECD spectrum, which is the standard outcome of the linear-response theory, e.g., by computing the ECD spectra from electronic excited states. [ABSTRACT FROM AUTHOR]

In this paper, we present a comprehensive account of quantum dissipation theories with the quadratic environment couplings. The theoretical development includes the Brownian solvation mode embedded hierarchical quantum master equations, a core-system hierarchy construction that verifies the extended dissipaton equation of motion (DEOM) formalism [R. X. Xu et al., J. Chem. Phys. 148, 114103 (2018)]. Developed are also the quadratic imaginary-time DEOM for equilibrium and the λ(t)-DEOM for nonequilibrium thermodynamics problems. Both the celebrated Jarzynski equality and Crooks relation are accurately reproduced, which, in turn, confirms the rigorousness of the extended DEOM theories. While the extended DEOM is more numerically efficient, the core-system hierarchy quantum master equation is favorable for "visualizing" the correlated solvation dynamics. [ABSTRACT FROM AUTHOR]

The electrical manipulation of magnetization is appealing to the relevant experiment and spintronic device. In this paper, we focus on the electrical and magnetic properties of a thin film cleaved from the nonmagnetic half-Heusler compound CoTiSb. By means of the first-principles calculations, we find that the surface of TiSb termination possesses ferrimagnetism with a magnetic moment of 0.35 (0.49) μB per unit cell without (with) Hubbard U, which can persist below the Curie temperature of 48 (54) K. What is more, such a surface magnetism can be tuned to nonmagnetism by gate-induced hole doping with a concentration of 2.83 × 1014 (3.55 × 1014) cm−2. This magnetic tunability of the CoTiSb slab provides a platform to realize the pseudo-spin valve with both the magnetic electrodes and nonmagnetic space layer made of the same material without hetero-interfaces. [ABSTRACT FROM AUTHOR]

With an ever-increasing interest in water properties, many intermolecular force fields have been proposed to describe the behavior of water. Unfortunately, good models for liquid water usually cannot provide simultaneously an accurate melting point for ice. For this reason, the TIP4P/Ice model was developed for targeting the melting point and has become the preferred choice for simulating ice at coexistence. Unfortunately, available data for its dynamic properties in the liquid state are scarce. Therefore, we demonstrate a series of simulations aimed at the calculation of transport coefficients for the TIP4P/Ice model over a large range of thermodynamic conditions, ranging from T = 245 K to T = 350 K, for the temperature, and from p = 0 to p = 500 MPa, for the pressure. We have found that the self-diffusion (shear viscosity) exhibits smaller (increased) values than TIP4P/2005 and experiments. However, rescaling the temperature with respect to the triple point temperature, as in a corresponding states plot, we find that TIP4P/Ice compares very well with TIP4P/2005 and experiment. Such observations allow us to infer that despite the different original purposes of these two models examined here, one can benefit from a vast number of reports regarding the behavior of transport coefficients for the TIP4P/2005 model and utilize them following the routine described in this paper. [ABSTRACT FROM AUTHOR]

Superionic (SI) water ices—high-temperature, high-pressure phases of water in which oxygen ions occupy a regular crystal lattice whereas the protons flow in a liquid-like manner—have attracted a growing amount of attention over the past few years, in particular due to their possible role in the magnetic anomalies of the ice giants Neptune and Uranus. In this paper, we consider the calculation of the free energies of such phases, exploring hybrid reference systems consisting of a combination of an Einstein solid for the oxygen ions occupying a crystal lattice and a Uhlenbeck-Ford potential for the protonic fluid that avoids irregularities associated with possible particle overlaps. Applying this approach to a recent neural-network potential-energy landscape for SI water ice, we compute Gibbs free energies as a function of temperature for the SI fcc and liquid phases to determine the melting temperature Tm at 340 GPa. The results are consistent with previous estimates and indicate that the entropy difference between both phases is comparatively small, in particular due to the large amplitude of vibration of the oxygen ions in the fcc phase at the melting temperature. [ABSTRACT FROM AUTHOR]

Utilizing localized orbitals, local correlation theory can reduce the unphysically high system-size scaling of post-Hartree–Fock (post-HF) methods to linear scaling in insulating molecules. The sparsity of the four-index electron repulsion integral (ERI) tensor is central to achieving this reduction. For second-order Møller–Plesset theory (MP2), one of the simplest post-HF methods, only the (ia|jb) ERIs are needed, coupling occupied orbitals i, j and virtuals a, b. In this paper, we compare the numerical sparsity (called the "ragged list") and two other approaches revealing the low-rank sparsity of the ERI. The ragged list requires only one set of (localized) virtual orbitals, and we find that the orthogonal valence virtual-hard virtual set of virtuals originally proposed by Subotnik et al. gives the sparsest ERI tensor. To further compress the ERI tensor, the pair natural orbital (PNO) type representation uses different sets of virtual orbitals for different occupied orbital pairs, while the occupied-specific virtual (OSV) approach uses different virtuals for each occupied orbital. Our results indicate that while the low-rank PNO representation achieves significant rank reduction, it also requires more memory than the ragged list. The OSV approach requires similar memory to that of the ragged list, but it involves greater algorithmic complexity. An approximation (called the "fixed sparsity pattern") for solving the local MP2 equations using the numerically sparse ERI tensor is proposed and tested to be sufficiently accurate and to have highly controllable error. A low-scaling local MP2 algorithm based on the ragged list and the fixed sparsity pattern is therefore promising. [ABSTRACT FROM AUTHOR]

The excitation energy transfer (EET) process for photosynthetic antenna complexes consisting of subunits, each comprised of multiple chromophores, remains challenging to describe. The multichromophoric Förster resonance energy transfer theory is a popular method to describe the EET process in such systems. This paper presents a new time-domain method for calculating energy transfer based on the combination of multichromophoric Förster resonance energy transfer theory and the Numerical Integration of the Schrödinger Equation method. After validating the method on simple model systems, we apply it to the Light-Harvesting antenna 2 (LH2) complex, a light harvesting antenna found in purple bacteria. We use a simple model combining the overdamped Brownian oscillators to describe the dynamic disorder originating from the environmental fluctuations and the transition charge from the electrostatic potential coupling model to determine the interactions between chromophores. We demonstrate that with this model, both the calculated spectra and the EET rates between the two rings within the LH2 complex agree well with experimental results. We further find that the transfer between the strongly coupled rings of neighboring LH2 complexes can also be well described with our method. We conclude that our new method accurately describes the EET rate for biologically relevant multichromophoric systems, which are similar to the LH2 complex. Computationally, the new method is very tractable, especially for slow processes. We foresee that the method can be applied to efficiently calculate transfer in artificial systems as well and may pave the way for calculating multidimensional spectra of extensive multichromophoric systems in the future. [ABSTRACT FROM AUTHOR]

Over the last century, quantum theories have revolutionized our understanding of material properties. One of the most striking quantum phenomena occurring in heterogeneous media is the quantum tunneling effect, where carriers can tunnel through potential barriers even if the barrier height exceeds the carrier energy. Interestingly, the tunneling process can be accompanied by the absorption or emission of light. In most tunneling junctions made of noble metal electrodes, these optical phenomena are governed by plasmonic modes, i.e., light-driven collective oscillations of surface electrons. In the emission process, plasmon excitation via inelastic tunneling electrons can improve the efficiency of photon generation, resulting in bright nanoscale optical sources. On the other hand, the incident light can affect the tunneling behavior of plasmonic junctions as well, leading to phenomena such as optical rectification and induced photocurrent. Thus, plasmonic tunneling junctions provide a rich platform for investigating light–matter interactions, paving the way for various applications, including nanoscale light sources, sensors, and chemical reactors. In this paper, we will introduce recent research progress and promising applications based on plasmonic tunneling junctions. [ABSTRACT FROM AUTHOR]

Many phase change materials (PCMs) are found to crystallize without exhibiting a glass transition endotherm upon reheating. In this paper, we review experimental evidence revealing that these PCMs and likely other hyperquenched molecular and metallic systems can crystallize from the glassy state when reheated at a standard rate. Among these evidences, PCMs annealed below the glass transition temperature Tg exhibit slower crystallization kinetics despite an increase in the number of sub-critical nuclei that should promote the crystallization speed. Flash calorimetry uncovers the glass transition endotherm hidden by crystallization and reveals a distinct change in kinetics when crystallization switches from the glassy to the supercooled liquid state. The resulting Tg value also rationalizes the presence of the pre-Tg relaxation exotherm ubiquitous of hyperquenched systems. Finally, the shift in crystallization temperature during annealing exhibits a non-exponential decay that is characteristic of structural relaxation in the glass. Modeling using a modified Turnbull equation for nucleation rate supports the existence of sub-Tg fast crystallization and emphasizes the benefit of a fragile-to-strong transition for PCM applications due to a reduction in crystallization at low temperature (improved data retention) and increasing its speed at high temperature (faster computing). [ABSTRACT FROM AUTHOR]

Introducing nanorods into a polymeric matrix can enhance the physical and mechanical properties of the resulting material. In this paper, we focus on understanding the dispersion and orientation patterns of nanorods in an unentangled polymer melt, particularly as a function of nanorod concentration, using molecular dynamics simulations. The system is comprised of flexible polymer chains and multi-thread nanorods that are equilibrated in the NPT ensemble. All interactions are purely repulsive except for those between polymers and rods. Results with attractive vs repulsive polymer–rod interactions are compared and contrasted. The concentration of rods has a direct impact on the phase behavior of the system. At lower concentrations, rods phase separate into nematic clusters, whereas at higher concentrations more isotropic and less structured rod configurations are observed. A detailed examination of the conformation of the polymer chains near the rod surface shows extension of the chains along the director of the rods (especially within clusters). The dispersion and orientation of the nanorods are a result of the competition between depletion entropic forces responsible for the formation of rod clusters, the enthalpic effects that improve mixing of rods and polymer, and entropic losses of polymers interpenetrating rod clusters. [ABSTRACT FROM AUTHOR]

Electrochemical impedance spectroscopy (EIS) is a powerful characterization technique for the in-depth investigation of kinetic/transport parameters detection, reaction mechanism understanding, and degradation effects exploration in lithium-ion battery (LIB) systems. However, due to the lack of standardized criterion/paradigm, severe misinterpretations occur frequently during an EIS measurement. In this paper, the significance of instrumental accuracy is described and the character/principle of selection on the simulation model is illuminated/proposed, showing that an adequate precision device and an appropriate fitting model are a prerequisite for a correct EIS analysis. Moreover, the drawbacks of conventional two-electrode EIS experiments for typical coin-type cells are rigorously pointed out by comparison with the ideal three-electrode configuration, where the real impedance information of the cathode would be masked by the sum of both the anode film resistance response and the unavoidable inductive loop signal. The three-electrode case enables efficient accurate observations on individual electrodes, thus facilitating abundant and useful information acquisition. Consequently, devices with a sufficient accuracy, rational simulation models, and advanced three-electrode cells are distinctly illustrated as standardized criterion/paradigm for EIS characterizations, which are essentially important for electrode and interface modifications in LIBs. [ABSTRACT FROM AUTHOR]

Fewest-switches surface hopping (FSSH) has emerged as one of the leading methods for modeling the quantum dynamics of molecular systems. While its original formulation was limited to adiabatic populations, the growing interest in the application of FSSH to coherent phenomena prompts the question of how one should construct a complete density matrix based on FSSH trajectories. A straightforward solution is to define adiabatic coherences based on wavefunction coefficients. In this paper, we demonstrate that inconsistencies introduced in the density matrix through such treatment may lead to a violation of positivity. We furthermore show that a recently proposed coherent generalization of FSSH results in density matrices that satisfy positivity while yielding improved accuracy throughout much (but not all) of parameter space. [ABSTRACT FROM AUTHOR]

The transport of particles through channels is of paramount importance in physics, chemistry, and surface science due to its broad real world applications. Much insight can be gained by observing the transition paths of a particle through a channel and collecting statistics on the lifetimes in the channel or the escape probabilities from the channel. In this paper, we consider the diffusive transport through a narrow conical channel of a Brownian particle subject to intermittent dynamics, namely, stochastic resetting. As such, resetting brings the particle back to a desired location from where it resumes its diffusive phase. To this end, we extend the Fick–Jacobs theory of channel-facilitated diffusive transport to resetting-induced transport. Exact expressions for the conditional mean first passage times, escape probabilities, and the total average lifetime in the channel are obtained, and their behavior as a function of the resetting rate is highlighted. It is shown that resetting can expedite the transport through the channel—rigorous constraints for such conditions are then illustrated. Furthermore, we observe that a carefully chosen resetting rate can render the average lifetime of the particle inside the channel minimal. Interestingly, the optimal rate undergoes continuous and discontinuous transitions as some relevant system parameters are varied. The validity of our one-dimensional analysis and the corresponding theoretical predictions is supported by three-dimensional Brownian dynamics simulations. We thus believe that resetting can be useful to facilitate particle transport across biological membranes—a phenomenon that can spearhead further theoretical and experimental studies. [ABSTRACT FROM AUTHOR]

Using an approach based on the diffusion analog of the Cattaneo–Vernotte differential model, we find the exact analytical solution to the corresponding time-dependent linear hyperbolic initial boundary value problem, describing irreversible diffusion-controlled reactions under Smoluchowski's boundary condition on a spherical sink. By means of this solution, we extend exact analytical calculations for the time-dependent classical Smoluchowski rate coefficient to the case that includes the so-called inertial effects, occurring in the host media with finite relaxation times. We also present a brief survey of Smoluchowski's theory and its various subsequent refinements, including works devoted to the description of the short-time behavior of Brownian particles. In this paper, we managed to show that a known Rice's formula, commonly recognized earlier as an exact reaction rate coefficient for the case of hyperbolic diffusion, turned out to be only its approximation being a uniform upper bound of the exact value. Here, the obtained formula seems to be of great significance for bridging a known gap between an analytically estimated rate coefficient on the one hand and molecular dynamics simulations together with experimentally observed results for the short times regime on the other hand. A particular emphasis has been placed on the rigorous mathematical treatment and important properties of the relevant initial boundary value problems in parabolic and hyperbolic diffusion theories. [ABSTRACT FROM AUTHOR]

In this paper, we investigate prenucleation and heterogeneous nucleation in the liquid Pb/solid Al system as an example of systems with large lattice misfit using molecular dynamics simulation. Solid Pb and Al have a large positive lattice misfit (f) of 18.2% along the densely packed [110] direction. This study reveals that prenucleation occurs at 600 K (an undercooling of 15 K), and a 2-dimensional (2D) ordered structure forms at the interface with a coincidence site lattice (CSL) between the first Pb and first Al layers. The CSL accommodates the major part of the f, and only a small residual lattice misfit (fr) of 1.9% remains. The formation of the CSL transforms the original substrate into a considerably potent nucleant, where the first Pb layer becomes the new surface layer of the substrate. At an undercooling of about 22 K, nucleation proceeds by merging 2D ordered structure through structural templating: the second Pb layer is epitaxial to the CSL Pb layer, the third Pb layer largely accommodates the fr, and the fourth Pb layer is a nearly perfect crystalline plane. Further analysis indicates that the interface with the CSL has a lower interfacial energy than with a cube-to-cube orientation relationship. For the first time, we established that the CSL was an effective mechanism to accommodate the f for systems with a large positive misfits. Heterogeneous nucleation is governed not by a single mechanism (misfit dislocations in Turnbull's model), but instead by various mechanisms depending on f. This study sheds new light on the atomistic mechanism of heterogeneous nucleation. [ABSTRACT FROM AUTHOR]

Coarse-grained (CG) models facilitate an efficient exploration of complex systems by reducing the unnecessary degrees of freedom of the fine-grained (FG) system while recapitulating major structural correlations. Unlike structural properties, assessing dynamic properties in CG modeling is often unfeasible due to the accelerated dynamics of the CG models, which allows for more efficient structural sampling. Therefore, the ultimate goal of the present series of articles is to establish a better correspondence between the FG and CG dynamics. To assess and compare dynamical properties in the FG and the corresponding CG models, we utilize the excess entropy scaling relationship. For Paper I of this series, we provide evidence that the FG and the corresponding CG counterpart follow the same universal scaling relationship. By carefully reviewing and examining the literature, we develop a new theory to calculate excess entropies for the FG and CG systems while accounting for entropy representability. We demonstrate that the excess entropy scaling idea can be readily applied to liquid water and methanol systems at both the FG and CG resolutions. For both liquids, we reveal that the scaling exponents remain unchanged from the coarse-graining process, indicating that the scaling behavior is universal for the same underlying molecular systems. Combining this finding with the concept of mapping entropy in CG models, we show that the missing entropy plays an important role in accelerating the CG dynamics. [ABSTRACT FROM AUTHOR]

For many laser-based methods of material characterization and analysis, a tensor formulation of theory is necessary, especially in techniques that exploit nonlinear or chiral optics. The fundamental interactions that underpin such methods offer various levels of approach to theory, but the most rigorous often lead to equations of considerable complexity. To compute the values for individual material parameters frequently demands making assumptions of extreme simplicity, overly dependent on calculational method, yet still providing unsatisfactory results. A pragmatic and entirely rigorous symmetry-based approach to the irreducible tensorial structures circumvents many of these problems, securing reliable results and guiding the pathway to applications. Instead of focusing on individual tensor components, such an approach can rapidly determine the number of linearly independent quantities—and hence the number of operationally different setups necessary for full characterization. By such means, one can directly ascertain how variations of optical polarization and beam geometry can reliably capture the response of any material system. The use of an irreducible tensor method operates independently of any means that might be chosen to calculate material properties. It removes the need for common simplifying assumptions, such as the approximation of tensorial structure by a scalar representation, adoption of a two-state model, or disregarding near-resonance damping. It also obviates any dependence on a choice of simulation package or quantum-calculational software. In this paper, the principles are set down and illustrated by application to experiments of varying degrees of complexity, including interactions of growing significance in the realm of chiral nonlinear optics. Limitations of this approach are also critically assessed. [ABSTRACT FROM AUTHOR]