The growth of stringlets upon heating ultimately also leads to the "softening" of these excitations, and the boson peak frequency and shear modulus drop in concert with this softening. The growth of string-like collective motion upon heating in the fast-dynamics regime is further shown to be responsible for the growth in the intensity of the fast relaxation process. Relaxation in cooled liquids clearly involves a hierarchy of relaxation processes acting on rather different timescales and spatial scales.We report on a theoretical study of second-harmonic generation (SHG) in plasmonic nanostructures interacting with two-level quantum emitters (QEs) under incoherent energy pump. We generalize the driven-dissipative Tavis-Cummings model by introducing the anharmonic surface plasmon-polariton (SPP) mode coupled to QEs and examine physical properties of corresponding SPP-QE polariton states. https://www.selleckchem.com/products/md-224.html Our calculations of the SHG efficiency for strong QE-SPP coupling demonstrate orders of magnitude enhancement facilitated by the polariton gain. We further discuss time-domain numerical simulations of SHG in a square lattice comprising Ag nanopillars coupled to QEs utilizing a fully vectorial nonperturbative nonlinear hydrodynamic model for conduction electrons coupled to Maxwell-Bloch equations for QEs. The simulations support the idea of gain enhanced SHG and show orders of magnitude increase in the SHG efficiency as the QEs are tuned in resonance with the lattice plasmon mode and brought above the population inversion threshold by incoherent pumping. By varying pump frequency and tuning QEs to a localized plasmon mode, we demonstrate further enhancement of the SHG efficiency facilitated by strong local electric fields. The incident light polarization dependence of the SHG is examined and related to the symmetries of participating plasmon modes.In this study, we investigate the nuclear quantum effects (NQEs) on the acidity constant (pKA) of liquid water isotopologs under the ambient condition by path integral molecular dynamics (PIMD) simulations. We compared simulations using a fully explicit solvent model with a classical polarizable force field, density functional tight binding, and ab initio density functional theory, which correspond to empirical, semiempirical, and ab initio PIMD simulations, respectively. The centroid variable with respect to the proton coordination number of a water molecule was restrained to compute the gradient of the free energy, which measures the reversible work of the proton abstraction for the quantum mechanical system. The free energy curve obtained by thermodynamic integration was used to compute the pKA value based on probabilistic determination. This technique not only reproduces the pKA value of liquid D2O experimentally measured (14.86) but also allows for a theoretical prediction of the pKA values of liquid T2O and aqueous HDO and HTO, which are unknown due to their scarcity. It is also shown that the NQEs on the free energy curve can result in a downshift of 4.5 ± 0.9 pKA units in the case of liquid water, which indicates that the NQEs plays an indispensable role in the absolute determination of pKA. The results of this study can help inform further extensions into the calculation of the acidity constants of isotope substituted species with high accuracy.We performed ab initio molecular dynamics (AIMD) simulations to benchmark bulk liquid structures and to evaluate results from all-atom force field molecular dynamics (FFMD) simulations with the generalized Amber force field (GAFF) for organophosphorus (OP) and organochlorine (OC) compounds. Our work also addresses the current and important topic of force field validation, applied here to a set of nonaqueous organic liquids. Our approach differs from standard treatments, which validate force fields based on thermodynamic data. Utilizing radial distribution functions (RDFs), our results show that GAFF reproduces the AIMD-predicted asymmetric liquid structures moderately well for OP compounds that contain bulky alkyl groups. Among the OCs, RDFs obtained from FFMD overlap well with AIMD results, with some offsets in position and peak structuring. However, re-parameterization of GAFF for some OCs is needed to reproduce fully the liquid structures predicted by AIMD. The offsets between AIMD and FFMD peak positions suggest inconsistencies in the developed force fields, but, in general, GAFF is able to capture short-ranged and long-ranged interactions of OPs and OCs observed in AIMD. Along with the local coordination structure, we also compared enthalpies of vaporization. Overall, calculated bulk properties from FFMD compared reasonably well with experimental values, suggesting that small improvements within the FF should focus on parameters that adjust the bulk liquid structures of these compounds.The combination of Markov state modeling (MSM) and molecular dynamics (MD) simulations has been shown in recent years to be a valuable approach to unravel the slow processes of molecular systems with increasing complexity. While the algorithms for intermediate steps in the MSM workflow such as featurization and dimensionality reduction have been specifically adapted to MD datasets, conventional clustering methods are generally applied to the discretization step. This work adds to recent efforts to develop specialized density-based clustering algorithms for the Boltzmann-weighted data from MD simulations. We introduce the volume-scaled common nearest neighbor (vs-CNN) clustering that is an adapted version of the common nearest neighbor (CNN) algorithm. A major advantage of the proposed algorithm is that the introduced density-based criterion directly links to a free-energy notion via Boltzmann inversion. Such a free-energy perspective allows a straightforward hierarchical scheme to identify conformational clusters at different levels of a generally rugged free-energy landscape of complex molecular systems.The present work introduces a new form of explicitly correlated factor in the context of the transcorrelated methods. The new correlation factor is obtained from the r12 ≈ 0 mathematical analysis of the transcorrelated Hamiltonian, and its analytical form is obtained such that the leading order in 1/r12 of the scalar part of the effective two-electron potential reproduces the long-range interaction of the range-separated density functional theory. The resulting correlation factor exactly imposes the cusp and is tuned by a unique parameter μ, which controls both the depth of the coulomb hole and its typical range in r12. The transcorrelated Hamiltonian obtained with such a new correlation factor has a straightforward analytical expression depending on the same parameter μ, and its physical contents continuously change by varying μ One can change from a non-divergent repulsive Hamiltonian at large μ to a purely attractive one at small μ. We investigate the convergence of the ground state eigenvalues and right eigenvectors of such a new transcorrelated Hamiltonian as a function of the basis set and as a function of μ on a series of two-electron systems.