Direct experimental assessment of the strength of orientational correlations in polar liquids

Thermodynamic calculations using reverse Monte Carlo: Simultaneously tuning multiple short-range order parameters for 2D lattice adsorption problem

Lattice simulations are an important class of problems in crystalline solids, surface science, alloys, adsorption, absorption, separation, catalysis, to name a few. We describe a fast computational method for performing lattice thermodynamic calculations that is based on the use of the reverse Monte Carlo (RMC) technique and multiple short-range order (SRO) parameters. The approach is comparable in accuracy to the Metropolis Monte Carlo (MC) method. The equilibrium configuration is determined in 5-10 Newton-Raphson iterations by solving a system of coupled nonlinear algebraic SRO growth rate equations. This makes the RMC-based method computationally more efficient than MC, given that MC typically requires sampling of millions of configurations. The technique is applied to the interacting 2D adsorption problem. Unlike grand canonical MC, RMC is found to be adept at tackling geometric frustration, as it is able to quickly and correctly provide the ordered c(2 × 2) adlayer configuration for Cl adsorbed on a Cu (100) surface.

A probabilistic microkinetic modeling framework for catalytic surface reactions

We present a probabilistic microkinetic modeling (MKM) framework that incorporates the short-ranged order (SRO) evolution for adsorbed species (adspecies) on a catalyst surface. The resulting model consists of a system of ordinary differential equations. Adsorbate-adsorbate interactions, surface diffusion, adsorption, desorption, and catalytic reaction processes are included. Assuming that the adspecies ordering/arrangement is accurately described by the SRO parameters, we employ the reverse Monte Carlo (RMC) method to extract the relevant local environment probability distributions and pass them to the MKM. The reaction kinetics is faithfully captured as accurately as the kinetic Monte Carlo (KMC) method but with a computational time requirement of few seconds on a standard desktop computer. KMC, on the other hand, can require several days for the examples discussed. The framework presented here is expected to provide the basis for wider application of the RMC-MKM approach to problems in computational catalysis, electrocatalysis, and material science.

The structure of liquid water beyond the first hydration shell

Distance dependent excess entropy calculations reveal that water is tetrahedrally structured up to long distances.

On the positional and orientational order of water and methanol around indole: a study on the microscopic origin of solubility

The increase in solubility for indole in methanol water solutions relative to pure water is a result methanol −OH–π interactions. In addition, excess entropy calculations suggest that this process is enthalpically rather than entropically driven.

Hyper-Rayleigh scattering from correlated molecules

The polarization dependence of hyper-Rayleigh scattering has been calculated for spherical domains of orientation correlated molecules. Distributions with radial or azimuthal mean polar orientation of the molecules are found that give results consistent with experimental observations, and expressions for the polarization ratios in terms of the product of correlation strength and correlated domain size are derived for these distributions. Assuming a plausible correlation strength, it is estimated that the correlated domain size in typical polar liquids is of order 100 molecular diameters.

New microscopic mechanism for secondary relaxation in glasses

The dynamics of simple molecular systems showing glassy properties has been explored by dielectric spectroscopy and nuclear quadrupole resonance (NQR) on the halogenomethanes CBr2Cl2 and CBrCl3 in their low-temperature monoclinic phases. The dielectric spectra display features which correspond to alpha- and beta-relaxation processes, commonly observed in canonical glass formers. NQR experiments, also performed in the ergodic monoclinic phase of CCl4, enable the determination of the microscopic mechanism underlying the beta dynamics in these simple model glasses: Molecules that are nonequivalent with respect to their molecular environment perform reorientational jumps at different time scales. Thus our findings reveal another mechanism that can give rise to typical beta-relaxation behavior, raising some doubt about the existence of a universal explanation of this phenomenon.

Multiple crossed isopolymorphism: two-component systems CCl4 +CBr2Cl2 and CBrCl3 + CBr2Cl2; inference of a metastable rhombohedral phase of CBr2Cl2

