Search

Nonadiabatic quantum-classical mapping approaches have significantly gained in popularity over the past several decades because they have acceptable accuracy while remaining numerically tractable even for large system sizes. In the recent few years, several novel mapping approaches have been developed that display higher accuracy than the traditional Ehrenfest method, linearized semiclassical initial value representation (LSC-IVR), and Poisson bracket mapping equation (PBME) approaches. While various benchmarks have already demonstrated the advantages and limitations of those methods, unified theoretical justifications of their short-time accuracy are still demanded. In this article, we systematically examine the intra-electron correlation, as a statistical measure of electronic phase space, which has been first formally proposed for mapping approaches in the context of the generalized discrete truncated Wigner approximation and which is a key ingredient for the improvement in short-time accuracy of such mapping approaches. We rigorously establish the connection between short-time accuracy and intra-electron correlation for various widely used models. We find that LSC-IVR, PBME, and Ehrenfest methods fail to correctly reproduce the intra-electron correlation. While some of the traceless Meyer-Miller-Stock-Thoss (MMST) approaches, partially linearized density matrix (PLDM) approach, and spin partially linearized density matrix (spin-PLDM) approach are able to sample the intra-electron correlation correctly, the spin linearized semiclassical (spin-LSC) approach, which is a specific example of the classical mapping model, and the other traceless MMST approaches sample the intra-correlation faithfully only for two-level systems. Our theoretical analysis provides insights into the short-time accuracy of semiclassical methods and presents mathematical justifications for previous numerical benchmarks., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

This research paper presents a new fundamental approach for evaluating accurate ab initio quartic, sextic, and octic centrifugal distortion parameters of A-reduced rotational effective Hamiltonians of asymmetric top molecules. In this framework, the original Watson Hamiltonian, expanded up to sextic terms of kinetic and potential energies, is subjected to a series of vibrational and rotational operator unitary transformations, leading to reduced Watson effective Hamiltonians for the equilibrium configuration, ground state, and weakly perturbed vibrationally excited states. The proposed scheme is based on a numerical-analytic implementation of the sixth-order Van Vleck operator perturbation theory with the systematic normal ordering of vibrational rising and lowering operators (a†, a) and cylindrical angular momentum operators (Jz, J+, J-). The efficiency of the developed theoretical model is demonstrated by the juxtaposition of predicted centrifugal distortion parameters for several three to eight atomic molecules, including H2S, CH2O, C2H4, CH2D2, CH2F2, CH2Cl2, and B2H6, using the coupled-cluster single double triple/quadruple-ζ level of quantum chemistry. In comparison with the values derived using the customary analytic expressions, the calculated quartic and sextic parameters may improve by an order of magnitude in the fourth and sixth orders, respectively, reaching an accuracy of about 1%. Predicted octic constants can serve as an excellent starting point for fitting to experimental spectra., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

We recorded and analyzed the microwave spectra of 2,5-dimethylanisole using a pulsed molecular jet Fourier transform microwave spectrometer and the newly developed Passage And Resonant-Impulse Synergy spectrometer across a frequency range of 2-26.5 GHz with support from quantum chemical calculations. Only one conformer was predicted and observed, where the methoxy group and its adjacent methyl group adopt anti-positions. The two methyl groups, located at the ortho- and meta-positions of the anisole ring, exhibit internal rotation, resulting in quintet splitting of all rotational transitions. The low torsional barrier of the m-methyl group, amounting to 65.723 611(84) cm-1, combined with the intermediate barrier of 451.664(19) cm-1 for the o-methyl group, presented challenges in the spectral analysis. Using the XIAM and BELGI-Cs-2Tops programs, we successfully fitted 460 torsional intrastate rotational transitions, allowing for precise determination of molecular parameters and internal rotation characteristics. The torsional barriers are compared to those in the related isomers 2,3-, 2,4-, and 3,4-dimethylanisole as well as other o- and m-substituted toluene derivatives., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

We explore the effect of Couette flow on knotted linear polymer chains with extensive molecular dynamics simulations. Hydrodynamic interactions are accounted for using multi-particle collision dynamics. The polymer chain, originally containing a simple trefoil knot at rest, is described by a coarse-grained bead-spring model in a coil or globular state. We demonstrate that under shear existing loosely localized knots in polymer coils typically tighten to several segments beyond a certain shear rate threshold. At large shear rates, the polymer undergoes a tumbling-like motion during which knot sizes can fluctuate. In contrast, sheared knotted globules unwind into a convoluted pearl-necklace structure of sub-globules that folds back onto itself and in which knot types change over time., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

Critical Casimir forces emerge among particles or surfaces immersed in a near-critical fluid, with the sign of the force determined by surface properties and with its strength tunable by minute temperature changes. Here, we show how such forces can be used to trap a colloidal particle and levitate it above a substrate with a bull's-eye pattern consisting of a ring with surface properties opposite to the rest of the substrate. Using the Derjaguin approximation and mean-field calculations, we find a rich behavior of spherical colloids at such a patterned surface, including sedimentation toward the ring and levitation above the ring (ring levitation) or above the bull's-eye's center (point levitation). Within the Derjaguin approximation, we calculate a levitation diagram for point levitation showing the depth of the trapping potential and the height at which the colloid levitates, both depending on the pattern properties, the colloid size, and the solution temperature. Our calculations reveal that the parameter space associated with point levitation shrinks if the system is driven away from a critical point, while, surprisingly, the trapping force becomes stronger. We discuss the application of critical Casimir levitation for sorting colloids by size and for determining the thermodynamic distance to criticality. Our results show that critical Casimir forces provide rich opportunities for controlling the behavior of colloidal particles at patterned surfaces., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

This research elucidates the intricate nature of electronic coupling in the redox-active (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO), commonly utilized in organic radical batteries. This study employs a combination of classical molecular dynamics and various electronic coupling calculation schemes. Within the context of the generalized Mulliken-Hush method, the electronic couplings are investigated via the complete active space self-consistent field approach, in combination with n-electron valence state perturbation theory, to provide an accurate description of both static and dynamic electron correlation as well as using (time-dependent) density functional theory simulations. Furthermore, the electronic communication between redox-active sites is studied using the cost-efficient density functional theory (DFT)-based frontier molecular orbital (FMO) approach. Our study reveals the dependence of the electronic coupling on the distance and the relative orientation of the redox pairs (TEMPO and TEMPO+). Apart from the expected exponential distance dependence, we found pronounced orientation dependence, with coupling values varying up to 0.2 eV, which is reflected by a substantial basis set dependency of the couplings, in particular at short distances. In addition, our study highlights the limitations of the DFT-based FMO method, in particular at short intermolecular distances between the redox-active sites, which may lead to a mixing of the involved molecular orbitals. This comparison will provide us with the most cost-accuracy-effective method for calculating electronic couplings in TEMPO-TEMPO+ systems., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

