Search

Results on single-wavelength transient hole burning (SW-THB) developed in paper II [J. Chem. Phys. 115, 4223 (2001)] are applied to the dichroism experiments on anthracene in benzyl alcohol reported in paper I [J. Chem. Phys. 115, 4212 (2001)]. The intermediate component of the dichroism decay is assigned to a SW-THB effect caused by nonpolar electronic solvation. The presence of a solvation component in dichroism experiments has not been demonstrated previously. The sparseness of anthracene’s electronic spectrum eliminates vibrational dynamics from the solvation measurement. Because data collection is focused on a single dimension, the viscosity dependence of the nonpolar solvation is determined with greater accuracy than in our previous two-dimensional transient hole-burning studies. The solvation time is obtained as a function of viscosity/temperature from 14.4 to 2.7 cP (1–56 °C). The times show good agreement with a viscoelastic theory of the diffusive component of nonpolar solvation. Combining the results of this paper with those of paper I allows for comparison of solvation and rotation dynamics within a single system. A correlation between the ratio of diffusive solvation and rotation times and the magnitude of the inertial rotation is suggested. © 2001 American Institute of Physics. [ABSTRACT FROM AUTHOR]

A complete set of four polarization-dependent differential cross sections in the reactions of Cl + aligned-CHD3(v1 = 1, |jK) → DCl(v = 0) + CHD2(v1 = 1) is reported here for two different, rotationally polarized states with j = 1: specifically the |jK = |10 state prepared via the R(0) excitation and the |1 ± 1 state via Q(1). In stark contrast to the complicated situation of the HCl(v) + CD3(v = 0) channel reported in Paper-I, the stereo-requirement of this isotopic channel for both polarized reactants appears quite straightforward and consistent with a direct rebound mechanism. The extent of steric effects is moderate and relatively smaller than the alternative H-atom abstraction channel. All major findings reported here can qualitatively be understood by first noting that the present reaction invokes abstracting a D-atom, which is the spectator in the IR-excitation process. Next, it is recognized that the directional properties of two polarized states of CHD3(v1 = 1, |jK) should manifest primarily in the IR-excited C-H bond, leaving secondary imprints in the unexcited CD3-moiety. The stereo-specificity of the DCl + CHD2 product channel is further reduced by the fact that the abstraction can occur with any one of the three spatially distinct D-atoms. [ABSTRACT FROM AUTHOR]

