A Monte Carlo study of methanol clusters (CH3OH)N, N=5–256

Unraveling environmental effects in the absorption and fluorescence spectra of p-methoxyphenylpiperazine derivatives

The p-methoxyphenylpiperazine motif can be found in many biologically active molecules, including approved drugs. It is characterized by a relatively weak fluorescence, which can be employed in different types of studies involving molecules with this motif. In this work, a thorough analysis of the absorption, excitation and emission spectra of the diphenyl(aminomethyl)phosphine and tris(aminomethyl)phosphine derivatives of p-methoxyphenylpiperazine, supported by the DFT calculations (ωB97XD/6-311++G(d,p)) with NBO and QTAIM analysis also for different model molecules (e.g. 1-(4-methoxyphenyl)-4-methylpiperazine) enabled determination of the mechanisms underlying beneath the electronic transitions and allowed to rationalize mixed solvent effects observed in electronic spectra of the studied compounds. Electronic transition from the ground state to the first excited state can be regarded as the n,π → π* transition with no solvatochromic effects, however the hydrogen bonds between the HBD solvent molecules and the nitrogen atom bound directly to the aromatic ring (N(4)) are shifting strongly the 1st absorption or excitation band maxima to the higher energies. Fluorescence band, as a result of the electron transition from the equilibrated 1st excited state to the ground state, can be described as the π*→π with positive solvatochromism. N(4) in the excited states adopts a sp2 hybridization and is no longer able to form HBs. On the other hand, increased electron density on the aromatic ring makes the emission processes vulnerable to its direct environment.

Capturing the potential energy landscape of large size molecular clusters from atomic interactions up to a 4-body system using deep learning

We propose a new method that utilizes the database of stable conformers and borrow the fragmentation concept of many-body-expansion (MBE) methods in ab initio methods to train a deep-learning machine learning (ML) model using SchNet.

Solvation Energies of the Proton in Methanol

pKa's, proton affinities, and proton dissociation free energies characterize numerous properties of drugs and the antioxidant activity of some chemical compounds. Even with a higher computational level of theory, the uncertainty in the proton solvation free energy limits the accuracy of these parameters. We investigated the thermochemistry of the solvation of the proton in methanol within the cluster-continuum model. The scheme used involves up to nine explicit methanol molecules, using the IEF-PCM and the strategy based on thermodynamic cycles. All computations were performed at B3LYP/6-31++G(dp) and M062X/6-31++G(dp) levels of theory. It comes out from our calculations that the functional M062X is better than B3LYP, on the evaluation of gas phase clustering energies of protonated methanol clusters, per methanol stabilization of neutral methanol clusters and solvation energies of the proton in methanol. The solvation free energy and enthalpy of the proton in methanol were obtained after converging the partial solvation free energy of the proton in methanol and the clustering free energy of protonated methanol clusters, as the cluster size increases. Finally, the recommended values for the solvation free energy and enthalpy of the proton in methanol are -257 and -252 kcal/mol, respectively.

A Kirkwood-Buff force field for the aromatic amino acids

In a continuation of our efforts to develop a united atom non-polarizable protein force field based upon the solution theory of Kirkwood and Buff i.e., the Kirkwood-Buff Force Field (KBFF) approach, we present KBFF models for the side chains of phenylalanine, tyrosine, tryptophan, and histidine, including both tautomers of neutral histidine and doubly-protonated histidine. The force fields were specifically designed to reproduce the thermodynamic properties of mixtures over the full composition range in an attempt to provide an improved description of intermolecular interactions. The models were developed by careful parameterization of the solution phase partial charges to reproduce the experimental Kirkwood-Buff integrals for mixtures of solutes representative of the amino acid sidechains in solution. The KBFF parameters and simulated thermodynamic and structural properties are presented for the following eleven binary mixtures: benzene + methanol, benzene + toluene, toluene + methanol, toluene + phenol, toluene + p-cresol, pyrrole + methanol, indole + methanol, pyridine + methanol, pyridine + water, histidine + water, and histidine hydrochloride + water. It is argued that the present approach and models provide a reasonable description of intermolecular interactions which ensures that the required balance between solute-solute, solute-solvent, and solvent-solvent distributions is obtained.

Chemistry of Styrene (Water)n Clusters, n = 1−5: Spectroscopy and Structure of the Neutral Clusters, Deprotonation of Styrene Dimer Cation, and Implication to the Inhibition of Cationic Polymerization

The styrene-water binary clusters SW(n), with n = 1-5 have been studied by the (one-color) resonant two-photon ionization technique using the resonance of styrene. The structures and energetics of the neutral clusters are investigated using a search technique that employs Monte Carlo procedure. The strong tendency for water molecules to form cyclic hydrogen-bonded structures is clearly observed in the SW(n) structures starting from n =3. The results indicate that the spectral shifts correlate with the interaction energies between styrene and the water subcluster (W(n)) within the SW(n) clusters. Evidence is presented that points to (1) the formation of a covalent bonded styrene radical cation dimer following the 193 nm MPI of styrene neutral clusters, (2) proton transfer from the styrene dimer cation to the water or methanol subcluster, resulting in the formation of protonated water or methanol clusters and a styrene dimer radical, and (3) extensive solvation of the styrene dimer radical within the protonated solvent molecules. The proton-transfer reactions may explain the strong inhibition effects exerted by small concentrations of water or methanol on the cationic polymerization of styrene. These results provide a molecular level view of the inhibition mechanism exerted by protic solvents on the cationic polymerization of styrene.

Gas-Phase Ion Mobilities and Structures of Benzene Cluster Cations (C6H6)n+, n = 2−6

Benzene clusters are generated by pulsed supersonic beam expansion, ionized by electron impact, mass-selected and then injected into a drift cell for ion mobility measurements in a helium buffer gas. The measured collision cross sections and theoretical calculations are used to determine the structures of the cluster cations (C(6)H(6))(n)(+) with n = 2-6. Density functional theory calculation, at an all-electron level and without any symmetry constraint, predicts that the dimer cation has two nearly degenerate ground state structures with the sandwich configuration more stable than the T-configuration by only 0.07 eV. The ion mobility experiment indicates that only one structure is observed for the mass-selected dimer cation at room temperature. The calculated cross section for the sandwich structure agrees very well (within 2.4%) with the experimental value. For the n = 3-6 clusters, the experiments suggest the presence of at least two structural isomers for each cluster. A Monte Carlo minimum-energy search technique using the 12-site OPLS potential for benzene is used to determine the structures of the lowest-energy isomers. The calculated cross sections for the two lowest-energy isomers of the n = 3-6 clusters agree well with the experimental results. The clusters' structures reveal two different growth patterns involving a sandwich dimer core or a pancake trimer stack core. The lowest-energy isomers of the n = 3-6 clusters incorporate the pancake trimer stack as the cluster's core. The trimer stack allows the charge to hop between two dimers, thus maximizing charge resonance interaction in the clusters. For larger clusters, the appearance of magic numbers at n = 14, 20, 24, 27, and 30 is consistent with the incorporation of a sandwich dimer cation within icosahedral, double icosahedral, and interpenetrating icosahedral structures. On the basis of the ion mobility results and the structural calculations, the parallel-stacked motif among charged aromatic-aromatic interactions is expected to play a major role in determining the structures of multi aromatic components. This conclusion may provide new insights for experimental and theoretical studies of molecular design and recognition involving aromatic systems.