When we place conducting bodies in electrolyte solutions, their surface potential Φs appears to be much smaller in magnitude than the applied one Φ0 and normally does not obey the classical electrostatic boundary condition of a constant potential expected for conductors. In this paper, we demonstrate that an explanation of these observations can be obtained by postulating that diffuse ions condense at the "wall" due to the reduced permittivity of a solvent. For small values of Φ0, the surface potential responds linearly. On increasing Φ0 further, Φs augments nonlinearly and then saturates to a constant value. Analytical approximations for Φs derived for these three distinct modes show that it always adjusts to salt concentration, which is equivalent to a violation of the constant potential condition. The latter would be appropriate for highly dilute solutions but only if Φ0 is small. Surprisingly, when the plateau with high Φs is reached, the conductor surface switches to a constant charge density condition normally expected for insulators. Our results are directly relevant for conducting electrodes, mercury drops, colloidal metallic particles, and more., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

This study investigates the translational and rovibrational energy of vacuum-ultraviolet (VUV) photodesorbed CO molecules from a CO polycrystalline ice (15 K) at ∼8 eV. The electronic excitation was produced by a pulsed VUV laser, and the photodesorption of CO molecules in their ground and first vibrational states was observed using resonance enhanced multiphoton ionization. Time-of-flight and rotationally resolved spectra were measured, and the kinetic and internal energy distribution were obtained. Vibrationally cold CO molecules were observed, with little energy in rotation and translation (≤300 meV). Ab Initio Molecular Dynamics (AIMD) simulations focusing on the description of the vibrational energy redistribution within an aggregate of 50 CO molecules were performed. The measured energy distributions are in very good agreement with those predicted by AIMD simulations. The rotational energy was found to slightly increase with translational energy, a trend also predicted by theory. This confirms the validity of the indirect desorption mechanism triggered by the excitation of CO in a high vibrational state., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

Water interfaces play a crucial role in regulating interactions and energy flow. Vibrational sum-frequency generation (vSFG) spectroscopy provides structural and dynamic information on water molecules at interfaces. It has revealed, for instance, the presence of the hydrogen-bonded and free OH groups at the air-water interface. Here, using temperature-dependent, time-resolved vSFG, we focus on the vibrational energy relaxation dynamics of interfacial heavy water (D2O). We reveal that while the relaxation timescale for hydrogen-bonded OD stretch modes is temperature-independent, the lifetime of the free OD stretch mode decreases with increasing temperature. Our data, supported by simulations, suggest that both intramolecular energy transfer and rotational reorientation mechanisms jointly contribute to the energy relaxation process of the free OD, with temperature influencing these mechanisms differently., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

Methane (CH4) oxidation is an important reaction to reduce the greenhouse effect caused by incomplete combustion of CH4. Here, we explored the mechanism of CH4 oxidation catalyzed by CeO2 and Ni-doped CeO2, focusing on the redox properties of these catalyst surfaces, using density functional theory (DFT). We found that the barriers for CH4* activation and H2O* formation are correlated with the surface redox capacity, which is enhanced by Ni doping. Furthermore, the complete reaction mechanism is explored by DFT calculations and microkinetic simulations on bare and Ni-doped CeO2 surfaces. Our calculations suggest that the doping of Ni leads to a much higher overall reactivity, due to a balance between the CH4* activation and H2O* formation steps. These results provide insights into the CH4 oxidation mechanism and the intrinsic relationship between redox properties and the activity of CeO2 surfaces., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

Unlike with the energy, which is a scalar property, machine learning (ML) prediction of vector or tensor properties poses the additional challenge of achieving proper invariance (covariance) with respect to molecular rotation. For the energy gradients needed in molecular dynamics (MD), this symmetry is automatically fulfilled when taking analytic derivative of the energy, which is a scalar invariant (using properly invariant molecular descriptors). However, if the properties cannot be obtained by differentiation, other appropriate methods should be applied to retain the covariance. Several approaches have been suggested to properly treat this issue. For nonadiabatic couplings and polarizabilities, for example, it was possible to construct virtual quantities from which the above tensorial properties are obtained by differentiation and thus guarantee the covariance. Another possible solution is to build the rotational equivariance into the design of a neural network employed in the model. Here, we propose a simpler alternative technique, which does not require construction of auxiliary properties or application of special equivariant ML techniques. We suggest a three-step approach, using the molecular tensor of inertia. In the first step, the molecule is rotated using the eigenvectors of this tensor to its principal axes. In the second step, the ML procedure predicts the vector property relative to this orientation, based on a training set where all vector properties were in this same coordinate system. As the third step, it remains to transform the ML estimate of the vector property back to the original orientation. This rotate-predict-rotate (RPR) procedure should thus guarantee proper covariance of a vector property and is trivially extensible also to tensors such as polarizability. The RPR procedure has an advantage that the accurate models can be trained very fast for thousands of molecular configurations, which might be beneficial where many training sets are required (e.g., in active learning). We have implemented the RPR technique, using the MLatom and Newton-X programs for ML and MD, and performed its assessment on the dipole moment along MD trajectories of 1,2-dichloroethane., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

We investigate attachment of slow electrons (0-10 eV) to naphthalene (Np) clusters in a crossed beam experiment. Supersonic expansions under different conditions using different buffer gases generate the clusters: in He, Ne, and low pressure Ar, neat (Np)N clusters are formed, while we also observe mixed clusters of naphthalene with rare-gas atoms in co-expansion with Ar above 0.5 bar and with Kr. Negatively charged (Np)n- and Rgm(Np)n- (Rg = Ar, Kr) clusters are analyzed by mass spectrometry, and electron energy dependent ion yields are measured. We show that the smallest stable naphthalene complex with an excess electron, the dimer (Np)2- anion, cannot be formed in a binary collision of a free electron with (Np)2 dimer, nor with (Np)3 trimer. Evaporation of a weakly bound Ar atom(s) from a mixed ArM(Np)2 cluster following electron attachment leads to the dimer (Np)2- anion. Larger (Np)n-, n > 3, transient cluster anions decay via evaporation of an Np unit on a timescale of tens of microseconds. The self-scavenging process opens around 6 eV, where a naphthalene unit is electronically excited by the incoming electron, which is slowed down and trapped. However, the transient negative ion is efficiently stabilized only in the mixed clusters, from which Ar atom(s) can be evaporated., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

