Development of an intermolecular potential function for interactions in formamide clusters based on ab initio calculations

Construction of a Gaussian Process Regression Model of Formamide for Use in Molecular Simulations

FFLUX, a novel force field based on quantum chemical topology, can perform molecular dynamics simulations with flexible multipole moments that change with geometry. This is enabled by Gaussian process regression machine learning models, which accurately predict atomic energies and multipole moments up to the hexadecapole. We have constructed a model of the formamide monomer at the B3LYP/aug-cc-pVTZ level of theory capable of sub-kJ mol^-1 accuracy, with the maximum prediction error for the molecule being 0.8 kJ mol^-1. This model was used in FFLUX simulations along with Lennard-Jones parameters to successfully optimize the geometry of formamide dimers with errors smaller than 0.1 Å compared to those obtained with D3-corrected B3LYP/aug-cc-pVTZ. Comparisons were also made to a force field constructed with static multipole moments and Lennard-Jones parameters. FFLUX recovers the expected energy ranking of dimers compared to the literature, and changes in C═O and C-N bond lengths associated with hydrogen bonding were found to be consistent with density functional theory.

Electrostatic penetration effects stand at the heart of aromatic π interactions

The nature of the interaction in benzene-containing dimers has been analysed by means of Symmetry Adapted Perturbation Theory (SAPT). The total interaction energy and the preference for the dimers to adopt slipped structures are, apparently, consequence of the balance between repulsion and dispersion. However, our results indicate that this only holds when trends are analysed using fixed intermolecular distances. Employing the most favourable separations between rings it turns out that the changes on the total interaction energy are mostly controlled by electrostatics, while repulsion and dispersion cancel each other to a great extent. Most of the electrostatic contribution is accounted for by electrostatic penetration, so a description based on multipoles should not be employed to rationalise the interaction in benzene-containing dimers. The changes on the interaction energy in benzene-containing dimers are steered by electrostatic penetration which, though often overlooked, plays an essential role for the description of aromatic π interactions.

Computational study of interactions and nuclear magnetic shielding constants in linear chains of formamide clusters

We investigated the energetic, structural, dielectric, and nuclear magnetic shielding properties of linear n-formamide clusters, with n up to 6, to quantitatively characterize cooperative effects in model biological systems. The geometries of the complexes were optimized at the MP2 and DFT/B3LYP levels by using the pc-2 and pc-3 basis sets, while the nuclear magnetic shielding constants were calculated by employing pcS-n type basis sets, which have been optimized specifically for density functional calculations of these properties. The interaction energies show the cooperative effect, which favors the successive addition of monomers. In addition, by analyzing structural changes in the intermolecular C=O, C-N and hydrogen O⋯H bonds, as well as in the average dipole moments as cluster size increases, we found that the cooperative interaction far exceeds that expected for electrostatic interactions. Such non-pairwise-additive effects are also reflected in the changes of the nuclear magnetic shielding constants. In particular, the negativity of O shielding decreases around 23% from the monomer to the 6-formamide chain. It is possible to note the decrease in the shielding of H and in the deshielding of O as a result of their hydrogen bonding. However, the results obtained show that these variations in the extremes of formamide chains tend to zero, and the respective shielding values tend to stabilize as the number of monomers increases in the chain. Also, the cooperative effect increases in the middle of the chains, by decreasing the shielding for all atoms except that of O, which decreases its deshielding. These results could serve to guide improvements in current conventional models for simulating hydrogen bonded systems.

Structural investigations of N-methylformamide-water mixtures at various concentrations

Structural investigations of N-methylformamide-water mixtures (NMF-water) are performed at room temperature and atmospheric pressure for two water molar fractions x w = 0.66 and x w = 0.75 . This paper extends our recent study on the equimolar system. H-bond networks are preferentially formed between NMF and water molecules. Among a large variety of DFT optimized models, X-ray scattering data shows that the local order of each mixture is better described by a tetramer where one NMF molecule is connected to three water molecules. No self-association is observed in the considered systems. The effect of hydration is compared to the temperature and pressure effects in some hydrogen-bonded liquids.

Kinetically controlled formation of formamide trimer from first principles

