The Journal of Physical Chemistry A最新文献

The internally contracted multireference coupled-cluster method with single, double and perturbative triple excitations, icMRCCSD(T), was tested for its performance in the context of computational high-accuracy thermochemistry. The results were gauged against the standard single-reference coupled-cluster hierarchy with up to 5-fold excitations. The test set comprised of a selection of first-row dinuclear compounds and the three 3d-transition metal compounds MnH, FeH, and CoH. The results revealed two problems with the current formulation of icMRCCSD(T). First, the choice of the Dyall Hamiltonian as the zeroth-order Hamiltonian, which leads to a biased description of the different orbital subspaces and particularly poor results for the atomic correlation energies, and second, the tendency to overestimate the perturbative correction for triply excited clusters, in particular in the presence of open shells and correspondingly low orbital-energy gaps. The two problems could be solved by resorting to the effective Fock operator as zeroth-order Hamiltonian and by adopting a modified amplitude equation that includes terms quadratic in the pair clusters. A similar modification was recently proposed by Masios et al. (Phys. Rev. Lett. 2023, 131, 186401) in the context of applying single-reference coupled-cluster theory to systems with small or vanishing band gaps and we chose the acronym ‘(cT*) correction’ in analogy to that work. In contrast to the work of Masios et al., additional terms including single excitation clusters were omitted, as these again lead to an overestimation of correlation effects in more difficult cases. We also tested another alternative for the zeroth-order Hamiltonian and additional higher-order corrections for the correlation energy. These extensions did not significantly improve the results and were also computationally more demanding. The improved icMRCCSD(cT*)_F method yields very accurate results with errors, relative to accurate benchmarks, better than 2 kJ/mol for total energies and atomization energies for the entire set of examples considered in this work.

We report a high-resolution crossed molecular beam experiment investigating the reaction dynamics of the F + CH₄ → HF + CH₃ reaction across a broad range of collision energies (3.1–13.8 kcal/mol). By using time-sliced velocity map imaging within the crossed molecular beam apparatus, we obtain correlated angular distributions and branching ratios for various product pairs (CH₃(ν), HF(ν′)). The resolved reactive rainbow-like features that display a distinct bulge in the angular distribution in different channels reveal the formation of vibrationally ground and excited states of the methyl radical accompanied by different rovibrationally excited HF products. These findings suggest distinct reaction dynamics for channels leading to vibrationally excited methyl radicals compared to those forming ground-state methyl radicals.

Halogenated volatile organic compounds (HVOCs) pose significant bioaccumulative and toxicological risks, necessitating effective strategies for their removal. Here, we show, through a computational study employing density functional theory and coupled cluster methods, the detailed mechanism and kinetic properties of Cl-initiated degradation reactions of 2-chloropropane (2-CP, (CH₃)₂CHCl) and 2-methylpropanoyl halide ((CH₃)₂CHCOX, X = Cl, Br, F). The reaction rate constants of all the channels were calculated by the canonical variational transition state theory (CVT) with the correction of the small curvature tunneling effect (SCT) at 200–1000 K. The subsequent transformation pathways of the major radical products of (CH₃)₂CHCl and (CH₃)₂CHCOCl in the presence of O₂, NO, and HO₂ radical were investigated. The results elucidate the reaction pathways and rate constants, which are in excellent agreement with the experimental data at 296 K. We further explore the atmospheric implications of these reactions by assessing the atmospheric lifetime (τ) and ozone depletion potential (ODP). Additionally, we delve into the aquatic toxicity and bioaccumulation potential of the reactants and their transformation products. This study not only advances our knowledge of the atmospheric fate of halogenated hydrocarbons but also underscores the importance of considering the environmental and toxicological impacts in the development of HVOC mitigation strategies.

Machine learning (ML) potentials such as the Gaussian approximation potential (GAP) have demonstrated impressive capabilities in mapping structure to properties across diverse systems. Here, we introduce a GAP model for low-dimensional Ni nanoclusters and demonstrate its flexibility and effectiveness in capturing the energetics, structural diversity, and thermodynamic properties of Ni nanoclusters across a broad size range. Through a systematic approach encompassing model development, validation, and application, we evaluate the model's efficacy in representing energetics and configurational features in low-dimensional regimes while also examining its extrapolative nature to vastly different spatiotemporal regimes. Our analysis and discussion shed light on the data quality required to effectively train such models. Trajectories from large-scale MD simulations using the GAP model analyzed with data-driven models like graph neural networks reveal intriguing insights into the size-dependent phase behavior and thermomechanical stability characteristics of porous Ni nanoparticles. Overall, our work underscores the potential of ML models, which coupled with data-driven approaches serve as versatile tools for studying low-dimensional systems and complex material dynamics.

Liquid-jet photoelectron spectroscopy (LJ-PES) and electronic-structure theory were employed to investigate the chemical and structural properties of the amino acid l-proline in aqueous solution for its three ionized states (protonated, zwitterionic, and deprotonated). This is the first PES study of this amino acid in its biologically relevant environment. Proline's structure in the aqueous phase under neutral conditions is zwitterionic, distinctly different from the nonionic neutral form in the gas phase. By analyzing the carbon 1s and nitrogen 1s core levels as well as the valence spectra of aqueous-phase proline, we found that the electronic structure is dominated by the protonation state of each constituent molecular site (the carboxyl and amine groups) with small yet noticeable interference across the molecule. The site-specific nature of the core-level spectra enables the probing of individual molecular constituents. The valence photoelectron spectra are more difficult to interpret because of the overlapping signals of proline with the solvent and pH-adjusting agents (HCl and NaOH). Yet, we are able to reveal subtle effects of specific (hydrogen-bonding) interaction with the solvent on the electronic structure. We also demonstrate that the relevant conformational space is much smaller for aqueous-phase proline than for its gas-phase analogue. This study suggests that caution must be taken when comparing photoelectron spectra for gaseous- and aqueous-phase molecules, particularly if those molecules are readily protonated/deprotonated in solution.