When encountering complex fluorescence decays that deviate from exponentiality, a very appealing approach is to use lifetime or rate constant distributions. These are related by Laplace transform to the sum of exponential functions, stretched exponentials, Becquerel's decay function, and others. However, the limitations of this approach have not been sufficiently discussed in the literature. In particular, the time-independent probability distributions of the rate constants or decay times are occasionally used to describe bimolecular quenching. We show that in such a case, this mathematical formalism has a clear physical interpretation only when the fluorophore and quencher molecules are immobile, as in the solid state. However, such an interpretation is no longer possible once we consider the motion of fluorophores with respect to quenchers. Therefore, for systems in which the relative motion of fluorophores and quenchers cannot be neglected, it is not appropriate to use the time-independent rate or decay time distributions to describe, fit, or rationalize experimental results on fluorescence decay., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

Simulating electrolyte-electrode systems poses challenges due to the need to account for the electrode's response to ion movements in order to maintain a constant electrode potential, which slows down the simulations. To circumvent this, computationally more efficient constant charge (CC) simulations are sometimes employed. However, the accuracy of CC simulations in capturing the behavior of electrolyte-electrode systems remains unclear, especially for microporous electrodes. Herein, we consider electrolyte-filled slit nanopores and systematically analyze the in-pore ion structure and diffusivity using CC and constant potential simulations. Our results indicate that CC simulations provide comparable pore occupancies at high bulk ion densities and for highly charged pores, but they fail to accurately describe the ion structure and dynamics, particularly in quasi-2D (single-layer) pores and at low ion densities. We attribute these results to the superionic state emerging in conducting nanoconfinement and its interplay with excluded volume interactions., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

The hydrodynamics of thin films is typically described using macroscopic models whose connection to the microscopic particle dynamics is a subject of ongoing research. Existing methods based on density functional theory provide a good description of static thin films but are not sufficient for understanding nonequilibrium dynamics. In this work, we present a microscopic derivation of the thin film equation using the Mori-Zwanzig projection operator formalism. This method allows to directly obtain the correct gradient dynamics structure along with microscopic expressions for mobility and free energy. Our results are verified against molecular dynamics simulations for both simple fluids and polymers., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

The diffusion of cations in organic solvent solutions is important for the performance of metal-ion batteries. In this article, pulsed field gradient nuclear magnetic resonance experiments and fully atomistic molecular dynamic simulations were employed to study the temperature-dependent diffusive behavior of various liquid electrolytes representing 1M propylene carbonate solutions of metal salts with bis(trifluoromethylsulfonyl)imide (TFSI-) or hexafluorophosphate (PF6-) anions commonly used in lithium-ion batteries and beyond. The experimental studies revealed the temperature dependence of the diffusion coefficients for the propylene carbonate (PC) solvent and for the anions following an Arrhenius type of behavior. It was observed that the PC molecules are the faster species. For the monovalent cations (Li+, Na+, K+), the PC solvent diffusion was enhanced as the cation size increased, while for the divalent cations (Mg2+, Ca2+, Sr2+, Ba2+), the opposite trend was observed, i.e., the diffusion coefficients decreased as the cation size increased. The anion diffusion in LiTFSI and NaTFSI solutions was found to be similar, while in electrolytes with divalent cations, a decrease in anion diffusion with increasing cation size was observed. It was shown that non-polarizable charge-scaled force fields could correspond perfectly to the experimental values of the anion and PC solvent diffusion coefficients in salt solutions of both monovalent (Li+, Na+, K+) and divalent (Mg2+, Ca2+, Sr2+, Ba2+) cations at a range of operational temperatures. Finally, after calculating the radial distribution functions between cations, anions, and solvent molecules, the increase in the PC diffusion coefficient established with the increase in cation size for monovalent cations was clearly explained by the large hydration shell of small Li+ cations, due to their strong interaction with the PC solvent. In solutions with larger monovalent cations, such as Na+, and with a smaller solvation shell of PC, the PC diffusion is faster due to more liberated solvent molecules. In the salt solutions with divalent cations, both the anion and the PC diffusion coefficients decreased as the cation size increased due to an enhanced cation-anion coordination, which was accompanied by an increase in the amount of PC in the cation solvation shell due to the presence of anions., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

Some of the authors of the present research group have previously reported mass spectrometric detection of [PdAu9(PPh3)8(CN)]2+ (PdAu9CN) by atmospheric pressure plasma (APP) irradiation of [MAu8(PPh3)8]2+ (PdAu8) in methanol and proposed based on density functional theory (DFT) calculations that PdAu9CN is constructed by inserting a CNAu or NCAu unit into the Au-PPh3 bond of PdAu8 [Emori et al., J. Chem. Phys. 155, 124312 (2021)]. In this follow-up study, we revisited the structure of PdAu9CN by high-resolution ion mobility spectrometry on an isolated sample of PdAu9CN with the help of dispersion-corrected DFT calculation. In contradiction to the previous proposal, we conclude that isomers in which an AuCN unit is directly bonded to the central Pd atom of PdAu8 are better candidates. This assignment was supported by Fourier transform infrared and ultraviolet-visible spectroscopies of isolated PdAu9CN. The simultaneous formation of [Au(PPh3)2]+ and PdAu9CN suggests that the AuCN species are formed by APP irradiation at the expense of a portion of PdAu8. These results indicate that APP may offer a unique method for transforming metal clusters into novel ones by generating in situ active species that were not originally added to the solution., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

Electrocatalytic CO2 reduction reaction (CO2RR) for CH4 production presents a promising strategy to address carbon neutrality, and the incorporation of a second metal has been proven effective in enhancing catalyst performance. Nevertheless, there remains limited comprehension regarding the fundamental factors responsible for the improved performance. Herein, the critical role of Pd in electrocatalytic CO2 reduction to CH4 on Cu-based catalysts has been revealed at a molecular level using in situ surface-enhanced Raman spectroscopy (SERS). A "borrowing" SERS strategy has been developed by depositing Cu-Pd overlayers on plasmonic Au nanoparticles to achieve the in situ monitoring of the dynamic change of the intermediate during CO2RR. Electrochemical tests demonstrate that Pd incorporation significantly enhances selectivity toward CH4 production, and the Faradaic efficiency (FE) of CH4 is more than two times higher than that for the catalysts without Pd. The key intermediates, including *CO2-, *CO, and *OH, have been directly identified under CO2RR conditions, and their evolution with the electrochemical environments has been determined. It is found that Pd incorporation promotes the activation of both CO2 and H2O molecules and accelerates the formation of abundant active *CO and hydrogen species, thus enhancing the CH4 selectivity. This work offers fundamental insights into the understanding of the molecular mechanism of CO2RR and opens up possibilities for designing more efficient electrocatalysts., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