The formation of formamide trimers was simulated using Car-Parrinello molecular dynamics. A variety of different initial setups were compared to study the effects of spatial arrangement, concentration, and temperature on the trimer product distribution. A total of nine different trimer species were obtained by simulation. Although a triangular initial arrangement of the three monomers is found to favour a less energetically stable chain-like product at high concentration, the more compact global minimum structure is expected to be the most abundant species overall in experiment. This is because there is evidence of a low activation barrier for conversion of the chain-like trimer to the lowest-energy structure. For one, this process is observed upon increasing the length of the trajectories. Furthermore, a slight rise in temperature drastically reduces the number of chain-like trimers. With regard to the intermolecular forces driving the aggregation dynamics, dispersion corrections to the underlying density functional theory description have a strong effect on the product distribution, again favouring the global minimum species. Certain local minimum structures are significantly destabilised relative to the global minimum by dispersion correction while the relative energies of the majority of species are practically unchanged.

Formamide Dimers: A Computational and Matrix Isolation Study

The dimerization of formamide (FMA) has been investigated by matrix isolation spectroscopy, static ab initio calculations, and ab initio molecular dynamics (AIMD) simulations. Comparison of the experimental matrix IR spectra with the ab initio calculations reveals that two types of dimers A and C are predominantly formed, with two and one strong NH...O hydrogen bonds, respectively. This is in accordance with previously published experiments. In addition, there is also experimental evidence for the formation of the thermally labile dimer B after deposition of high concentrations of FMA in solid xenon. The AIMD simulations of the aggregation process show that in all cases dimer C is initially formed, but rearrangement to the more stable doubly hydrogen-bonded structures A or B occurs for a fraction of collisions on the sub-picosecond time scale.

Electron localization in negatively charged formamide clusters studied by photodetachment spectroscopy

Size-dependent features of the electron localization in negatively charged formamide clusters (FAn-, n = 5-21) have been studied by photodetachment spectroscopy. In the photoelectron spectra for all the sizes studied, two types of bands due to different isomers of anions were found. The low binding energy band peaking around 1 eV is assigned to the solvated electron state by relative photodetachment cross-section measurements in the near-infrared region. It is suggested that nascent electron trapping is dominated by formation of the solvated electron. The higher energy band originates from the covalent anion state generated after a significant relaxation process, which exhibits a rapid increase of electron binding energy as a function of the cluster size. A unique behavior showing a remarkable band intensity of the higher energy band was found only for n = 9.

The Nature of Interactions in the Ionic Crystal of 3-Pentenenitrile, 2-Nitro-5-oxo, Ion(−1), Sodium

The hybrid variation -- perturbation many-body interaction energy decomposition scheme has been applied to analyze the physical nature of interactions in the ionic 3-pentenenitrile, 2-nitro-5-oxo, ion(-1), sodium crystal, which can be regarded as a model for a large group of aromatic quaternary nitrogen salts. In the crystal structure the sodium ions and water molecules of adjacent unit cells form a positively charged "inorganic layer" with the sodium ions clustered together along the ab faces with the organic (negative) part in between. This puzzling crystal packing is due to a strong favorable interaction between the water molecule and the sodium ions and a substantial charge transfer from the carbanions that balances out the destabilizing sodium-sodium ion repulsion. Although the majority of cohesion energy of the crystal structure comes from the electrostatic interactions of ions, the resulting net stabilization also depends heavily on the nonadditive delocalization components, due to a counterbalance between the two-body delocalization and exchange effects. The estimated nonadditivity of interactions varies between 12% and 22%.

Inter- and intramolecular potential for the N-formylglycinamide-water system. A comparison between theoretical modeling and empirical force fields

An intramolecular NEMO potential is presented for the N-formylglycinamide molecule together with an intermolecular potential for the N-formylglycinamide-water system. The intramolecular N-formylglycinamide potential can be used as a building block for the backbone of polypeptides and proteins. Two intramolecular minima have been obtained. One, denoted as C5, is stabilized by a hydrogen bonded five member ring, and the other, denoted as C7, corresponds to a seven membered ring. The interaction between one water molecule and the N-formylglycinamide system is also studied and compared with Hartree-Fock SCF calculations and with the results obtained for some of the more commonly used force fields. The agreement between the NEMO and SCF energies for the complexes is in general superior to that of the other force fields. In the C7 region the surfaces obtained from the intramolecular part of the commonly used force fields are too flat compared to the NEMO potential and the ab initio calculations. We further analyze the possibility of using a charge distribution obtained from one conformation to describe the charge distribution of other conformations. We have found that the use of polarizabilities and generic dipoles can model most of the changes in charge density due to the different geometry of the new conformations, but that one can expect additional errors in the interaction energies that are of the order of 1 kcal/mol.