Quantitative characterization of the surface of microdroplets is important to understanding and predicting numerous chemical and physical processes, such as cloud droplet formation and accelerated chemistry in microdroplets. However, it is increasingly appreciated that the surface compositions of microdroplets do not necessarily match those of macroscale solution due to their large surface-area-to-volume (SA–V) ratios and confined volumes. In this work, we explore how both droplet size and composition affect the surface composition of microdroplets by measuring the equilibrium surface tensions of levitated microdroplets containing a single surfactant. We measure the critical micelle concentrations (CMCs) for surfactants of various strengths (macroscale CMC values ranging from 0.02 to 10 mM) in microdroplets with radii ranging from 5 to 25 μm. We accurately model the surface tensions of microdroplets using an equilibrium partitioning model that only requires droplet size and adsorption parameters from macroscale measurements as inputs. Our model predicts that surfactants have an “effective CMC” in microdroplets that is always larger in value than the corresponding macroscale CMC. In some instances, the effective CMC of a surfactant in microdroplets is several orders of magnitude larger than both its macroscale CMC and its macroscale solubility limit. We present a simple expression for the effective CMC in microdroplets that depends on both the macroscale CMC of a surfactant and the SA–V ratio of the microdroplet. Ultimately, our experimental results and model can be used broadly to predict microdroplet surface compositions when investigating surface-driven accelerated chemistry in microdroplets or estimating cloud droplet activation.

The thermal unimolecular decay of ethoxy is important in high-temperature combustion environments where the ethoxy radical is a key reactive intermediate. Two dissociation pathways of ethoxy, including the β-C–C scission to yield CH₃ + CH₂O and the H-elimination to make H + CH₃CHO, were characterized using a high-level coupled-cluster-based composite quantum chemical method (mHEAT-345(Q_Λ)). The former route is found to be dominant while the latter is insignificant, in agreement with previous experimental and theoretical studies. Thermal rate coefficients are calculated for P = 0.001–658 atm (of air) and T = 300–2500 K using semiclassical transition state theory (SCTST) in combination with a pragmatic two-dimensional E,J-resolved master equation (2DME). The effects of tunneling and anharmonicity on the calculated rate constants are also examined. The tunneling factor is found to be inversely dependent on pressure, contrary to previous observations of pressure-dependent tunneling in entrance channels.

I and Br atoms are placed on opposite ends of a n-butyl group, with each allowed to form a halogen bond (XB) with NH₃. DFT calculations show that the intrinsic preference of the nucleophile for the heavier I over Br can be reversed by the proper placement of substituents on the alkyl chain. A similar reversal occurs for NH₂ and OH groups on the alkyl chain, where substituents make the O a better electron donor than N in an XB to an electrophilic ICCH. The highly mobile π-electron cloud of an aromatic ring makes such reversals much more difficult when the pair of competing atoms are placed on, or within, such a ring.

Transition paths refer to rare events in physics, chemistry, and biology where the molecules cross barriers separating stable molecular conformations. The conventional analysis of the transition path times employs a diffusive and memoryless transition over a smooth potential barrier. However, it is widely acknowledged that the free energy profile between two minima in biomolecular processes is inherently not smooth. In this article, we discuss a theoretical model with a parabolic rough potential barrier and obtain analytical results of the transition path distribution and mean transition path times by incorporating absorbing boundary conditions across the boundaries under the driving of Gaussian white noise. Further, the influence of anomalous dynamics in rough potential driven by a power-law memory kernel is analyzed by deriving a time-dependent scaled diffusion coefficient that coarse-grains the effects of roughness, and the system’s dynamics is reduced to a scaled diffusion on a smooth potential. Our theoretical results are tested and validated against numerical simulations. The findings of our study show the influence of the boundary conditions, barrier height, barrier roughness, and memory effect on the transition path time distributions in a rough potential, and the validity of the scaling diffusion coefficient has been discussed.

When trying to reach convergence of quantum chemical calculations toward the complete basis set limit, crystalline solids generally prove to be more challenging than molecules. This is due both to the closer packing of atoms─hence, to linear dependencies─and to the problematic behavior of Ewald techniques used for dealing with the infinite character of Coulomb sums. Thus, a dual basis set approach is even more desirable for periodic systems than for molecules. In such an approach, the self-consistent procedure is implemented in a small basis set, and the effect of the enlargement of the basis set is estimated a posteriori. In this paper, we extend to crystalline solids our previous coupled perturbed dual basis set approach [J. Chem. Theory Comput. 2020, 16, 1, 340–353] in which the basis set enlargement is treated as a perturbation. Among the notable features of this approach are (i) the possibility of obtaining not only a correction to the energy but also to energy bands and electron density; (ii) the absence of a diagonalization step for the full Fock matrix in the large basis set; and (iii) the possibility of extrapolating low order perturbation energy corrections to infinite order. We also present here the first periodic implementation of the dual basis set method of Liang and Head-Gordon [J. Phys. Chem. A 2004, 108, 3206–3210]. The effectiveness of both approaches is, then, compared on a small, but representative, set of solids.