We present a simple approximation to estimate the largest charge that a given molecule can hold until fragmentation into smaller charged species becomes more energetically favorable. This approximation solely relies on the ionization potentials, electron affinities of the parent and fragment species, and also on the neutral parent's dissociation energy. By parameterizing these quantities, it is possible to obtain analytical phase diagrams of polycationic stability. We demonstrate the applicability of this approach by discussing the maximal charge dependence on the size of the molecular system. A numerical demonstration for linear polyenes, monocyclic annulenes, and helium clusters is provided., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

Ice XIX and ice XV are both partly hydrogen-ordered counterparts to disordered ice VI. The ice XIX → XV transition represents the only order-to-order transition in ice physics. Using Raman and dielectric spectroscopies, we investigate the ambient-pressure kinetics of the two individual steps in this transition in real time (of hours), that is, ice XIX → transient ice VI (the latter called VI‡) and ice VI‡ → ice XV. Hydrogen-disordered ice VI‡ appears intermittent between 101 and 120 K, as inferred from the appearance and subsequent disappearance of the ice VI Raman marker bands. A comparison of the rate constants for the H2O ices reported here with those from D2O samples [Thoeny et al., J. Chem. Phys. 156, 154507 (2022)] reveals a large kinetic isotope effect for the ice XIX decay, but a much smaller one for the ice XV buildup. An enhancement of the classical overbarrier rate through quantum tunneling for the former can provide a possible explanation for this finding. The activation barriers for both transitions are in the 18-24 kJ/mol range, which corresponds to the energy required to break a single hydrogen bond. These barriers do not show an H/D isotope effect and are the same, no matter whether they are derived from Raman scattering or from dielectric spectroscopy. These findings favor the notion that a dipolar reorientation, involving the breakage of a hydrogen bond, is the rate determining step at the order-to-order transition., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

We have investigated pressure-induced amorphization (PIA) of an alcohol clathrate hydrate (CH) of cubic structure type I (sI) in the presence of NH4F utilizing dilatometry and x-ray powder diffraction. PIA occurs at 0.98 GPa at 77 K, which is at a much lower pressure than for other CHs of the same structure type. The amorphized CH also shows remarkable resistance against crystallization upon decompression. While amorphized sI CHs could not be recovered previously at all, this is possible in the present case. By contrast to other CHs, the recovery of the amorphized CHs to ambient pressure does not even require a high-pressure annealing step, where recovery without any loss of amorphicity is possible at 120 K and below. Furthermore, PIA is accessible upon compression at unusually high temperatures of up to 140 K, where it reaches the highest degree of amorphicity. Molecular dynamics simulations confirm that polar alcoholic guests, as opposed to non-polar guests, induce cage deformation at lower pressure. The substitution of NH4F into the host-lattice stabilizes the collapsed state more than the crystalline state, thereby enhancing the collapse kinetics and lowering the pressure of collapse., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

We determine the surface tension of aqueous electrolyte solutions in contact with non-polar dielectric media using a thermomechanical approach, which involves deriving the stress tensor from the thermodynamic potential of an inhomogeneous fluid. To obtain the surface tension, we calculate both the normal and tangential pressures using the components of the stress tensor, recently derived by us [Y. A. Budkov and P. E. Brandyshev, J. Chem. Phys. 159, 174103 (2023)] within the framework of Wang's variational field theory. Using this approach, we derive an analytical expression for the surface tension in the linear approximation. At low ionic concentrations, this expression represents the classical Onsager-Samaras limiting law. By utilizing only one fitting parameter, which is related to the affinity of anions to the dielectric boundary, we successfully approximated experimental data on the surface tension of several aqueous electrolyte solutions. This approximation applies to both the solution-air and solution-dodecane interfaces, covering a wide range of electrolyte concentrations., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

The pervasive presence of nanoplastics (NPs) in the environment has gained increasing attention due to their accumulation in living organisms. These emerging contaminants inevitably interact with extracellular polymeric substances along respiratory or gastrointestinal tracts, and diverse organic coating on the surface of NPs, known as bio- or eco-corona, is formed. Although its impact on altering the NP properties and potential cell internalization has been extensively examined, studies on its role in NP partitioning in the cell membrane are elusive yet. In this work, molecular dynamics is used to investigate the formation of chitosan (CT) corona centered on a polyvinyl chloride (PVC) nanoparticle and the uptake of the resulting complex onto lipid membranes. Coarse-grained models compatible with the newly developed Martini 3.0 force field are implemented for the two polymers employing the atomistic properties as targets in the parameterization. The reliability of the coarse-grained polymer models is demonstrated by reproducing the structural properties of the PVC melt and of solvated CT strands, as well as by determining the conformation adopted by the latter at the NP surface. Results show that the spontaneous binding of CT chains of high and intermediate protonation degrees led to the formation of soft and hard corona that modulates the interaction of PVC core with model membranes. The structural changes of the corona adsorbed at the lipid-water interface enable a subsequent transfer of the NP to the center of the saturated lipid membranes and a complete or partial transition to a snorkel conformation depending on the hydrophilic/hydrophobic balance in the CT-PVC complex. Overall, the computational investigation of the coarse-grained model system provides implications for understanding how the eco-corona development influences the uptake and implicit toxicology of NPs., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