Erratum: "Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM+][OAc-]" [J. Chem. Phys. We have recomputed the PMFs with AIMD, utilizing a consistent collective variable definition of Eq. (1). This was done using the plumed/ASE interface with a CP2K calculator.[[2], [5]] Corrected PMFs that employ the collective variable of Eq. (1) for all calculations are shown in the corrected Fig. [Extracted from the article]

In Paper I [Litman et al., J. Chem. Phys. (in press) (2022)], we presented the ring-polymer instanton with explicit friction (RPI-EF) method and showed how it can be connected to the ab initio electronic friction formalism. This framework allows for the calculation of tunneling reaction rates that incorporate the quantum nature of the nuclei and certain types of non-adiabatic effects (NAEs) present in metals. In this paper, we analyze the performance of RPI-EF on model potentials and apply it to realistic systems. For a 1D double-well model, we benchmark the method against numerically exact results obtained from multi-layer multi-configuration time-dependent Hartree calculations. We demonstrate that RPI-EF is accurate for medium and high friction strengths and less accurate for extremely low friction values. We also show quantitatively how the inclusion of NAEs lowers the crossover temperature into the deep tunneling regime, reduces the tunneling rates, and, in certain regimes, steers the quantum dynamics by modifying the tunneling pathways. As a showcase of the efficiency of this method, we present a study of hydrogen and deuterium hopping between neighboring interstitial sites in selected bulk metals. The results show that multidimensional vibrational coupling and nuclear quantum effects have a larger impact than NAEs on the tunneling rates of diffusion in metals. Together with Paper I [Litman et al., J. Chem. Phys. (in press) (2022)], these results advance the calculations of dissipative tunneling rates from first principles. [ABSTRACT FROM AUTHOR]

A Brownian bridge is a continuous random walk conditioned to end in a given region by adding an effective drift to guide paths toward the desired region of phase space. This idea has many applications in chemical science where one wants to control the endpoint of a stochastic process—e.g., polymer physics, chemical reaction pathways, heat/mass transfer, and Brownian dynamics simulations. Despite its broad applicability, the biggest limitation of the Brownian bridge technique is that it is often difficult to determine the effective drift as it comes from a solution of a Backward Fokker–Planck (BFP) equation that is infeasible to compute for complex or high-dimensional systems. This paper introduces a fast approximation method to generate a Brownian bridge process without solving the BFP equation explicitly. Specifically, this paper uses the asymptotic properties of the BFP equation to generate an approximate drift and determine ways to correct (i.e., re-weight) any errors incurred from this approximation. Because such a procedure avoids the solution of the BFP equation, we show that it drastically accelerates the generation of conditioned random walks. We also show that this approach offers reasonable improvement compared to other sampling approaches using simple bias potentials. [ABSTRACT FROM AUTHOR]

Performing full-resolution atomistic simulations of nucleic acid folding has remained a challenge for biomolecular modeling. Understanding how nucleic acids fold and how they transition between different folded structures as they unfold and refold has important implications for biology. This paper reports a theoretical model and computer simulation of the ab initio folding of DNA inverted repeat sequences. The formulation is based on an all-atom conformational model of the sugar-phosphate backbone via chain closure, and it incorporates three major molecular-level driving forces—base stacking, counterion-induced backbone self-interactions, and base pairing—via separate analytical theories designed to capture and reproduce the effects of the solvent without requiring explicit water and ions in the simulation. To accelerate computational throughput, a mixed numerical/analytical algorithm for the calculation of the backbone conformational volume is incorporated into the Monte Carlo simulation, and special stochastic sampling techniques were employed to achieve the computational efficiency needed to fold nucleic acids from scratch. This paper describes implementation details, benchmark results, and the advantages and technical challenges with this approach. [ABSTRACT FROM AUTHOR]

The rotationally resolved infrared (IR) spectrum of the He– H 3 + complex has been measured in a cryogenic ion trap experiment at a nominal temperature of 4 K. Predissociation of the stored complex has been invoked by excitation of the degenerate ν2 mode of the H 3 + sub-unit using a pulsed optical parametric oscillator system. An assignment of the experimental spectrum became possible through one-to-one correlations with bands of the spectrum theoretically predicted in Paper I [Harding et al., J. Chem. Phys. 156, 144307 (2022)]. 19 bands have been assigned and analyzed, and the energy term diagram of the lower states of this floppy molecular complex has been derived from combination differences (CDs) in the experimental spectrum. Ground state combination differences (GSCDs) reveal a large part of the energy term diagram for the He– H 3 + complex in its vibrational ground state, v = 0. Experimental and theoretical term energies agree within experimental accuracy for the rotational fine structure associated with the total angular momentum quantum number J and the parity e/f as well as for the coarse spacing of the lowest K states of the complex. This favorable comparison shows that the potential energy surface (PES) calculated in Paper I is accurate. The barriers between the three equivalent global minima in this PES are relatively low and the He– H 3 + complex is extremely floppy, with nearly unhindered internal rotation of the H 3 + sub-unit. The resulting Coriolis interactions couple the internal and end-over-end rotation of the complex and contribute significantly to the energy terms. They are observed both in experiment and theory and are, e.g., the origin of different rotational constants for states of e and f parity. Also in this respect, experiment and theory agree very well. Despite the assignment and analysis of many bands of the extremely rich IR spectrum of He– H 3 + , higher levels of excitation, including the complex stretching mode, need further attention. [ABSTRACT FROM AUTHOR]

Vibrational spectroscopy is a useful technique for probing chemical environments. The development of models that can reproduce the spectra of nitriles and azides is valuable because these probes are uniquely suited for investigating complex systems. Empirical vibrational spectroscopic maps are commonly employed to obtain the instantaneous vibrational frequencies during molecular dynamics simulations but often fail to adequately describe the behavior of these probes, especially in its transferability to a diverse range of environments. In this paper, we demonstrate several reasons for the difficulty in constructing a general-purpose vibrational map for methyl thiocyanate (MeSCN), a model for cyanylated biological probes. In particular, we found that electrostatics alone are not a sufficient metric to categorize the environments of different solvents, and the dominant features in intermolecular interactions in the energy landscape vary from solvent to solvent. Consequently, common vibrational mapping schemes do not cover all essential interaction terms adequately, especially in the treatment of van der Waals interactions. Quantum vibrational perturbation (QVP) theory, along with a combined quantum mechanical and molecular mechanical potential for solute–solvent interactions, is an alternative and efficient modeling technique, which is compared in this paper, to yield spectroscopic results in good agreement with experimental FTIR. QVP has been used to analyze the computational data, revealing the shortcomings of the vibrational maps for MeSCN in different solvents. The results indicate that insights from QVP analysis can be used to enhance the transferability of vibrational maps in future studies. [ABSTRACT FROM AUTHOR]

Chemical shakers are active particles with zero propulsion velocity whose activity derives from chemical reactions on portions of their surfaces. Although they do not move, except through Brownian motion, the nonequilibrium concentration and velocity fields that they generate endow them with properties that differ from their equilibrium counterparts. In particular, collections of such shakers can actively move, reorient, and self-assemble into various cluster states, which are the subject of this paper. Elongated chemical shakers constructed from linked catalytic and noncatalytic spheres are considered, and it is shown how hydrodynamic, chemotactic, and shape-dependent interactions give rise to various self-assembled shaker structures. The chemical forces responsible for cluster formation are described in terms of a model based on pair-wise additive contributions. The forms of the self-assembled structures can be varied by changing the chemostat concentrations that control the nonequilibrium state. The resulting structures and their manipulation through chemical means suggest ways to construct a class of active materials for applications. [ABSTRACT FROM AUTHOR]

Packaging insulation materials with high thermal conductivity and excellent dielectric properties are favorable to meet the high demand and rapid development of third generation power semiconductors. In this study, we propose to improve the thermal conductivity of epoxy resin (EP) by incorporating a three-dimensional boron nitride thermally conductive network. Detailedly, polyurethane foam (PU) was used as a supporter, and boron nitride nanosheets (BNNSs) were loaded onto the PU supporter through chemical bonding (BNNS@PU). After immersing BNNS@PU into the EP resin, EP-based thermally conductive composites were prepared by vacuum-assisted impregnation. Fourier transform infrared spectrometer and scanning electron microscope were used to characterize the chemical bonding and morphological structure of BNNS@PU, respectively. The content of BNNS in BNNS@PU/EP composites was quantitatively analyzed by TGA. The results show that the thermal conductivity of the BNNS@PU/EP composites reaches 0.521 W/m K with an enhancement rate η of 30.89 at an ultra-low BNNS filler content (5.93 wt. %). Additionally, the BNNS@PU/EP composites have excellent dielectric properties with the frequency range from 101 to 106 Hz. This paper provides an interesting idea for developing high thermal conductivity insulating materials used for power semiconductor packaging. [ABSTRACT FROM AUTHOR]

This paper explores the unintentional contamination of Surface-Enhanced Raman Scattering (SERS) substrates by ambient hydrocarbon contaminants and their contribution to SERS spectra. Previous studies have identified amorphous carbon as a potential complicating factor in data analysis in SERS experiments, although its origin has been elusive. Our work showed that ambient hydrocarbon contamination and its decomposition products can be detected by SERS on a gold substrate. We propose that ambient air itself is a source of amorphous carbon contamination on SERS substrates. This understanding is crucial for the correct interpretation of SERS data and highlights the need for careful consideration of potential environmental contaminants in SERS analysis. [ABSTRACT FROM AUTHOR]

Semiconductor nanocrystals (NCs) can function as efficient gain materials with chemical versatility because of their surface ligands. Because the properties of NCs in solution are sensitive to ligand–environment interactions, local chemical changes can result in changes in the optical response. However, amplification of the optical response is technically challenging because of colloidal instability at NC concentrations needed for sufficient gain to overcome losses. This paper demonstrates liquid lasing from plasmonic lattice cavities integrated with ligand-engineered CdZnS/ZnS NCs dispersed in toluene and water. By taking advantage of calcium ion-induced aggregation of NCs in aqueous solutions, we show how lasing threshold can be used as a transduction signal for ion detection. Our work highlights how NC solutions and plasmonic lattices with open cavity architectures can serve as a biosensing platform for lab-on-chip devices. [ABSTRACT FROM AUTHOR]

In this paper, the responses of assembled behaviors of block polyelectrolytes (PEs) to the strength of electrostatic interactions are studied through molecular dynamic simulations. The results show that the assembled structures closely depend on the electrostatic strength. It should be noted that PE coacervation can outweigh the nucleation of hydrophobic blocks and invert the micelle structures at strong electrostatic strengths, leading to the formation of inverted micelles of PE cores and hydrophobic coronas. In the poor solvent condition for neutral block, diverse anisotropic micelles are presented; candy-like conventional micelles of hydrophobic cores and PE patches coexist with inverted candy-like micelles of PE cores and hydrophobic patches and with Janus micelles of semi-neutral aggregate and semi-PE cluster in the presence of divalent and trivalent counterions. The formation of conventional or inverted micelle is largely determined by the type of micellar fusion, which results from the nucleation competition between electrostatic correlation and hydrophobic interaction. The merge of micelles mediated by hydrophobic attraction leads to conventional hydrophobic cores, and the fusion induced by electrostatic correlations results in PE cores micelles. At strong electrostatic strengths, the PE chains exhibit rich conformations at trivalent counterions, ranging from a fully collapsed state to a rod-like state, and parallel alignment of PE chains is found. [ABSTRACT FROM AUTHOR]

Ion channels are proteins that play a significant role in physiological processes, including neuronal excitability and signal transduction. However, the precise mechanisms by which these proteins facilitate ion diffusion through cell membranes are not well understood. This is because experimental techniques to characterize ion channel activity operate on a time scale too large to understand the role of the various protein conformations on diffusion. Meanwhile, computational approaches operate on a time scale too short to rationalize the observed behavior at the microscopic scale. In this paper, we present a continuous-time random walk model that aims to bridge the scales between the atomistic models of ion channels and the experimental measurement of their conductance. We show how diffusion slows down in complex systems by using 3D lattices that map out the pore geometry of two channels: Nav1.7 and gramicidin. We also introduce spatial and dynamic site disorder to account for system heterogeneity beyond the mean field approximation. Computed diffusion coefficients show that an increase in spatial disorder slows down diffusion kinetics, while dynamic disorder has the opposite effect. Our results imply that microscopic or phenomenological models based on the potential of mean force data overlook the functional importance of protein dynamics on ion diffusion through channels. [ABSTRACT FROM AUTHOR]

Conical intersection (CI) leads to fast electronic energy transfer. However, Hamm and Stock [Phys. Rev. Lett. 109, 173201 (2012)] showed the existence of a vibrational CI and its role in vibrational energy relaxation. In this paper, we further investigate the vibrational energy relaxation using an isolated model Hamiltonian system of four vibrational modes with two distinctively different timescales (two fast modes and two slow modes). We show that the excitation of the slow modes plays a crucial role in the energy relaxation mechanism. We also analyze the system from a mixed quantum-classical (surface hopping method) and a completely classical point of view. Notably, surface hopping and even classical simulations also capture fast energy relaxation, which is a signature of CI's existence. [ABSTRACT FROM AUTHOR]

Alchemical transformations showed that perturbation theory can be applied also to changes in the atomic nuclear charges of a molecule. The alchemical path that connects two different chemical species involves the conceptualization of a non-physical system in which an atom possess a non-integer nuclear charge. A correct quantum mechanical treatment of these systems is limited by the fact that finite size atomic basis sets do not define exponents and contraction coefficients for fractional charge atoms. This paper proposes a solution to this problem and shows that a smooth interpolation of the atomic orbital coefficients and exponents across the periodic table is a convenient way to produce accurate alchemical predictions, even using small size basis sets. [ABSTRACT FROM AUTHOR]

This paper introduces a spin-free formulation of the supporting subspace factorization [C. Song and T. J. Martínez, J. Chem. Phys. 149, 044108 (2018)], enabling a reduction in the computational scaling of the extended multi-state complete active space second-order perturbation (XMS-CASPT2) method for arbitrary spins. Compared to the original formulation that is defined in the spin orbitals and is limited to singlet states, the spin-free formulation in this work treats different spin states equivalently, thus naturally generalizing the idea beyond singlet states. In addition, we will present a new way of deriving the supporting subspace factorization with the purpose of understanding its physical interpretation. In this new derivation, we separate the sources that make CASPT2 difficult into the "same-site interactions" and "inter-site interactions." We will first show how the Kronecker sum can be used to remove the same-site interactions in the absence of inter-site interactions, leading to MP2 energy in dressed orbitals. We will then show how the inter-site interactions can be exactly recovered using Löwdin partition, where the supporting subspace concept will naturally arise. The new spin-free formulation maintains the main advantage of the supporting subspace factorization, i.e., allowing XMS-CASPT2 energies to be computed using highly optimized MP2 energy codes and Fock build codes, thus reducing the scaling of XMS-CASPT2 to the same scaling as MP2. We will present and discuss results that benchmark the accuracy and performance of the new method. To demonstrate how the new method can be useful in studying real photochemical systems, the supporting subspace XMS-CASPT2 is applied to a photoreaction sensitive to magnetic field effects. The new spin-free formulation makes it possible to calculate the doublet and quartet states required in this particular photoreaction mechanism. [ABSTRACT FROM AUTHOR]

Within the context of fewest-switch surface hopping (FSSH) dynamics, one often wishes to remove the angular component of the derivative coupling between states J and K . In a previous set of papers, Shu et al. [J. Phys. Chem. Lett. 11, 1135–1140 (2020)] posited one approach for such a removal based on direct projection, while we isolated a second approach by constructing and differentiating a rotationally invariant basis. Unfortunately, neither approach was able to demonstrate a one-electron operator O ̂ whose matrix element J O ̂ K was the angular component of the derivative coupling. Here, we show that a one-electron operator can, in fact, be constructed efficiently in a semi-local fashion. The present results yield physical insight into designing new surface hopping algorithms and are of immediate use for FSSH calculations. [ABSTRACT FROM AUTHOR]

To address the challenge of performance portability and facilitate the implementation of electronic structure solvers, we developed the basic matrix library (BML) and Parallel, Rapid O(N), and Graph-based Recursive Electronic Structure Solver (PROGRESS) library. The BML implements linear algebra operations necessary for electronic structure kernels using a unified user interface for various matrix formats (dense and sparse) and architectures (CPUs and GPUs). Focusing on density functional theory and tight-binding models, PROGRESS implements several solvers for computing the single-particle density matrix and relies on BML. In this paper, we describe the general strategies used for these implementations on various computer architectures, using OpenMP target functionalities on GPUs, in conjunction with third-party libraries to handle performance critical numerical kernels. We demonstrate the portability of this approach and its performance in benchmark problems. [ABSTRACT FROM AUTHOR]

The diffusion of dilute molecular penetrants within polymers plays a crucial role in the advancement of material engineering for applications such as coatings and membrane separations. The potential of highly cross-linked polymer networks in these applications stems from their capacity to adjust the size and shape selectivity through subtle changes in network structures. In this paper, we use molecular dynamics simulation to understand the role of penetrant shape (aspect ratios) and its interaction with polymer networks on its diffusivity. We characterize both local penetrant hopping and the long-time diffusive motion for penetrants and consider different aspect ratios and penetrant–network interaction strengths at a variety of cross-link densities and temperatures. The shape affects the coupling of penetrant motion to the cross-link density- and temperature-dependent structural relaxation of networks and also affects the way a penetrant experiences the confinement from the network meshes. The attractive interaction between the penetrant and network primarily affects the former since only the system of dilute limit is of present interest. These results offer fundamental insights into the intricate interplay between penetrant characteristics and polymer network properties and also suggest future directions for manipulating polymer design to enhance the separation efficiency. [ABSTRACT FROM AUTHOR]

In this paper, we investigate the performance of different machine learning potentials (MLPs) in predicting key thermodynamic properties of water using RPBE + D3. Specifically, we scrutinize kernel-based regression and high-dimensional neural networks trained on a highly accurate dataset consisting of about 1500 structures, as well as a smaller dataset, about half the size, obtained using only on-the-fly learning. This study reveals that despite minor differences between the MLPs, their agreement on observables such as the diffusion constant and pair-correlation functions is excellent, especially for the large training dataset. Variations in the predicted density isobars, albeit somewhat larger, are also acceptable, particularly given the errors inherent to approximate density functional theory. Overall, this study emphasizes the relevance of the database over the fitting method. Finally, this study underscores the limitations of root mean square errors and the need for comprehensive testing, advocating the use of multiple MLPs for enhanced certainty, particularly when simulating complex thermodynamic properties that may not be fully captured by simpler tests. [ABSTRACT FROM AUTHOR]

This article describes a method for calculating higher order or nonadiabatic corrections in Born–Oppenheimer theory and its interaction with the translational degrees of freedom. The method uses the Wigner–Weyl correspondence to map nuclear operators into functions on the classical phase space and the Moyal star product to represent operator multiplication on those functions. These are explained in the body of the paper. The result is a power series in κ2, where κ = (m/M)1/4 is the usual Born–Oppenheimer parameter. The lowest order term is the usual Born–Oppenheimer approximation, while higher order terms are nonadiabatic corrections. These are needed in calculations of electronic currents, momenta, and densities. The separation of nuclear and electronic degrees of freedom takes place in the context of the exact symmetries (for an isolated molecule) of translations and rotations, and these, especially translations, are explicitly incorporated into our discussion. This article presents an independent derivation of the Moyal expansion in molecular Born–Oppenheimer theory. We show how electronic currents and momenta can be calculated within the framework of Moyal perturbation theory; we derive the transformation laws of the electronic Hamiltonian, the electronic eigenstates, and the derivative couplings under translations; we discuss in detail the rectilinear motion of the molecular center of mass in the Born–Oppenheimer representation; and we show how the elimination of the translational components of the derivative couplings leads to a unitary transformation that has the effect of exactly separating the translational degrees of freedom. [ABSTRACT FROM AUTHOR]

We introduce a generalized micro–macro Markov chain Monte Carlo (mM-MCMC) method with pseudo-marginal approximation to the free energy that is able to accelerate sampling of the microscopic Gibbs distributions when there is a time-scale separation between the macroscopic dynamics of a reaction coordinate and the remaining microscopic degrees of freedom. The mM-MCMC method attains this efficiency by iterating four steps: (i) propose a new value of the reaction coordinate, (ii) accept or reject the macroscopic sample, (iii) run a biased simulation that creates a microscopic molecular instance that lies close to the newly sampled macroscopic reaction coordinate value, and (iv) microscopic accept/reject step for the new microscopic sample. In the present paper, we eliminate the main computational bottleneck of earlier versions of this method: the necessity to have an accurate approximation of free energy. We show that the introduction of a pseudo-marginal approximation significantly reduces the computational cost of the microscopic accept/reject step while still providing unbiased samples. We illustrate the method's behavior on several molecular systems with low-dimensional reaction coordinates. [ABSTRACT FROM AUTHOR]

Protein–calixarenes binding plays an increasingly central role in many applications, spanning from molecular recognition to drug delivery strategies and protein inhibition. These ligands obey a specific bio-supramolecular chemistry, which can be revealed by computational approaches, such as molecular dynamics simulations. In this paper, we rely on all-atom, explicit-solvent molecular dynamics simulations to capture the electrostatically driven association of a phosphonated calix-[4]-arene with cytochome-C, which critically relies on surface-exposed paired lysines. Beyond two binding sites identified in direct agreement with the x-ray structure, the association has a larger structural impact on the protein dynamics. Then, our simulations allow a direct comparison to analogous calixarenes, namely, sulfonato, similarly reported as "molecular glue." Our work can contribute to a robust in silico predictive tool to assess binding sites for any given protein of interest for crystallization, with the specificity of a macromolecular cage whose endo/exo orientation plays a role in the binding. [ABSTRACT FROM AUTHOR]

In this paper, we develop and analyze a number of perturbation series that target the coupled cluster singles and doubles (CCSD) ground state energy. We show how classical Møller–Plesset perturbation theory series can be restructured to target the CCSD energy based on a reference CCS calculation and how the corresponding cluster perturbation series differs from the classical Møller–Plesset perturbation series. Subsequently, we reformulate these series using the coupled cluster Lagrangian framework to obtain series, where fourth and fifth order energies are determined only using parameters through second order. To test the methods, we perform a series of test calculations on molecular photoswitches of both total energies and reaction energies. We find that the fifth order reaction energies are of CCSD quality and that they are of comparable accuracy to state-of-the-art approximations to the CCSD energy based on local pair natural orbitals. The advantage of the present approach over local correlation methods is the absence of user defined threshold parameters for neglecting or approximating contributions to the correlation energy. Fixed threshold parameters lead to discontinuous energy surfaces, although this effect is often small enough to be ignored, but the present approach has a differentiable energy that will facilitate derivation and implementation of gradients and higher derivatives. A further advantage is that the calculation of the perturbation correction is non-iterative and can, therefore, be calculated in parallel, leading to a short time-to-solution. [ABSTRACT FROM AUTHOR]

Time-Resolved Serial Femtosecond Crystallography (TR-SFX) conducted at X-ray Free Electron Lasers (XFELs) has become a powerful tool for capturing macromolecular structural movies of light-initiated processes. As the capabilities of XFELs advance, we anticipate that a new range of coherent control and structural Raman measurements will become achievable. Shorter optical and x-ray pulse durations and increasingly more exotic pulse regimes are becoming available at free electron lasers. Moreover, with high repetition enabled by the superconducting technology of European XFEL (EuXFEL) and Linac Coherent Light Source (LCLS-II) , it will be possible to improve the signal-to-noise ratio of the light-induced differences, allowing for the observation of vibronic motion on the sub-Angstrom level. To predict and assign this coherent motion, which is measurable with a structural technique, new theoretical approaches must be developed. In this paper, we present a theoretical density matrix approach to model the various population and coherent dynamics of a system, which considers molecular system parameters and excitation conditions. We emphasize the use of the Wigner transform of the time-dependent density matrix, which provides a phase space representation that can be directly compared to the experimental positional displacements measured in a TR-SFX experiment. Here, we extend the results from simple models to include more realistic schemes that include large relaxation terms. We explore a variety of pulse schemes using multiple model systems using realistic parameters. An open-source software package is provided to perform the density matrix simulation and Wigner transformations. The open-source software allows us to define any arbitrary level schemes as well as any arbitrary electric field in the interaction Hamiltonian. [ABSTRACT FROM AUTHOR]

In order to prepare highly heat-resistant packaging insulation materials, in this paper, bismaleimide/epoxy resin (BMI/EP55) composites with different contents of BMI were prepared by melt blending BMI into amino tetrafunctional and phenolic epoxy resin (at a ratio of 5:5). The microstructures and thermal and electrical properties of the composites were tested. The electrostatic potential distribution, energy level distribution, and molecular orbitals of BMI were calculated using Gaussian. The results showed that the carbonyl group in BMI is highly electronegative, implying that the carbonyl group has a strong electron trapping ability. The thermal decomposition temperature of the composites gradually increased with the increase of BMI content, and the 20% BMI/EP55 composites had the highest heat-resistance index, along with a glass transition temperature (Tg) of >250 °C. At different test temperatures, with increase in the BMI content, the conductivity of epoxy resin composites showed a tendency to first decrease and then increase, the breakdown field strength showed a tendency to first increase and then decrease, and the dielectric constant was gradually decreased. Two trap centers were present simultaneously in the composites, where the shallow trap energy level is the deepest in 20% BMI/EP composites and the deep trap energy level is the deepest in 10% BMI/EP55 composites. Correspondingly, the 10% BMI/EP55 composite had a slower charge decay rate, while the 20% BMI/EP55 had a faster charge decay rate. In summary, the BMI/EP55 composites with high heat resistance and insulating properties were prepared in this study, which provided ideas for preparing high-temperature packaging insulating materials. [ABSTRACT FROM AUTHOR]

The large discrepancy among the nucleation kinetics extracted from experimental measurements and computer simulations and the prediction of the classical nucleation theory (CNT) has stimulated intense arguments about its origin in the past decades, which is crucially relevant to the validity of the CNT. In this paper, we investigate the atomistic mechanism of the nucleation in liquid Al in contact with amorphous substrates with atomic-level smooth/rough surfaces, using molecular dynamics (MD) simulations. This study reveals that the slightly distorted local fcc/hcp structures in amorphous substrates with smooth surfaces can promote heterogeneous nucleation through a structural templating mechanism, and on the other hand, homogeneous nucleation will occur at a larger undercooling through a fluctuation mechanism if the surface is rough. Thus, some impurities, previously thought to be impotent, could be activated in the homogeneous nucleation experiments. We further find that the initial growth of the nucleus on smooth surfaces of amorphous substrates is one order of magnitude faster than that in homogeneous nucleation. Both these factors could significantly contribute to the discrepancy in the nucleation kinetics. This study is also supported by a recent study of the synthesis of high-entropy alloy nanoparticles assisted with the liquid metal Ga [Cao et al., Nature 619, 73 (2023)]. In this study, we established that the boundary existed between homogeneous and heterogeneous nucleation, i.e., the structural templating is a general mechanism for heterogeneous nucleation, and in its absence, homogeneous nucleation will occur through the fluctuation mechanism. This study provides an in-depth understanding of the nucleation theory and experiments. [ABSTRACT FROM AUTHOR]

Methyl nitrite has two stable conformational isomers resulting from rotation about the primary C–O–N–O dihedral angle: cis-CH3ONO and trans-CH3ONO, with cis being more stable by ∼5 kJ/mol. The barrier to rotational interconversion (∼45 kJ/mol) is too large for isomerization to occur under ambient conditions. This paper presents evidence of a change in conformer abundance when dilute CH3ONO is deposited onto a cold substrate; the relative population of the freshly deposited cis conformer is seen to increase compared to its gas-phase abundance, measured by in situ infrared spectroscopy. We observe abundance changes depending on the identity of the bath gas (N2, Ar, and Xe) and deposition angle. The observations indicate that the surface properties of the growing matrix influence conformer abundance—contrary to the widely held assumption that conformer abundance in matrices reflects gas-phase abundance. We posit that differences in the angle-dependent host-gas deposition dynamics affect the growing surfaces, causing changes in conformer abundances. Quantum chemistry calculations of the binding energies between CH3ONO and a single bath-gas component reveal that significant energetic stabilization is not observed in 1:1 complexes of N2:CH3ONO, Ar:CH3ONO, or Xe:CH3ONO. From our results, we conclude that the growing surface plays a significant role in trapping cis-CH3ONO more effectively than trans-CH3ONO, likely because cis-CH3ONO is more compact. Taken together, the observations highlight the necessity for careful characterization of conformers in matrix-isolated systems, emphasizing a need for further study into the deposition dynamics and surface structure of chemically inert matrices. [ABSTRACT FROM AUTHOR]

In this paper, we report on the implementation of CC2 and CC3 in the context of molecules in finite magnetic fields. The methods are applied to the investigation of atoms and molecules through spectroscopic predictions and geometry optimizations for the study of the atmosphere of highly magnetized White Dwarf stars. We show that ground-state finite-field (ff) CC2 is a reasonable alternative to CCSD for energies and, in particular, for geometrical properties. For excited states, ff-CC2 is shown to perform well for states with predominant single-excitation character. Yet, for cases in which the excited state wavefunction has double-excitation character with respect to the reference, ff-CC2 can easily lead to completely unphysical results. Ff-CC3, however, is shown to reproduce the CCSDT behavior very well and enables the treatment of larger systems at a high accuracy. [ABSTRACT FROM AUTHOR]

Conical intersections are known to play a vital role in many photochemical processes. The breakdown of the Born–Oppenheimer approximation in the vicinity of a conical intersection causes exciting phenomena, such as the ultrafast radiationless decay of excited states. The passage of a molecule through a conical intersection creates a coherent superposition of electronic states via nonadiabatic couplings. Detecting this coherent superposition may serve as a direct probe of the conical intersection. In this paper, we theoretically demonstrate the use of coherent spontaneous emission in samples with long-range order for probing the occurrence of a conical intersection in a molecule. Our simulations show that the spectrum contains clear signatures of the created coherent superposition of electronic states. We investigate the bandwidth requirements for the x-ray probes, which influence the observation of coherent superposition generated by the conical intersection. [ABSTRACT FROM AUTHOR]

Semiflexible slender filaments are ubiquitous in nature and cell biology, including in the cytoskeleton, where reorganization of actin filaments allows the cell to move and divide. Most methods for simulating semiflexible inextensible fibers/polymers are based on discrete (bead-link or blob-link) models, which become prohibitively expensive in the slender limit when hydrodynamics is accounted for. In this paper, we develop a novel coarse-grained approach for simulating fluctuating slender filaments with hydrodynamic interactions. Our approach is tailored to relatively stiff fibers whose persistence length is comparable to or larger than their length and is based on three major contributions. First, we discretize the filament centerline using a coarse non-uniform Chebyshev grid, on which we formulate a discrete constrained Gibbs–Boltzmann (GB) equilibrium distribution and overdamped Langevin equation for the evolution of unit-length tangent vectors. Second, we define the hydrodynamic mobility at each point on the filament as an integral of the Rotne–Prager–Yamakawa kernel along the centerline and apply a spectrally accurate "slender-body" quadrature to accurately resolve the hydrodynamics. Third, we propose a novel midpoint temporal integrator, which can correctly capture the Ito drift terms that arise in the overdamped Langevin equation. For two separate examples, we verify that the equilibrium distribution for the Chebyshev grid is a good approximation of the blob-link one and that our temporal integrator for overdamped Langevin dynamics samples the equilibrium GB distribution for sufficiently small time step sizes. We also study the dynamics of relaxation of an initially straight filament and find that as few as 12 Chebyshev nodes provide a good approximation to the dynamics while allowing a time step size two orders of magnitude larger than a resolved blob-link simulation. We conclude by applying our approach to a suspension of cross-linked semiflexible fibers (neglecting hydrodynamic interactions between fibers), where we study how semiflexible fluctuations affect bundling dynamics. We find that semiflexible filaments bundle faster than rigid filaments even when the persistence length is large, but show that semiflexible bending fluctuations only further accelerate agglomeration when the persistence length and fiber length are of the same order. [ABSTRACT FROM AUTHOR]

Second-order N-electron valence state perturbation theory (NEVPT2) is an exactly size-consistent and intruder-state-free multi-reference theory. To accelerate the NEVPT2 computation, Guo and Neese combined it with the local pair-natural orbital (PNO) method using the projected atomic orbitals (PAOs) as the underlying local basis [Guo et al., J. Chem. Phys. 144, 094111 (2016)]. In this paper, we report the further development of the PNO-NEVPT2 method using the orthonormal and non-redundant localized virtual molecular orbitals (LVMOs) instead of PAOs. The LVMOs were previously considered to perform relatively poor compared to PAOs because the resulting orbital domains were unacceptably large. Our prior work, however, showed that this drawback can be remedied by re-forming the domain construction scheme using differential overlap integrals [Saitow et al., J. Chem. Phys. 157, 084101 (2022)]. In this work, we develop further refinements to enhance the feasibility of using LVMOs. We first developed a two-level semi-local approach for screening out so-called weak-pairs to select or truncate the pairs for PNO constructions more flexibly. As a refinement specific to the Pipek–Mezey localization for LVMOs, we introduced an iterative scheme to truncate the Givens rotations using varying thresholds. We assessed the LVMO-based PNO-NEVPT2 method through benchmark calculations for linear phenyl alkanes, which demonstrate that it performs comparably well relative to the PAO-based approach. In addition, we evaluated the Co–C bond dissociation energies for the cobalamin derivatives composed of 200 or more atoms, which confirms that the LVMO-based method can recover more than 99.85% of the canonical NEVPT2 correlation energy. [ABSTRACT FROM AUTHOR]

In two recent papers [A. D. Becke, J. Chem. Phys. 156, 214101 (2022) and A. D. Becke, J. Chem. Phys. 157, 234102 (2022)], we compared two Kohn–Sham density functionals based on physical modeling and theory with the best density-functional power-series fits in the literature. The best error statistics reported to date for a hybrid functional on the general main-group thermochemistry, kinetics, and noncovalent interactions (GMTKN55) chemical database of Goerigk et al. [Phys. Chem. Chem. Phys. 19, 32184 (2017)] were obtained. In the present work, additional second-order perturbation-theory terms are considered. The result is a 12-parameter double-hybrid density functional with the lowest GMTKN55 WTMAD2 "weighted total mean absolute deviation" error (1.76 kcal/mol) yet seen for any hybrid or double-hybrid density-functional approximation. We call it "DH23." [ABSTRACT FROM AUTHOR]

This paper aims at the theoretical study of the CO2 magnetic-dipole ν2 + ν3 rovibrational absorption band that was recently detected in the Martian atmosphere. Specific characteristics of the magnetic dipole operator are carefully examined. Our evaluation of the magnetic-dipole line intensities is based on the variational calculations and the use of molecular properties is determined through specially performed ab initio quantum chemical calculations. The comparison of our simulated magnetic-dipole spectrum with available laboratory taken data also requires the knowledge of line intensities in the quadrupole band, which partially overlaps with that magnetic-dipole. Quadrupole intensities, once reconsidered, are permitted to correct previously reported values of the integrated intensity as well as the intensity of selected branches. The sum of our calculated magnetic-dipole and quadrupole rovibrational lines is shown to be in good agreement with both sets of presently available data from FTIR and OFCEAS laboratory observations. [ABSTRACT FROM AUTHOR]

Metal oxide nanoparticle (NP) supports of both good conductivity and stability have the potential to enhance both the reaction activity and stability of the loaded electrocatalysts. In this paper, a facile two-step approach to disperse Pt nanoparticles on the surface of an IrO2 NP support (Pt/IrO2) was developed. Physical characterization by x-ray diffraction spectroscopy and transmission/scanning electron microscopy suggests a good dispersion of the Pt NPs. The temperature effect (from 293 to 353 K) of oxygen reduction reaction on Pt/IrO2 was studied by using a rotating ring disk electrode The results show that although the kinetic current density on Pt/IrO2 is close to that on commercial Pt/C at room temperature, the apparent activation energy (Ea,app) in the former case is much lower, suggesting a much higher activity at elevated temperatures. The superiority in Ea,app is attributed to the electron interaction between Pt and the IrO2 support, as supported by the change of surface chemical state given by x-ray photo-electron spectroscopy. [ABSTRACT FROM AUTHOR]

We report further investigations to aid the development of a Laplace MP2 (second-order Møller Plesset) method with a range separated Coulomb potential partitioned into short- and long-range parts. The implementation of the method extensively uses sparse matrix algebra, density fitting techniques for the short-range part, and a Fourier transformation in spherical coordinates for the long-range part of the potential. Localized molecular orbitals are employed for the occupied space, whereas virtual space is described by orbital specific virtual orbitals (OSVs) associated with localized molecular orbitals. The Fourier transform is deficient for very large distances between localized occupied orbitals, and a multipole expansion for widely separated pairs is introduced for the direct MP2 contribution, which is applicable also to non-Coulombic potentials that do not satisfy the Laplace equation. For the exchange contribution, an efficient screening of contributing localized occupied pairs is employed, which is discussed more completely here. To mitigate errors due to the truncation of OSVs, a simple and efficient extrapolation procedure is used to obtain results close to MP2 for the full basis set of atomic orbitals Using a suitable set of default parameters, the accuracy of the approach is demonstrated. The current implementation of the approach is not very efficient, and the aim of this paper is to introduce and critically discuss ideas that can have more general applicability beyond MP2 calculations for large molecules. [ABSTRACT FROM AUTHOR]

The complexity of mathematical models in biology has rendered model reduction an essential tool in the quantitative biologist's toolkit. For stochastic reaction networks described using the Chemical Master Equation, commonly used methods include time-scale separation, Linear Mapping Approximation, and state-space lumping. Despite the success of these techniques, they appear to be rather disparate, and at present, no general-purpose approach to model reduction for stochastic reaction networks is known. In this paper, we show that most common model reduction approaches for the Chemical Master Equation can be seen as minimizing a well-known information-theoretic quantity between the full model and its reduction, the Kullback–Leibler divergence defined on the space of trajectories. This allows us to recast the task of model reduction as a variational problem that can be tackled using standard numerical optimization approaches. In addition, we derive general expressions for propensities of a reduced system that generalize those found using classical methods. We show that the Kullback–Leibler divergence is a useful metric to assess model discrepancy and to compare different model reduction techniques using three examples from the literature: an autoregulatory feedback loop, the Michaelis–Menten enzyme system, and a genetic oscillator. [ABSTRACT FROM AUTHOR]

This paper concerns the representation of angular momentum operators in the Born–Oppenheimer theory of polyatomic molecules and the various forms of the associated conservation laws. Topics addressed include the question of whether these conservation laws are exactly equivalent or only to some order of the Born–Oppenheimer parameter κ = (m/M)1/4 and what the correlation is between angular momentum quantum numbers in the various representations. These questions are addressed in both problems involving a single potential energy surface and those with multiple, strongly coupled surfaces and in both the electrostatic model and those for which fine structure and electron spin are important. The analysis leads to an examination of the transformation laws under rotations of the electronic Hamiltonian; of the basis states, both adiabatic and diabatic, along with their phase conventions; of the potential energy matrix; and of the derivative couplings. These transformation laws are placed in the geometrical context of the structures in the nuclear configuration space that are induced by rotations, which include the rotational orbits or fibers, the surfaces upon which the orientation of the molecule changes but not its shape, and the section, an initial value surface that cuts transversally through the fibers. Finally, it is suggested that the usual Born–Oppenheimer approximation can be replaced by a dressing transformation, that is, a sequence of unitary transformations that block-diagonalize the Hamiltonian. When the dressing transformation is carried out, we find that the angular momentum operator does not change. This is a part of a system of exact equivalences among various representations of angular momentum operators in Born–Oppenheimer theory. Our analysis accommodates large-amplitude motions and is not dependent on small-amplitude expansions about an equilibrium position. Our analysis applies to noncollinear configurations of a polyatomic molecule; this covers all but a subset of measure zero (the collinear configurations) in the nuclear configuration space. [ABSTRACT FROM AUTHOR]

The interactions between molecules and electrode surfaces play a key role in electrochemical processes and are a subject of extensive research, both experimental and theoretical. In this paper, we address the water dissociation reaction on a Pd(111) electrode surface, modeled as a slab embedded in an external electric field. We aim at unraveling the relationship between surface charge and zero-point energy in aiding or hindering this reaction. We calculate the energy barriers with dispersion-corrected density-functional theory and an efficient parallel implementation of the nudged-elastic-band method. We show that the lowest dissociation barrier and consequently the highest reaction rate take place when the field reaches a strength where two different geometries of the water molecule in the reactant state are equally stable. The zero-point energy contributions to this reaction, on the other hand, remain nearly constant across a wide range of electric field strengths, despite significant changes in the reactant state. Interestingly, we show that the application of electric fields that induce a negative charge on the surface can make nuclear tunneling more significant for these reactions. [ABSTRACT FROM AUTHOR]

Anti-freeze proteins (AFPs) protect organisms at freezing conditions by attaching to the ice surface and arresting its growth. Each adsorbed AFP locally pins the ice surface, resulting in a metastable dimple for which the interfacial forces counteract the driving force for growth. As supercooling increases, these metastable dimples become deeper, until metastability is lost in an engulfment event where the ice irreversibly swallows the AFP. Engulfment resembles nucleation in some respects, and this paper develops a model for the "critical profile" and free energy barrier for the engulfment process. Specifically, we variationally optimize the ice–water interface and estimate the free energy barrier as a function of the supercooling, the AFP footprint size, and the distance to neighboring AFPs on the ice surface. Finally, we use symbolic regression to derive a simple closed-form expression for the free energy barrier as a function of two physically interpretable, dimensionless parameters. [ABSTRACT FROM AUTHOR]

Trapping diffusive particles at surfaces is a key step in many systems in chemical and biological physics. Trapping often occurs via reactive patches on the surface and/or the particle. The theory of boundary homogenization has been used in many prior works to estimate the effective trapping rate for such a system in the case that either (i) the surface is patchy and the particle is uniformly reactive or (ii) the particle is patchy and the surface is uniformly reactive. In this paper, we estimate the trapping rate for the case that the surface and the particle are both patchy. In particular, the particle diffuses translationally and rotationally and reacts with the surface when a patch on the particle contacts a patch on the surface. We first formulate a stochastic model and derive a five-dimensional partial differential equation describing the reaction time. We then use matched asymptotic analysis to derive the effective trapping rate, assuming that the patches are roughly evenly distributed and occupy a small fraction of the surface and the particle. This trapping rate involves the electrostatic capacitance of a four-dimensional duocylinder, which we compute using a kinetic Monte Carlo algorithm. We further use Brownian local time theory to derive a simple heuristic estimate of the trapping rate and show that it is remarkably close to the asymptotic estimate. Finally, we develop a kinetic Monte Carlo algorithm to simulate the full stochastic system and then use these simulations to confirm the accuracy of our trapping rate estimates and homogenization theory. [ABSTRACT FROM AUTHOR]

In this paper, we perform the exact diagonalization of a light–matter strongly coupled system taking into account arbitrary losses via both energy dissipation in the optically active material and photon escape out of the resonator. This allows us to naturally treat the cases of couplings with structured reservoirs, which can strongly impact the polaritonic response via frequency-dependent losses or discrete-to-continuum strong coupling. We discuss the emergent gauge freedom of the resulting theory and provide analytical expressions for all the gauge-invariant observables in both the Power–Zienau–Woolley and the Coulomb representations. In order to exemplify the results, the theory is finally specialized to two specific cases. In the first one, both light and matter resonances are characterized by Lorentzian linewidths, and in the second one, a fixed absorption band is also present. The analytical expressions derived in this paper can be used to predict, fit, and interpret results from polaritonic experiments with arbitrary values of the light–matter coupling and with losses of arbitrary intensity and spectral shape in both the light and matter channels. A Matlab code implementing our results is provided. [ABSTRACT FROM AUTHOR]

Transient absorption (TA) and time-resolved photoluminescence (PL) spectroscopies have been used to elucidate the hole tunneling and Auger dynamics in biexcitons and negative trions in high-quality InP/ZnSe/ZnS quantum dots (QDs). In a previous paper [Nguyen et al., J. Phys. Chem. C 125, 15405–15414 (2021)], we showed that under high-intensity photoexcitation, two types of biexcitons are formed: those having two conduction band electrons and two valence band holes (designated as an XX state) and those having two conduction band electrons, one valence band hole, and an additional trapped hole (designated as an XT state). In the present paper, we show that both types of biexcitons can undergo Auger processes, with those of the XT state being a factor of four to five slower than those of the XX state. In addition, the trapped holes can undergo tunneling into the valence band, converting an XT state to an XX state. The relative amplitudes of the fast (XX) and slow (XT) components are different in the TA and PL kinetics, and these differences can be quantitatively understood in terms of oscillator strengths and electron–hole overlap integrals of each state. XT to XX hole tunneling rates are obtained from the comparison of the XT state lifetimes with those of the negative trions. This comparison shows that the tunneling times decrease with decreasing core size and shell thickness. These times are about 2 ns for the thinnest shell red-emitting QDs and decrease to 330 ps for QDs that luminesce in the yellow. [ABSTRACT FROM AUTHOR]

This paper reports the energies and charge and spin distributions of both the vertically excited and fully relaxed GR1 states of the neutral singlet vacancy in diamond obtained from direct Δ-SCF calculations used previously to describe the low-lying excited states in AFII NiO and α-Al2O3. The calculations are based on the B3LYP functional in its standard form, with a C basis set that is identical to that which was used previously in numerous calculations of the ground state properties of defective diamond. Both the vertically excited and thermally relaxed GR1 states are predicted to be excitonic and insulating, with extensive re-distribution of charge and spin density and back-donation to the donor site. The present calculations suggest that the triplet state makes no contribution to the GR1 excitation. The predicted energy of the zero phonon line (1.57 eV) compares with the observed value of 1.67 eV, which also suggests that the GR1 state is neutral. The bandgaps lead to an estimate of the next higher (GR2) excited state energy, which is close to that found in the observed spectra. Similar calculations are used to predict the energies of the higher gap states at (5.0–5.5) eV, including the bulk value of 7.3 eV, which compares with the experimental value of (7.3–7.4) eV. An explanation is suggested as to why only the GR1 luminescence is observed. This paper also suggests an alternative channel for the recovery of the ground state in photoluminescence studies. [ABSTRACT FROM AUTHOR]

It is well-known experimentally that the positively charged muon and the muonium atom may bind to molecules and solids, and through muon's magnetic interaction with unpaired electrons, valuable information on the local environment surrounding the muon is deduced. Theoretical understanding of the structure and properties of resulting muonic species requires accurate and efficient quantum mechanical computational methodologies. In this paper, the two-component density functional theory (TC-DFT), as a first principles method, which treats electrons and the positive muon on an equal footing as quantum particles, is introduced and implemented computationally. The main ingredient of this theory, apart from the electronic exchange–correlation functional, is the electron–positive muon correlation functional that is foreign to the purely electronic DFT. A Wigner-type local electron–positive muon correlation functional, termed eμc-1, is proposed in this paper and its capability is demonstrated through its computational application to a benchmark set of muonic organic molecules. The TC-DFT equations containing eμc-1 are not only capable of predicting the muon's binding site correctly, but they also reproduce muon's zero-point vibrational energies and the muonic densities much more accurately than the TC-DFT equations lacking eμc-1. Thus, this study sets the stage for developing accurate electron–positive muon functionals, which can be used within the context of the TC-DFT to elucidate the intricate interaction of the positive muon with complex molecular systems. [ABSTRACT FROM AUTHOR]

We obtain the phase diagram of the hard core lattice gas with third nearest neighbor exclusion on the triangular lattice using Monte Carlo simulations that are based on a rejection-free flat histogram algorithm. In a recent paper [Darjani et al., J. Chem. Phys. 151, 104702 (2019)], it was claimed that the lattice gas with third nearest neighbor exclusion undergoes two phase transitions with increasing density with the phase at intermediate densities exhibiting hexatic order with continuously varying exponents. Although a hexatic phase is expected when the exclusion range is large, it has not been seen earlier in hard core lattice gases with short range exclusion. In this paper, by numerically determining the entropies for all densities, we show that there is only a single phase transition in the system between a low-density fluid phase and a high density ordered sublattice phase and that a hexatic phase is absent. The transition is shown to be first order in nature, and the critical parameters are determined accurately. [ABSTRACT FROM AUTHOR]

An error in the theoretical treatment by Hadinger et al., using the experimental data on HCl of Coxon and Ogilvie, is pointed out. [ABSTRACT FROM AUTHOR]

Presents a response to the commentary on the article 'Further Contributions to the Energy Levels of a Perturbed Anharmonic Oscillator: Application to Adiabatic Corrections.'

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]