The phases diagrams of the two-component systems CCl4 +CBr2Cl2 and CBrCl3 + CBr2Cl2 have been determined by means of X-ray powder diffraction and thermal analysis techniques from the low-temperature ordered phase to the liquid state. The isomorphism relationship between the stable orientationally disordered (OD) face-centered cubic (FCC) phases of CBrCl3 and CBr2Cl2 and the metastable OD FCC phase (monotropic behavior with respect to the OD rhombohedral stable phase) of CCl4 has been put into evidence throughout the continuous evolution of the lattice parameters and the existence of the two-phase equilibrium [FCC + L] for the whole range of composition in both two-component systems. This equilibrium interferes, for the CCl4 +CBr2Cl2 system, with a rhombohedral (R) plus liquid ([R + L]) equilibrium giving rise to a peritectic invariant. In addition, whatever the system, [R + FCC] equilibrium also interferes with the low-temperature equilibria between the low-temperature monoclinic (C2/c) phase and the OD R and FCC phases. In regards to the low-temperature monoclinic phases, isomorphism is evidenced, and by means of Rietveld profile refinement, any ordering of the molecules by varying the fractional occupancy of the halogen sites has been detected. The thermodynamic assessment, conducted by means of the concept of crossed isopolymorphism, coherently reproduces all the involved equilibria and provides a coherent set of data for the thermodynamic properties of nonexperimentally available phase transitions of pure compound CBr2Cl2 which enables us to obtain the topological properties of its pressure-temperature phase diagram and to infer the existence of a high-pressure R phase for such a compound.

Quantitative characterization of orientational order in liquid carbon tetrachloride

A simple geometrical construct is proposed for a clear-cut classification of the relative orientation between two tetrahedral molecules in terms of six orientational classes. When applied to sort out configurations from condensed phase simulations, it leads to a quantitative characterization of orientational order: A definite percentage for each class is obtained as a function of the distance between molecular centers. The basic picture that emerges, for liquid carbon tetrachloride, is that the dominant configuration for each distance is such that the number of chlorines in between both carbons diminishes with increasing separation, with a configuration here termed edge-to-face being the dominant one at contact. Regarding the range of orientational order, remnants are still noticeable at approximately 20 A, i.e., up to the fourth solvation shell. Beyond this distance the distributions are hardly distinguishable from the analytical predictions for random orientation. The analysis of the small fluctuations at such long distances shows that there are no significant differences between the ranges of positional and orientational order.

Two-component system CCl4 + (CH3)3CBr: extrema in equilibria involving orientationally disordered phases

Phase equilibria involving orientationally disordered (OD) and liquid phases of the two-component system between carbon tetrachloride (CCl4) and 2-methyl-2-bromomethane ((CH3)3CBr) have been determined by means of X-ray powder diffraction and thermal analysis techniques from 210 K up to the liquid state. The isomorphism relation between the OD stable face-centered cubic (FCC) phase of (CH3)3CBr and the metastable FCC phase of CCl4 has been demonstrated throughout the continuous evolution of the lattice parameters and the existence of the two-phase equilibrium [FCC + L] for the whole range of composition, despite the monotropy of the FCC phase for the CCl4 component with respect to its OD rhombohedral (R) stable phase. A continuous series of OD R mixed crystals is found, which confirms the R lattice symmetry of the OD phase II of (CH3)3CBr, for which the crystallographic results have been long-time misinterpreted. X-ray patterns of such a phase were indexed according to the recent single-crystal results obtained by Rudman (Rudman, R. J. Mol. Struct. 2001, 569, 157). In addition, some experimental evidences are given to confirm the number of molecules per unit cell (Z = 21). The thermodynamic assessment reproduces coherently the phase diagram for the stable [R + L] and [R + FCC] two-phase equilibria as well as for the partially metastable [FCC + L] two-phase equilibrium and provides a set of data for the thermodynamic properties of nonexperimentally available phase transitions of pure components. Surprisingly, the phase equilibrium involving R and FCC OD phases appears as one of the very few showing a solid-solid equilibrium with two extremes.

α and β relaxation dynamics of a fragile plastic crystal

We present a thorough dielectric investigation of the relaxation dynamics of plastic crystalline Freon112, which exhibits freezing of the orientational degrees of freedom into a glassy crystal below 90 K. Among other plastic crystals, Freon112 stands out by being relatively fragile within Angell's [Relaxations in Complex Systems, edited by K. L. Ngai and G. B. Wright (NRL, Washington, DC, 1985), p. 3] classification scheme and by showing an unusually strong beta relaxation. Comparing the results to those on Freon112a, having only a single molecular conformation, points to the importance of the presence of two molecular conformations in Freon112 for the explanation of its unusual properties.