A nitrogen K-edge x-ray absorption near-edge structure (XANES) survey is presented for tetrapyrido[3,2-a:2',3'-c:3″,2″-h:2‴,3‴-j]phenazine (tpphz)-bridged bimetallic assemblies that couple chromophore and catalyst transition metal complexes for light driven catalysis, as well as their individual molecular constituents. We demonstrate the high N site sensitivity of the N pre-edge XANES features, which are energetically well-separated for the phenazine bridge N atoms and for the individual metal-bound N atoms of the inner coordination sphere ligands. By comparison with the time-dependent density functional theory calculated spectra, we determine the origins of these distinguishable spectral features. We find that metal coordination generates large shifts toward higher energy for the metal-bound N atoms, with increasing shift for 3d < 4d < 5d metal bonding. This is attributed to increasing ligand-to-metal σ donation that increases the effective charge of the bound N atoms and stabilizes the N 1s core electrons. In contrast, the phenazine bridge N pre-edge peak is found at a lower energy due to stabilization of the low energy electron accepting orbital localized on the phenazine motif. While no sensitivity to ground state electronic coupling between the individual molecular subunits was observed, the spectra are sensitive to structural distortions of the tpphz bridge. These results demonstrate N K-edge XANES as a local probe of electronic structure in large bridging ligand motifs, able to distinctly investigate the ligand-centered orbitals involved in metal-to-ligand and ligand-to-ligand electron transfer following light absorption., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

Bromine atom (Br) reactions lead to ozone depletion in the troposphere and stratosphere. Photodegradation of bromocarbons is one of the main sources of bromine atoms in the atmosphere. Here, we use high-level ab initio methods, including spin-orbit effects, to study the photodissociation of the CH2Br radical. All possible fragmentation pathways, namely CH2Br + hν → CH2 + Br, HCBr + H, and CBr + H2, have been analyzed. Potential-energy curves of the ground and several excited electronic states along the corresponding dissociating bond distance of each pathway have been calculated. Considering the actinic fluxes of solar irradiation in the troposphere and in the stratosphere in the relevant range of frequencies, it is found that the first five excited states of CH2Br can be accessed from the ground state. Analysis of the potential curves shows that the pathways producing CH2 + Br and HCBr + H can proceed through a fast direct dissociation mechanism, while the pathway leading to CBr + H2 involves much slower dissociation mechanisms like internal conversion between electronic states, predissociation, or tunneling through exit barriers. The main implications are that the two faster channels are predicted to be dominant, and the slower pathway is expected to be less relevant. The tropospheric and stratospheric solar actinic fluxes also allow for further dissociation of the HCBr and CBr fragments, generating additional Br atoms, provided that they survive possible collisions with other atmospheric reagents. Finally, we discuss the possible effect of each of the three CH2Br dissociation pathways on the depletion of atmospheric ozone., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

A route to assess non-linear light-matter interactions from the increasingly popular GW-Bethe-Salpeter equation (GW-BSE) method is outlined. In the present work, the necessary analytic expressions within the static-screened exchange approximation of the BSE are derived. This enables a straightforward implementation of the computation of the first hyperpolarizability as well as two-photon absorption processes for molecular systems. Benchmark calculations on small molecular systems reveal that the GW-BSE method is intriguingly accurate for predicting both first hyperpolarizabilities and two-photon absorption strengths. Using state-of-the-art Kohn-Sham references as a starting point, the accuracy of the GW-BSE method rivals that of the coupled-cluster singles-and-doubles method, outperforming both second-order coupled-cluster and time-dependent density-functional theory., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

In the past decade, the quantum chemical version of the density matrix renormalization group method has established itself as the method of choice for strongly correlated molecular systems. However, despite its favorable scaling, in practice, it is not suitable for computations of dynamic correlation. Several approaches to include that in post-DMRG methods exist; in our group, we focused on the tailored coupled cluster (TCC) approach. This method works well in many situations; however, in exactly degenerate cases (with two or more determinants of equal weight), it exhibits a bias toward the reference determinant representing the Fermi vacuum. Although sometimes it is possible to use a compensation scheme to avoid this bias for energy differences, it is certainly a drawback. In order to overcome this bias of the TCC method, we have developed a Hilbert-space multireference version of tailored CC, which can treat several determinants on an equal footing. We have implemented and compared the performance of three Hilbert-space multireference coupled cluster (MRCC) variants-the state universal one and the Brillouin-Wigner and Mukherjee's state specific ones. We have assessed these approaches on the cyclobutadiene and tetramethyleneethane molecules, which are both diradicals with exactly degenerate determinants at a certain geometry. We have also investigated the sensitivity of the results on the orbital rotation of the highest occupied and lowest unoccupied molecular orbital (HOMO-LUMO) pair, as it is well known that Hilbert-space MRCC methods are not invariant to such transformations., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

