Neutron Scattering of Aromatic and Aliphatic Liquids

Ab initio electronic absorption spectra of para-nitroaniline in different solvents: Intramolecular charge transfer effects

Intramolecular charge transfer (ICT) effects of para-nitroaniline (pNA) in eight solvents (cyclohexane, toluene, acetic acid, dichloroethane, acetone, acetonitrile, dimethylsulfoxide, and water) are investigated extensively. The second-order algebraic diagrammatic construction, ADC(2), ab initio wave function is employed with the COSMO implicit and discrete multiscale solvation methods. We found a decreasing amine group torsion angle with increased solvent polarity and a linear correlation between the polarity and ADC(2) transition energies. The first absorption band involves π → π* transitions with ICT from the amine and the benzene ring to the nitro group, increased by 4%-11% for different solvation models of water compared to the vacuum. A second band of pNA is characterized for the first time. This band is primarily a local excitation on the nitro group, including some ICT from the amine group to the benzene ring that decreases with the solvent polarity. For cyclohexane, the COSMO implicit solvent model shows the best agreement with the experiment, while the explicit model has the best agreement for water.

Neutron Scattering Studies of Heterogeneous Catalysis

Understanding the structural dynamics/evolution of catalysts and the related surface chemistry is essential for establishing structure-catalysis relationships, where spectroscopic and scattering tools play a crucial role. Among many such tools, neutron scattering, though less-known, has a unique power for investigating catalytic phenomena. Since neutrons interact with the nuclei of matter, the neutron-nucleon interaction provides unique information on light elements (mainly hydrogen), neighboring elements, and isotopes, which are complementary to X-ray and photon-based techniques. Neutron vibrational spectroscopy has been the most utilized neutron scattering approach for heterogeneous catalysis research by providing chemical information on surface/bulk species (mostly H-containing) and reaction chemistry. Neutron diffraction and quasielastic neutron scattering can also supply important information on catalyst structures and dynamics of surface species. Other neutron approaches, such as small angle neutron scattering and neutron imaging, have been much less used but still give distinctive catalytic information. This review provides a comprehensive overview of recent advances in neutron scattering investigations of heterogeneous catalysis, focusing on surface adsorbates, reaction mechanisms, and catalyst structural changes revealed by neutron spectroscopy, diffraction, quasielastic neutron scattering, and other neutron techniques. Perspectives are also provided on the challenges and future opportunities in neutron scattering studies of heterogeneous catalysis.

Solvation Effects on Polarizability of Aromatic Fluids

Elucidating solvation effects on polarizability in condensed phases is important for the description of the optical and dielectric behavior of high-refractive-index molecular materials. We study these effects via the polarizability model combining electronic, solvation, and vibrational contributions. The method is applied to well-characterized highly polarizable liquid precursors: benzene, naphthalene, and phenanthrene. We find that the solvation and vibrational terms are of opposite signs and cancel almost exactly for benzene, but for naphthalene and phenanthrene, a 2.5 and 5.0% decrease relative to the equilibrium electronic polarizability of the respective monomer, α₁^e, is predicted, respectively. The increase in electronic polarizability affects interaction polarizability of all contacts, which is the main reason for the increasing importance of solvation contribution. The calculated refractive indices agree very well with experiment for all three systems.

A variational formulation of the Harris functional as a correction to approximate Kohn-Sham density functional theory

Accurate descriptions of intermolecular interactions are of great importance in simulations of molecular liquids. We present an electronic structure method that combines the accuracy of the Harris functional approach with the computational efficiency of approximately linear-scaling density functional theory (DFT). The non-variational nature of the Harris functional has been addressed by constructing a Lagrangian energy functional, which restores the variational condition by imposing stationarity with respect to the reference density. The associated linear response equations may be treated with linear-scaling efficiency in an atomic orbital based scheme. Key ingredients to describe the structural and dynamical properties of molecular systems are the forces acting on the atoms and the stress tensor. These first-order derivatives of the Harris Lagrangian have been derived and implemented in consistence with the energy correction. The proposed method allows for simulations with accuracies close to the Kohn-Sham DFT reference. Embedded in the CP2K program package, the method is designed to enable ab initio molecular dynamics simulations of molecular solutions for system sizes of several thousand atoms. Available subsystem DFT methods may be used to provide the reference density required for the energy correction at near linear-scaling efficiency. As an example of production applications, we applied the method to molecular dynamics simulations of the binary mixtures cyclohexane-methanol and toluene-methanol, performed within the isobaric-isothermal ensemble, to investigate the hydrogen bonding network in these non-ideal mixtures.

The structure of protic ionic liquids based on sulfuric acid, doped with excess of sulfuric acid or with water

