{"title":"Simplification of solvation shell with water clusters in the simulation of electrochemical nitrogen reduction reaction.","authors":"Han-Xuan Liu, Zheng-Qing Huang, Xin Gao, Guidong Yang, Chun-Ran Chang","doi":"10.1063/5.0230137","DOIUrl":"https://doi.org/10.1063/5.0230137","url":null,"abstract":"<p><p>Electrochemical methods for nitrogen reduction have received extensive attention due to the mild reaction conditions. In order to gain an insight into the mechanism of the electrochemical nitrogen reduction process, theoretical simulations are necessary. However, current simulation studies contain many imprecise approximations that may hinder the real recognition of the reaction process. Although solvation methods have recently been developed to provide efficient descriptions, their further applications are hindered by controversy over modeling of solvents and the enormous computational effort required. In this work, we simplify the solvation conditions by using water clusters and compare them with an accurate water layer model. The results demonstrate that the simplified water clusters can effectively capture protons, simulate the surface electric environment, and enable the calculation of the activation energy of the reaction. This method offers an affordable approach for simulating the surface potentials and solvents and provides a new reference for the theoretical study of the electrochemical nitrogen reduction reaction.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 18","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142620932","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Felix Kahlert, Valentin Link, Richard Hartmann, Walter T Strunz
{"title":"Simulating the Landau-Zener sweep in deeply sub-Ohmic environments.","authors":"Felix Kahlert, Valentin Link, Richard Hartmann, Walter T Strunz","doi":"10.1063/5.0235741","DOIUrl":"https://doi.org/10.1063/5.0235741","url":null,"abstract":"<p><p>With the goal to study dissipative Landau-Zener (LZ) sweeps in realistic solid-state qubits, we utilize novel methods from non-Markovian open quantum system dynamics that enable reliable long-time simulations for sub-Ohmic environments. In particular, we combine a novel representation of the dynamical propagator, the uniform time evolving matrix product operator method, with a stochastic realization of finite temperature fluctuations. The latter greatly reduces the computational cost for the matrix product operator approach, enabling convergence in the experimentally relevant deeply sub-Ohmic regime. Our method allows the exact simulation of dynamical protocols with long operation times, such as the LZ sweep, in challenging parameter regimes that are realized in current experimental platforms.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 18","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142620935","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A computer simulation study of a chiral active ring polymer.","authors":"Shalabh K Anand","doi":"10.1063/5.0232538","DOIUrl":"https://doi.org/10.1063/5.0232538","url":null,"abstract":"<p><p>We investigate a ring polymer under the influence of chiral active Brownian forces in two dimensions using coarse-grained computer simulations. We observe a non-monotonic behavior of the radius of gyration of an active Brownian ring as a function of active force. However, the shrinkage of the ring in the intermediate strength of active forces becomes more pronounced in the presence of chiral active forces, and the shrinkage is monotonic at a given activity level as a function of the angular frequency controlling the direction of the active force. The distribution of radius of gyration, inter-monomer distance, and radial distribution suggest that the monomers come close to each other, eventually leading to the shrinkage of the ring. Moreover, the bond-correlation suggests that the chirality introduces a local folding of the monomers. Furthermore, using the diameter correlation function, we show that the ring performs tank-treading motion with a frequency following power-law relation with active force with exponent 3/2. The mean squared displacement of the monomers further assists the tank-treading dynamics by exhibiting oscillatory behavior.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 18","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142604650","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Hydrogen diffusion on Ni(100): A combined machine-learning, ring polymer molecular dynamics, and kinetic Monte Carlo study.","authors":"J Steffen, A Alibakhshi","doi":"10.1063/5.0226314","DOIUrl":"https://doi.org/10.1063/5.0226314","url":null,"abstract":"<p><p>We introduce a methodological framework coupling machine-learning potentials, ring polymer molecular dynamics (RPMD), and kinetic Monte Carlo (kMC) to draw a comprehensive physical picture of the collective diffusion of hydrogen atoms on metal surfaces. For the benchmark case of hydrogen diffusion on a Ni(100) surface, the hydrogen adsorption and diffusion energetics and its dependence on the local coverage is described via a neural-network potential, where the training data are computed via periodic density functional theory (DFT) and include all relevant optimized diffusion and desorption paths, sampled by nudged elastic band optimizations and molecular dynamics simulations. Nuclear quantum effects, being crucial for processes involving hydrogen at low temperatures, are treated by RPMD. The diffusion rate constants are calculated with a combination of umbrella samplings employed to map the free energy profile and separate samplings of recrossing trajectories to obtain the transmission coefficient. The calculated diffusion rates for different temperatures and local environments are then combined and fitted into a kMC model allowing access to larger time and length scales. Our results demonstrate an outstanding performance for the trained neural network potential in reproducing reference DFT energies and forces. We report the effective diffusion rates for different temperatures and hydrogen surface coverages obtained via this recipe in good agreement with the experimental results. The method combination proposed in this study can be instrumental for a wide range of applications in materials science.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 18","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142621407","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Influence of temperature and crack-tip speed on crack propagation in elastic solids.","authors":"B N J Persson","doi":"10.1063/5.0236699","DOIUrl":"https://doi.org/10.1063/5.0236699","url":null,"abstract":"<p><p>I study the influence of temperature and the crack-tip velocity of bond breaking at the crack tip in rubber-like materials. Bond breaking is considered as a stress-aided thermally activated process and results in an effective crack propagation energy, which increases strongly with decreasing temperature or increasing crack-tip speed. This effect is particularly important for adhesive (interfacial) crack propagation but less important for cohesive (bulk) crack propagation owing to the much larger bond-breaking energies in the latter case. For adhesive cracks, the theory results are consistent with adhesion measurements for silicone rubber polydimethylsiloxane (PDMS) in contact with silica glass surfaces. For cohesive cracks, the theory agrees well with experimental results PDMS films chemically bound to silanized glass.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 18","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142621409","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Propagation gap for shear waves in binary liquids: Analytical and simulation study.","authors":"Taras Bryk, Maria Kopcha, Ihor Yidak","doi":"10.1063/5.0236047","DOIUrl":"https://doi.org/10.1063/5.0236047","url":null,"abstract":"<p><p>Transverse collective excitations in 50-50 and 80-20 Lennard-Jones binary liquid mixtures are studied for different mass ratio of components R at fixed numerical densities. Increasing the mass ratio results in a growing difference between frequencies of shear waves and transverse optic modes. We report an increase in the propagation gap width for shear waves with mass ratio of components and compare it to the gap width expression, known from the transverse dynamics of simple liquids. A four-variable dynamic model of transverse dynamics in binary liquids with an account of cross correlations between total-mass and mass-concentration transverse current fluctuations is solved analytically in the long-wavelength region. An equation for the propagation gap of shear waves for binary liquids is reported and analyzed.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 18","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142620485","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Reduced density matrix dynamics in multistate harmonic models via time-convolution and time-convolutionless quantum master equations with quantum-mechanical and semiclassical kernels.","authors":"Xiang Sun, Zengkui Liu","doi":"10.1063/5.0231561","DOIUrl":"https://doi.org/10.1063/5.0231561","url":null,"abstract":"<p><p>In this work, we explore the electronic reduced density matrix (RDM) dynamics using time-convolution (TC) and time-convolutionless (TCL) quantum master equations (QMEs) that are based on perturbative electronic couplings within the framework of multistate harmonic (MSH) models. The MSH model Hamiltonian consistently incorporates the electronic-vibrational correlations between all pairs of states by satisfying the pairwise reorganization energies directly obtained from all-atom simulations, representing the globally heterogeneous environments that couple to the multiple states differently. We derive the exact quantum-mechanical and a hierarchy of semiclassical approximate expressions for the kernels in TC and TCL QMEs that project the full RDM for general shifted harmonic systems, including the MSH model. These QMEs are applied to simulate RDM dynamics of photoinduced charge transfer (PICT) in organic photovoltaic carotenoid-porphyrin-fullerene triad solvated in tetrahydrofuran solution and the excitation energy transfer (EET) dynamics in photosynthetic Fenna-Matthews-Olson complexes from C. tepidum and P. aestuarii. Our results show that while both TC and TCL QMEs capture similar phenomena in PICT and EET processes, TC QME generally provides more accurate results than TCL QME, particularly in the initial oscillation of EET population dynamics. This study highlights the effectiveness of the TC and TCL QMEs in modeling RDM dynamics of nonadiabatic processes, offering insights for realistic condensed phase systems.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 18","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142620614","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Kai Töpfer, Shyamsunder Erramilli, Lawrence D Ziegler, Markus Meuwly
{"title":"Energy relaxation of N2O in gaseous, supercritical, and liquid xenon and SF6.","authors":"Kai Töpfer, Shyamsunder Erramilli, Lawrence D Ziegler, Markus Meuwly","doi":"10.1063/5.0235760","DOIUrl":"https://doi.org/10.1063/5.0235760","url":null,"abstract":"<p><p>Rotational and vibrational energy relaxation (RER and VER) of N2O embedded in xenon and SF6 environments ranging from the gas phase to the liquid, including the supercritical regime, is studied at a molecular level. Calibrated intermolecular interactions from high-level electronic structure calculations, validated against experiments for the pure solvents, were used to carry out classical molecular dynamics simulations corresponding to experimental state points for near-critical isotherms. The computed RER rates in low-density solvents of krotXe=(3.67±0.25)×1010 s-1 M-1 and krotSF6=(1.25±0.12)×1011 s-1 M-1 compare well with the rates determined by the analysis of two-dimensional infrared experiments. Simulations find that an isolated binary collision description is successful up to solvent concentrations of ∼4 M. For higher densities, including the supercritical regime, the simulations do not correctly describe RER, probably due to the neglect of solvent-solute coupling in the analysis of the rotational motion. For VER, the near-quantitative agreement between simulations and pump-probe experiments captures the solvent density-dependent trends.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 18","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142621405","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Improvement on the screening of nonlinear commutator operations in selective coupled-cluster using Lagrangian.","authors":"Suhita Basumallick, Enhua Xu, Seiichiro L Ten-No","doi":"10.1063/5.0231739","DOIUrl":"https://doi.org/10.1063/5.0231739","url":null,"abstract":"<p><p>We introduce a Lagrangian implementation of the full coupled-cluster reduction [Xu et al., Phys. Rev. Lett. 121, 113001 (2018)], that is, a selected coupled-cluster (CC) based on an arbitrary-order full CC expansion using direct commutator expansions. In this method, the screening for the products of cluster amplitudes plays a central role to reduce the computational cost for the nonlinear commutator operations, while the convergence of the total energy in the standard energy expression is not rapid with tightening the threshold. The new implementation using Lagrangian is robust, containing error only quadratic to those of amplitudes, allowing a much larger screening threshold. We demonstrate the performance of the new implementation by investigating the calculations of N2 and C6H6. The accuracy and applicability are also demonstrated for the potential energy curve of H2O in comparison with conventional quantum chemical methods.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 18","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142621408","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Learning glass transition temperatures via dimensionality reduction with data from computer simulations: Polymers as the pilot case.","authors":"Artem Glova, Mikko Karttunen","doi":"10.1063/5.0229161","DOIUrl":"https://doi.org/10.1063/5.0229161","url":null,"abstract":"<p><p>Machine learning methods provide an advanced means for understanding inherent patterns within large and complex datasets. Here, we employ the principal component analysis (PCA) and the diffusion map (DM) techniques to evaluate the glass transition temperature (Tg) from low-dimensional representations of all-atom molecular dynamic simulations of polylactide (PLA) and poly(3-hydroxybutyrate) (PHB). Four molecular descriptors were considered: radial distribution functions (RDFs), mean square displacements (MSDs), relative square displacements (RSDs), and dihedral angles (DAs). By applying Gaussian Mixture Models (GMMs) to analyze the PCA and DM projections and by quantifying their log-likelihoods as a density-based metric, a distinct separation into two populations corresponding to melt and glass states was revealed. This separation enabled the Tg evaluation from a cooling-induced sharp increase in the overlap between log-likelihood distributions at different temperatures. Tg values derived from the RDF and MSD descriptors using DM closely matched the standard computer simulation-based dilatometric and dynamic Tg values for both PLA and PHB models. This was not the case for PCA. The DM-transformed DA and RSD data resulted in Tg values in agreement with experimental ones. Overall, the fusion of atomistic simulations and DMs complemented with the GMMs presents a promising framework for computing Tg and studying the glass transition in a unified way across various molecular descriptors for glass-forming materials.</p>","PeriodicalId":15313,"journal":{"name":"Journal of Chemical Physics","volume":"161 18","pages":""},"PeriodicalIF":3.1,"publicationDate":"2024-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142604685","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}