Equilibrium thermodynamics describes the energy exchange of a body with its environment. Here, we describe the global energy exchange of an ideal gas in the Coutte flow in a thermodynamic-like manner. We derive a fundamental relation between internal energy as a function of parameters of state. We analyze a non-equilibrium transition in the system and postulate the extremum principle, which determines stable steady states in the system. The steady-state thermodynamic framework resembles equilibrium thermodynamics., (© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

When electrolyte solutions are confined in micro- and nanochannels their conductivity is significantly different from those in a bulk phase. Here we revisit the theory of this phenomenon by focusing attention on the reduction in the ion mobility with the concentration of salt and a consequent impact to the conductivity of a monovalent solution, from bulk to confined in a narrow slit. We first give a systematic treatment of electrophoresis of ions and obtain equations for their zeta potentials and mobilities. The latter are then used to obtain a simple expression for a bulk conductivity, which is valid in a concentration range up to a few molars and more accurate than prior analytic theories. By extending the formalism to the electrolyte solution in the charged channel the equations describing the conductivity in different modes are presented. They can be regarded as a generalization of prior work on the channel conductivity to a more realistic case of a nonzero reduction of the electrophoretic mobility of ions with salt concentration. Our analysis provides a framework for interpreting measurements on the conductivity of electrolyte solutions in the bulk and in narrow channels., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

We present an analysis of high-resolution quasi-elastic neutron scattering spectra of phosphoglycerate kinase which elucidates the influence of the enzymatic activity on the dynamics of the protein. We show that in the active state the inter-domain motions are amplified and the intra-domain asymptotic power-law relaxation ∝t-α is accelerated, with a reduced coefficient α. Employing an energy landscape picture of protein dynamics, this observation can be translated into a widening of the distribution of energy barriers separating conformational substates of the protein., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

In hyperpolarization experiments using parahydrogen, a partially negative line (PNL) of o-H2 is occasionally spotted in the nuclear magnetic resonance (NMR) spectra. This is a manifestation of the two-spin order (TSO) of the proton spins, appearing transiently in o-H2 molecules freshly derived from p-H2. For the TSO to be observable, the o-H2 NMR signal must be split into a doublet. In the literature, the splitting is believed to originate from a slow exchange of the dissolved dihydrogen with the dihydride moiety bound to a catalyst present in the reaction mixture. Because this hypothesis may be debatable, in this work a different splitting mechanism is proposed. It employs a residual dipolar coupling (RDC) between the hydrogen protons, originating from a partial orientation of the H2 molecules by the external magnetic field. The orientation effect is due to the anisotropic magnetic polarizability of H2. In a magnetic field of 11.74 T at room temperature, the currently predicted value of the RDC is -0.0024 Hz. Even such small RDC values are sufficient for the PNL effect to be clearly visible in NMR spectra for physically reasonable levels of the TSO in the o-H2 molecules. For RDC values much smaller than the natural linewidth of o-H2, the theoretical frequency distance between the minimum and maximum of PNL proves to be practically independent of the RDC and is of the order of the linewidth. The calculated amplitudes of the PNLs are proportional to the RDC values used in the calculations., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

We suggest a consistent framework for the embedding of reduced-space correlated vibrational wave functions in a potential of the remaining modes and generalize this concept to arbitrary many subspaces. We present an implementation of this framework for vibrational coupled-cluster theory and response treatments. For C=O stretches of small molecules, we show that the embedded treatment accelerates convergence for enlarging subsets. For the water dimer and trimer as well as a water wire in bacteriorhodopsin, we investigate different partitioning schemes for the embedding approach: In the local partitioning of the vibrations, the modes dominated by motions in the same spatial region are correlated, whereas in the energy-based partitioning, modes of similar fundamental frequencies are correlated. In most cases, we obtain better agreement with superset reference results for the local partitioning than for energy-based partitioning. This work represents an important step toward multi-level methodologies in vibrational-structure theory required for its application to sizable (bio-)molecular systems., (© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

The work W to form a nucleus (also known as the critical nucleus) is a key quantity in the description of nucleation phenomena because of its exponentially strong effect on the nucleation rate. The present study provides a general approximate expression for W, which comprises a hierarchy of approximations to the dependence of W on the experimentally controlled overpressure Δp of a nucleating multicomponent phase. This general expression is used to derive explicit formulas for the lowest-order members of the W(Δp) hierarchy as well as for the respective lowest-order approximations to the Δp dependences of the nucleus surface tension, the nucleus radius, the Gibbs-Tolman length, and the stationary nucleation rate. The second-order and the third-order approximations to the W(Δp) dependence are confronted with available W(Δp) data, and the latter is found to agree very well with the data. The results obtained are applicable to homogeneous single-component or multicomponent nucleation from the binodal to the spinodal of the old phase, i.e., in the entire range of the old-phase metastability., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

We investigate the role of Compton ionization in ultrafast non-resonant x-ray scattering using a molecular model system, which includes the ionization continuum via an orthonormalized plane wave ansatz. Elastic and inelastic components of the scattering signal, as well as coherent-mixed scattering that arises from electron dynamics, are calculated. By virtue of a near-quantitative distinction between scattering related to electronic transitions into bound and continuum states, we demonstrate how Compton ionization contributes to the coherent-mixed component. Analogous to inelastic scattering, the contribution to the coherent-mixed signal is significant and particularly manifests at intermediate and high-momentum transfers. Strikingly, for molecules with inversion symmetry, the exclusion of bound or continuum transitions may lead to the prediction of spurious coherent-mixed signals. We conclude that qualitative and quantitative accuracies of predicted scattering signals on detectors without energy resolution require that elements of the two-electron density operator are used. This approach inherently accounts for all accessible electronic transitions, including ionization., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

A two-component contour deformation (CD) based GW method that employs frequency sampling to drastically reduce the computational effort when assessing quasiparticle states far away from the Fermi level is outlined. Compared to the canonical CD-GW method, computational scaling is reduced by an order of magnitude without sacrificing accuracy. This allows for an efficient calculation of core ionization energies. The improved computational efficiency is used to provide benchmarks for core ionized states, comparing the performance of 15 density functional approximations as Kohn-Sham starting points for GW calculations on a set of 65 core ionization energies of 32 small molecules. Contrary to valence states, GW calculations on core states prefer functionals with only a moderate amount of Hartree-Fock exchange. Moreover, modern ab initio local hybrid functionals are also shown to provide excellent generalized Kohn-Sham references for core GW calculations. Furthermore, the core-valence separated Bethe-Salpeter equation (CVS-BSE) is outlined. CVS-BSE is a convenient tool to probe core excited states. The latter is tested on a set of 40 core excitations of eight small inorganic molecules. Results from the CVS-BSE method for excitation energies and the corresponding absorption cross sections are found to be in excellent agreement with those of reference damped response BSE calculations., (© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

Excited states are the key species in photocatalysis, while the critical parameters that govern their applications are (i) excitation energy, (ii) accessibility, and (iii) lifetime. However, in molecular transition metal-based photosensitizers, there is a design tension between the creation of long-lived excited (triplet), e.g., metal-to-ligand charge transfer (3MLCT) states and the population of such states. Long-lived triplet states have low spin-orbit coupling (SOC) and hence their population is low. Thus, a long-lived triplet state can be populated but inefficiently. If the SOC is increased, the triplet state population efficiency is improved-coming at the cost of decreasing the lifetime. A promising strategy to isolate the triplet excited state away from the metal after intersystem crossing (ISC) involves the combination of transition metal complex and an organic donor/acceptor group. Here, we elucidate the excited state branching processes in a series of Ru(II)-terpyridyl push-pull triads by quantum chemical simulations. Scalar-relativistic time-dependent density theory simulations reveal that efficient ISC takes place along 1/3MLCT gateway states. Subsequently, competitive electron transfer (ET) pathways involving the organic chromophore, i.e., 10-methylphenothiazinyl and the terpyridyl ligands are available. The kinetics of the underlying ET processes were investigated within the semiclassical Marcus picture and along efficient internal reaction coordinates that connect the respective photoredox intermediates. The key parameter that governs the population transfer away from the metal toward the organic chromophore either by means of ligand-to-ligand (3LLCT; weakly coupled) or intra-ligand charge transfer (3ILCT; strongly coupled) states was determined to be the magnitude of the involved electronic coupling., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

Tensor algebra operations such as contractions in computational chemistry consume a significant fraction of the computing time on large-scale computing platforms. The widespread use of tensor contractions between large multi-dimensional tensors in describing electronic structure theory has motivated the development of multiple tensor algebra frameworks targeting heterogeneous computing platforms. In this paper, we present Tensor Algebra for Many-body Methods (TAMM), a framework for productive and performance-portable development of scalable computational chemistry methods. TAMM decouples the specification of the computation from the execution of these operations on available high-performance computing systems. With this design choice, the scientific application developers (domain scientists) can focus on the algorithmic requirements using the tensor algebra interface provided by TAMM, whereas high-performance computing developers can direct their attention to various optimizations on the underlying constructs, such as efficient data distribution, optimized scheduling algorithms, and efficient use of intra-node resources (e.g., graphics processing units). The modular structure of TAMM allows it to support different hardware architectures and incorporate new algorithmic advances. We describe the TAMM framework and our approach to the sustainable development of scalable ground- and excited-state electronic structure methods. We present case studies highlighting the ease of use, including the performance and productivity gains compared to other frameworks., (© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

A versatile multifunctional laboratory-based near ambient pressure x-ray photoelectron spectroscopy (XPS) instrument is presented. The entire device is highly customized regarding geometry, exchangeable manipulators and sample stages for liquid- and solid-state electrochemistry, cryochemistry, and heterogeneous catalysis. It therefore delivers novel and unique access to a variety of experimental approaches toward a broad choice of functional materials and their specific surface processes. The high-temperature (electro)catalysis manipulator is designed for probing solid state/gas phase interactions for heterogeneous catalysts including solid electrolyzer/fuel cell electrocatalysts at pressures up to 15 mbar and temperatures from room temperature to 1000 °C. The liquid electrochemistry manipulator is specifically designed for in situ spectroscopic investigations of polarized solid/liquid interfaces using aqueous electrolytes and the third one for experiments for ice and ice-like materials at cryogenic temperatures to approximately -190 °C. The flexible and modular combination of these setups provides the opportunity to address a broad spectrum of in situ and operando XPS experiments on a laboratory-based system, circumventing the limited accessibility of experiments at synchrotron facilities., (© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

An intermolecular potential is reported for molecular hydrogen that combines two-body interactions from ab initio data with three-body interactions. The accuracy of the two-body potential is validated by comparison with experimental second virial coefficient data. Experimental pressure-density-temperature data are used to validate the addition of three-body interactions, often yielding very accurate predictions. Classical Monte Carlo simulations that neglect quantum effects are reported for the vapor-liquid equilibria (VLE), critical properties, and the triple point. A comparison with experimental data indicates that the effect of quantum interactions is to narrow the VLE phase envelope and to lower the critical temperature. The three-body interactions have a considerable influence on the phase behavior, resulting in good agreement with the experimental density. The critical properties of the two-body + three-body potential for hydrogen provide an alternative set of input parameters to improve the accuracy of theoretical predictions at temperatures above 100 K. In the vicinity of the critical point, the coexistence densities do not obey the law of rectilinear diameters, which is a feature that has largely been overlooked in both experimental data and reference equations of state., (© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

We here review mostly experimental and some computational work devoted to nucleation in amorphous ices. In fact, there are only a handful of studies in which nucleation and growth in amorphous ices are investigated as two separate processes. In most studies, crystallization temperatures T x or crystallization rates R JG are accessed for the combined process. Our Review deals with different amorphous ices, namely, vapor-deposited amorphous solid water (ASW) encountered in many astrophysical environments; hyperquenched glassy water (HGW) produced from μm-droplets of liquid water; and low density amorphous (LDA), high density amorphous (HDA), and very high density amorphous (VHDA) ices produced via pressure-induced amorphization of ice I or from high-pressure polymorphs. We cover the pressure range of up to about 6 GPa and the temperature range of up to 270 K, where only the presence of salts allows for the observation of amorphous ices at such high temperatures. In the case of ASW, its microporosity and very high internal surface to volume ratio are the key factors determining its crystallization kinetics. For HGW, the role of interfaces between individual glassy droplets is crucial but mostly neglected in nucleation or crystallization studies. In the case of LDA, HDA, and VHDA, parallel crystallization kinetics to different ice phases is observed, where the fraction of crystallized ices is controlled by the heating rate. A key aspect here is that in different experiments, amorphous ices of different "purities" are obtained, where "purity" here means the "absence of crystalline nuclei." For this reason, "preseeded amorphous ice" and "nuclei-free amorphous ice" should be distinguished carefully, which has not been done properly in most studies. This makes a direct comparison of results obtained in different laboratories very hard, and even results obtained in the same laboratory are affected by very small changes in the preparation protocol. In terms of mechanism, the results are consistent with amorphous ices turning into an ultraviscous, deeply supercooled liquid prior to nucleation. However, especially in preseeded amorphous ices, crystallization from the preexisting nuclei takes place simultaneously. To separate the time scales of crystallization from the time scale of structure relaxation cleanly, the goal needs to be to produce amorphous ices free from crystalline ice nuclei. Such ices have only been produced in very few studies.

The increasingly popular GW method is becoming a convenient tool to determine vertical ionization energies in molecular systems. However, depending on the formalism used and the range of orbitals investigated, it may be hampered by a steep computational scaling. To alleviate this issue, correlated natural virtual orbitals (NVOs) based on second-order Møller-Plesset (MP2) and direct MP2 correlation energies are implemented, and the resulting correlated NVOs are tested on GW quasiparticle energies. Test cases include the popular GW variants G 0 W 0 and evGW 0 as well as more elaborate vertex corrections. We find that for increasingly larger molecular systems and basis sets, NVOs considerably improve efficiency. Furthermore, we test the performance of the truncated (frozen) NVO ansatz on the GW100 test set. For the latter, it is demonstrated that, using a carefully chosen truncation threshold, NVOs lead to a negligible loss in accuracy while providing speedups of one order of magnitude. Furthermore, we compare the resulting quasiparticle energies to very accurate vertical ionization energies obtained from coupled-cluster theory with singles, doubles, and noniterative triples [CCSD(T)], confirming that the loss in accuracy introduced by truncating the NVOs is negligible compared to the methodical errors in the GW approximation. It is also demonstrated that the choice of basis set impacts results far more than using a suitably truncated NVO space. Therefore, at the same computational expense, more accurate results can be obtained using NVOs. Finally, we provide improved reference CCSD(T) values for the GW100 test set, which have been obtained using the def2-QZVPP basis set.

Dispersion forces start to play role in modern micro/nanoelectromechanical devices, but the methods to measure these forces at distances close to contact (<50 nm) suffer from pull-in instability. The method of adhered cantilever proposed recently has no instability and is able to make measurements at short separations. To measure the force at the average distance between surfaces in contact, one has to know the shape of an elastic beam with one end fixed at a height of 1-10 μm and the other end adhered to the substrate. The maximum contribution to this shape from the dispersion forces is in a range of 30-100 nm, which is well measurable by the interferometric methods. This paper describes the instrument, measurements, and data processing that make possible the reconstruction of the beam shape with an accuracy of 1 nm in a height range of at least 5000 nm. Critical steps of the fabrication procedure of cantilevers that are 12 mm long, 200 μm wide, and 10 μm thick are described. The interferometer measures the shape based on the differential interference-contrast method; the scanning is realized by a stage with a step of 0.1 μm. The signal recorded from the adhered cantilever has a noise level of 0.33 nm at a maximum sensitivity in a frequency band of 20 MHz. It is concluded that the instrument and data processing algorithm can be used to measure the dispersion forces and adhesion energies between rough surfaces in unloaded contact., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

Coupled-perturbed equations for degenerate orbitals were implemented for third order density-functional tight binding, which allowed the use of Mulliken charges as reaction coordinates. The method was applied to proton-coupled electron transfer (PCET) reactions in a model system and thoroughly tested for QM and QM/MM setups (i.e., coupled quantum and molecular mechanics). The performed enhanced sampling simulations were stable, and the obtained potentials of the mean force were able to address the thermodynamic and kinetic features of the reactions by showing the expected topography and energy barriers. Hence, this method has the potential to distinguish between concerted and sequential mechanisms and could next be applied to proton-coupled electron transfer reactions in more complex systems like proteins.

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.

The structure and scaling properties of inwardly curved polymer brushes, tethered under good solvent conditions to the inner surface of spherical shells such as membranes and vesicles, are studied by extensive molecular dynamics simulations and compared with earlier scaling and self-consistent field theory predictions for different molecular weights of the polymer chains N and grafting densities σ g in the case of strong surface curvature, R -1 . We examine the variation of the critical radius R*(σ g ), separating the regimes of weak concave brushes and compressed brushes, predicted earlier by Manghi et al. [Eur. Phys. J. E 5, 519-530 (2001)], as well as various structural properties such as the radial monomer- and chain-end density profiles, orientation of bonds, and brush thickness. The impact of chain stiffness, κ, on concave brush conformations is briefly considered as well. Eventually, we present the radial profiles of the local pressure normal, P N , and tangential, P T , to the grafting surface, and the surface tension γ(σ g ), for soft and rigid brushes, and find a new scaling relationship P N (R)∝σ g , independent of the degree of chain stiffness.4 , independent of the degree of chain stiffness.

The chemical interactions of two types of graphite and two types of carbon black (CB) with acetone, toluene, and phenol were studied in order to evaluate the influence of chemical treatment on the structure and morphology of the carbon phases. The experimental treatment of carbon phases was carried out at room temperature for 1 hour. The chemical and phase composition were studied by x-ray photoelectron (XP) and Raman spectroscopies, while the morphology and structure were determined by powder x-ray diffraction, as well as transmission electron microscopy techniques. To shed light on the most probable explanation of the observed results, we performed simulations and calculations of the binding energies of acetone, toluene, and phenol with model carbon phases: a perfect graphene sheet and a defective graphene sheet containing various structural defects (vacancies as well as zigzag and armchair edges). Simulations show that all non-covalent and most covalent coupling reactions are exothermic, with acetone coupling having the higher calorimetric effect. Based on the results of the simulations and the XP spectroscopy measurements, the probable reactions taking place during the respective treatments are outlined. The conducted studies (both theoretical and experimental) show that the treatment of graphite powders and CB with acetone, toluene, or phenol can be used as a preliminary stage of their modification and/or functionalization, including their conversion into graphene-like (defective graphene, reduced graphene oxide, and/or graphene oxide) phases. For example, the treatment of SPHERON 5000 with acetone significantly facilitates their subsequent modification with laser radiation to graphene-like phases.

A multiconfigurational adiabatic connection (AC) formalism is an attractive approach to compute the dynamic correlation within the complete active space self-consistent field and density matrix renormalization group (DMRG) models. Practical realizations of AC have been based on two approximations: (i) fixing one- and two-electron reduced density matrices (1- and 2-RDMs) at the zero-coupling constant limit and (ii) extended random phase approximation (ERPA). This work investigates the effect of removing the "fixed-RDM" approximation in AC. The analysis is carried out for two electronic Hamiltonian partitionings: the group product function- and the Dyall Hamiltonians. Exact reference AC integrands are generated from the DMRG full configuration interaction solver. Two AC models are investigated, employing either exact 1- and 2-RDMs or their second-order expansions in the coupling constant in the ERPA equations. Calculations for model molecules indicate that lifting the fixed-RDM approximation is a viable way toward improving the accuracy of existing AC approximations.

C 4 F 7 N is a promising candidate for the replacement of sulfur hexafluoride as an insulating medium, and it is important to understand the chemical changes initiated in the molecule by collision with free electrons, specifically the formation of neutral fragments. The first step of neutral fragmentation is electronic excitation, yet neither the absorption spectrum in the vacuum ultraviolet (VUV) region nor the electron energy loss spectrum have previously been reported. Here, we experimentally probed the excited states by VUV photoabsorption spectroscopy and electron energy loss spectroscopy (EELS). We found that the distribution of states populated upon electron impact with low-energy electrons is significantly different from that following photoabsorption. This difference was confirmed and interpreted with ab initio modeling of both VUV and EELS spectra. We propose here a new computational protocol for the simulation of EELS spectra combining the Born approximation with approximate forms of correlated wave functions, which allows us to calculate the (usually very expensive) scattering cross sections at a cost similar to the calculation of oscillator strengths. Finally, we perform semi-classical non-adiabatic dynamics simulations to investigate the possible neutral fragments of the molecule formed through electron-induced neutral dissociation. We show that the product distribution is highly non-statistical.

It is well established that an isolated benzene radical anion is not electronically stable. In the present study, we experimentally show that electron attachment to benzene clusters leads to weak albeit unequivocal occurrence of a C 6 H 6 moiety. We propose here-based on electronic structure calculation-that this moiety actually corresponds to linear structures formed by the opening of the benzene ring via electron attachment. The cluster environment is essential in this process since it quenches the internal energy released upon ring opening, which in the gas phase leads to further dissociation of this anion.- moiety. We propose here-based on electronic structure calculation-that this moiety actually corresponds to linear structures formed by the opening of the benzene ring via electron attachment. The cluster environment is essential in this process since it quenches the internal energy released upon ring opening, which in the gas phase leads to further dissociation of this anion.