Neutron scattering with isotopic substitution was used to study the structure of concentrated sulfuric acid, and two protic ionic liquids (PILs): a Brønsted-acidic PIL, synthesised using pyridine and excess of sulfuric acid, [Hpy][HSO₄]·H₂SO₄, and a hydrated PIL, in which an equimolar mixture of sulfuric acid and pyridine has been doped with water, [Hpy][HSO₄]·2H₂O. Brønsted acidic PILs are excellent solvents/catalysts for esterifications, driving reaction to completion by phase-separating water and ester products. Water-doped PILs are efficient solvents/antisolvents in biomass fractionation. This study was carried out to provide an insight into the relationship between the performance of PILs in the two respective processes and their liquid structure. It was found that a persistent sulfate/sulfuric acid/water network structure was retained through the transition from sulfuric acid to PILs, even in the presence of 2 moles (∼17 wt%) of water. Hydrogen sulfate PILs have the propensity to incorporate water into hydrogen-bonded anionic chains, with strong and directional hydrogen bonds, which essentially form a new water-in-salt solvent system, with its own distinct structure and physico-chemical properties. It is the properties of this hydrated PIL that can be credited both for the good performance in esterification and beneficial solvent/antisolvent behaviour in biomass fractionation.

Photoelectron Spectroscopy of Benzene in the Liquid Phase and Dissolved in Liquid Ammonia

We report valence band photoelectron spectroscopy measurements of gas-phase and liquid-phase benzene as well as those of benzene dissolved in liquid ammonia, complemented by electronic structure calculations. The origins of the sizable gas-to-liquid-phase shifts in electron binding energies deduced from the benzene valence band spectral features are quantitatively characterized in terms of the Born-Haber solvation model. This model also allows to rationalize the observation of almost identical shifts in liquid ammonia and benzene despite the fact that the former solvent is polar while the latter is not. For neutral solutes like benzene, it is the electronic polarization response determined by the high frequency dielectric constant of the solvent, which is practically the same in the two liquids, that primarily determines the observed gas-to-liquid shifts.

Bulk and Confined Benzene-Cyclohexane Mixtures Studied by an Integrated Total Neutron Scattering and NMR Method

Herein mixtures of cyclohexane and benzene have been investigated in both the bulk liquid phase and when confined in MCM-41 mesopores. The bulk mixtures have been studied using total neutron scattering (TNS), and the confined mixtures have been studied by a new flow-utilising, integrated TNS and NMR system (Flow NeuNMR), all systems have been analysed using empirical potential structure refinement (EPSR). The Flow NeuNMR setup provided precise time-resolved chemical sample composition through NMR, overcoming the difficulties of ensuring compositional consistency for computational simulation of data ordinarily found in TNS experiments of changing chemical composition—such as chemical reactions. Unique to the liquid mixtures, perpendicularly oriented benzene molecules have been found at short distances from the cyclohexane rings in the regions perpendicular to the carbon–carbon bonds. Upon confinement of the hydrocarbon mixtures, a stronger parallel orientational preference of unlike molecular dimers, at short distances, has been found. At longer first coordination shell distances, the like benzene molecular spatial organisation within the mixture has also found to be altered upon confinement.

In situ total scattering experiments of nucleation and crystallisation of tantalum-based oxides: from highly dilute solutions via cluster formation to nanoparticles

The exact formation mechanism of tantalum oxides (and in general, metal/mixed metal oxides) from alkoxide precursors is still not fully understood, particularly when forming cluster-like or amorphous materials. The structural evolution of Ta-based oxides was studied in detail using X-ray total scattering experiments along with subsequent pair distribution function (PDF) analyses. Starting from a tantalum alkoxide precursor (Ta₂(OEt)₁₀), the formation of hydrolysed Ta_xO_yH_z clusters in highly diluted aqueous solution was analysed. From the PDF data, the connectivity and arrangement of Ta_xO_y octahedra in the cluster could be deduced as well as the approximate size of the clusters (<1 nm). Construction of cluster models allowed for identification of common structural motifs in the Ta_xO_yH_z clusters, ruling out the formation of chain- or ring-like clusters. More likely, bulky clusters with a high number of corner-sharing octahedra are formed. After separation of the amorphous solid from the liquid, temperature-induced crystallisation processes were monitored via in situ total scattering experiments. Between room temperature and 600 °C, only small rearrangements of the amorphous structure are observed. At about 610 °C, amorphous Ta_xO_yH_z transforms directly into crystalline orthorhombic L-Ta₂O₅ without formation of any crystalline intermediate structures.

Nanoarchitecture through Strained Molecules: Cubane-Derived Scaffolds and the Smallest Carbon Nanothreads

Relative to the rich library of small-molecule organics, few examples of ordered extended (i.e., nonmolecular) hydrocarbon networks are known. In particular, sp³ bonded, diamond-like materials represent appealing targets because of their desirable mechanical, thermal, and optical properties. While many covalent organic frameworks (COFs)-extended, covalently bonded, and porous structures-have been realized through molecular architecture with exceptional control, the design and synthesis of dense, covalent extended solids has been a longstanding challenge. Here we report the preparation of a sp³-bonded, low-dimensional hydrocarbon synthesized via high-pressure, solid-state diradical polymerization of cubane (C₈H₈), which is a saturated, but immensely strained, cage-like molecule. Experimental measurements show that the obtained product is crystalline with three-dimensional order that appears to largely preserve the basic structural topology of the cubane molecular precursor and exhibits high hardness (comparable to fused quartz) and thermal stability up to 300 °C. Among the plausible theoretical candidate structures, one-dimensional carbon scaffolds comprising six- and four-membered rings that pack within a pseudosquare lattice provide the best agreement with experimental data. These diamond-like molecular rods with extraordinarily small thickness are among the smallest members in the carbon nanothread family, and calculations indicate one of the stiffest one-dimensional systems known. These results present opportunities for the synthesis of purely sp³-bonded extended solids formed through the strain release of saturated molecules, as opposed to only unsaturated precursors.

Probing the dynamics and structure of confined benzene in MCM-41 based catalysts

Combined MD simulations and QENS experiments on benzene in MCM-41 provide insight into the dynamics and structure of benzene

Frustrated Lewis pairs in ionic liquids and molecular solvents - a neutron scattering and NMR study of encounter complexes

The presence of the weakly-associated encounter complex in the model frustrated Lewis pair solution (FLP): tris(tert-butyl)phosphine (P(^tBu)₃) and tris(pentafluorophenyl)borane (BCF) in benzene, was confirmed via PB correlation analysis from neutron scattering data. On average, ca. 5% of dissolved FLP components were in the associated state. NMR spectra of the FLP in benzene gave no evidence of such association, in agreement with earlier reports and the transient nature of the encounter complex. In contrast, the corresponding FLP solution in the ionic liquid, 1-decyl-3-methylimidazolium bistriflamide, [C₁₀mim][NTf₂], generated NMR signals that can be attributed to formation of encounter complexes involving over 20% of the dissolved species. The low diffusivity characteristics of ionic liquids is suggested to enhance high populations of encounter complex. The FLP in the ionic liquid solution retained its ability to split hydrogen.

Impact of Molecular Structure on Properties of n-Hexadecane and Alkylbenzene Binary Mixtures

Because of the complexity of petroleum-based fuels, researchers typically use simplified mixtures, known as surrogates, to study combustion behavior and to attempt to identify how physical properties are related to combustion. The process of determining the surrogate composition to yield a desired set of thermophysical properties can be a complicated and time-consuming task. As a result, the use of computer simulations to narrow the number of possible surrogate compositions is beginning to be explored. Herein, molecular dynamics (MD) simulations are used to model binary mixtures of n-hexadecane with either benzene, toluene, n-ethylbenzene, n-propylbenzene, or n-butylbenzene. Calculated densities are in quantitative agreement with experimental values. With the exception of the mixtures containing benzene, simulated excess molar volumes are also in very good agreement with measured values. Isentropic bulk moduli are in qualitative agreement with experiment, and reproduce interesting trends observed in the experimental data. Specifically, minima in the bulk moduli at intermediate compositions of several of the alkylbenzenes are correctly reproduced. In addition, the structures of the fluids are also examined. For mixtures of n-hexadecane with alkylbenzenes with longer chains, the orientation of the aromatic rings is not substantially impacted by composition. In contrast, increasing n-hexadecane content increases the ratio of parallel to perpendicular arrangements of benzene and toluene molecules. In those mixtures, this change in orientation of the aromatic rings could be responsible for the minima observed in the bulk moduli data. These results show that MD simulations can assist in development of fuel surrogates, both by predicting thermophysical properties and by providing insight into how molecular structure and composition affect those properties.

Confinement Effects on the Benzene Orientational Structure

Liquids under confinement exhibit different properties compared with their corresponding bulk phases, for example, miscibility, phase transitions, and diffusion. The underlying cause is the local ordering of molecules, which is usually only studied using pure simulation methods. Herein, we derive experimentally the structure of benzene confined in MCM-41 using total neutron scattering measurements. The study reveals a layering of molecules across a pore, and four concentric cylindrical shells can be distinguished for a pore with the radius of 18 Å. The nanoscale confinement of the liquid has a major effect on the spatial and orientational correlations observed between the molecules, when compared with the structure of the bulk liquid. These differences are most marked for molecules in parallel configurations, and this suggests differences in chemical reactivity between the confined and bulk liquids.

Computational study of the DPAP molecular rotor in various environments: from force field development to molecular dynamics simulations and spectroscopic calculations

Fluorescent molecular rotors (FMRs) belong to an important class of environment-sensitive dyes capable of acting as nanoprobes in the measurement of viscosity and polarity of their micro-environment. FMRs have found widespread applications in various research fields, ranging from analytical to biochemical sciences, for example in intracellular imaging studies or in volatile organic compound detection. Here, a computational investigation of a recently proposed FMR, namely 4-(diphenylamino)phthalonitrile (DPAP), in various chemical environments is presented. A purposely developed molecular mechanics force field is proposed and then applied to simulate the rotor in a high- and low-polar solvent (i.e., acetonitrile, tetrahydrofuran, o-xylene and cyclohexane), a polymer matrix and a lipid membrane. Subtle effects of the molecular interactions with the embedding medium, the structural fluctuations of the rotor and its rotational dynamics are analyzed in some detail. The results correlate with a previous work, thus supporting the reliability of the model, and provide further insights into the environment-specific properties of the dye. In particular, it is shown how molecular diffusion and rotational correlation times of the FMR are affected by the surrounding medium and how the molecular orientation of the dye becomes anisotropic once immersed in the lipid bilayer. Moreover, a qualitative correlation between the FMR rotational dynamics and the fluorescence lifetime is detected, a result in line with the observed viscosity dependence of its emission. Finally, optical absorption spectra are computed and successfully compared with their experimental counterparts.

Dihydrogen vs. hydrogen bonding in the solvation of ammonia borane by tetrahydrofuran and liquid ammonia

The solvation structures of two systems rich in hydrogen and dihydrogen bonding interactions have been studied in detail experimentally through neutron diffraction with hydrogen/deuterium isotopic substitution.

An integrated total neutron scattering – NMR approach for the study of heterogeneous catalysis

By combining total neutron scattering with nuclear magnetic resonance (NeuNMR) in a single experimental apparatus, new insights into the kinetics and mechanisms of heterogeneous catalytic reactions occurring in situ, within the catalyst pore space, are possible.

The structures of liquid pyridine and naphthalene: the effects of heteroatoms and core size on aromatic interactions

Spatial and orientational structures of liquid naphthalene and pyridine revealed using neutron scattering combined with empirical potential structure refinement.

The Structure of Ethylbenzene, Styrene and Phenylacetylene Determined by Total Neutron Scattering

Organic solvents such as phenylacetylene, styrene and ethylbenzene are widely used in industrial processes, especially in the production of rubber or thermoplastics. Despite their important applications detailed knowledge about their structure is limited. In this paper the structures of these three aromatic solvents were investigated using neutron diffraction. The results show that many of their structural characteristics are similar, although the structure of phenylacetylene is more ordered and has a smaller solvation sphere than either ethylbenzene or styrene. Two regions within the first coordination sphere, in which the surrounding molecules show different preferable orientations with respect to the central molecule, were found for each liquid. Additionally, the localisation of the aliphatic chains reveals that they tend to favour closer interactions with each other than to the aromatic rings of the adjacent molecules.

Neutron Scattering of Aromatic and Aliphatic Liquids

Organic solvents, such as cyclohexane, cyclohexene, methylcyclohexane, benzene and toluene, are widely used as both reagents and solvents in industrial processes. Despite the ubiquity of these liquids, the local structures that govern the chemical properties have not been studied extensively. Herein, we report neutron diffraction measurements on liquid cyclohexane, cyclohexene, methylcyclohexane, benzene and toluene at 298 K to obtain a detailed description of the local structure in these compounds. The radial distribution functions of the centres of the molecules, as well as the partial distribution functions for the double bond for cyclohexene and methyl group for methylcyclohexane and toluene have been calculated. Additionally, probability density functions and angular radial distribution functions were extracted to provide a full description of the local structure within the chosen liquids. Structural motifs are discussed and compared for all liquids, referring specifically to the functional group and aromaticity present in the different liquids.

Determination of toluene hydrogenation kinetics with neutron diffraction

Total neutron scattering has been used to follow the hydrogenation of toluene-d₈ to methylcyclohexane-d₁₄ over 3 wt% platinum supported on highly ordered mesoporous silica (MCM-41) at 298 K and under 150–250 mbar D₂ pressure.