Abstract: We study the structure and phase behavior of fullerene-cubane C60·C8H8 by Monte Carlo simulation. Using a simple potential model capturing the icosahedral and cubic symmetries of its molecular constituents, we reproduce the experimentally observed phase transition from a cubic to an orthorhombic crystal lattice and the accompanying rotational freezing of the C60 molecules. We elaborate a scheme to identify the low-temperature orientations of individual molecules and to detect a pattern of orientational ordering similar to the arrangement of C60 molecules in solid C60. Our configuration of orientations supports a doubled periodicity along one of the crystal axes.

Abstract: Superexchange interactions in alkali fullerides AC(60) are derived for C-60 molecular ions separated by interstitial alkali-metal ions. We use a multiconfiguration approach which comprises the lowest molecular orbital states of the C-60 molecule and the excited s and d states of the alkali-metal atom A. Interactions are described by the valence bond (Heitler-London) method for a complex (C-60 – A – C-60) – with two valence electrons. The electronic charge transfer between the alkali-metal atom and a neighboring C-60 molecule is not complete. The occupation probability of excited d and s states of the alkali atom is not negligible. In correspondence with the relative positions of the C-60 molecules and A atoms in the polymer crystal, we consider 180degrees and 90degrees (angle) superexchange pathways. For the former case the ground state is found to be a spin singlet separated from a triplet at similar to20 K. For T < 20 K there appear strong spin correlations for the 180degrees superexchange pathway. The results are related to spin lattice relaxation experiments on CsC60 in the polymerized and in the quenched cubic phase. (C) 2005 American Institute of Physics.

Abstract: In order to increase the accuracy of the linearized augmented plane wave (LAPW) method, we present a new approach where the plane wave basis function is augmented by two different atomic radial components constructed at two different linearization energies corresponding to two different electron bands (or energy windows). We demonstrate that this case can be reduced to the standard treatment within the LAPW paradigm where the usual basis set is enriched by the basis functions of the tight binding type, which go to zero with zero derivative at the sphere boundary. We show that the task is closely related with the problem of extended core states which is currently solved by applying the LAPW method with local orbitals (LAPW+LO). In comparison with LAPW+LO, the number of supplemented basis functions in our approach is doubled, which opens up a new channel for the extension of the LAPW and LAPW+LO basis sets. The appearance of new supplemented basis functions absent in the LAPW+LO treatment is closely related with the existence of the ul-component in the canonical LAPW method. We discuss properties of additional tight binding basis functions and apply the extended basis set for computation of electron energy bands of lanthanum (face and body centered structures) and hexagonal close packed lattice of cadmium. We demonstrate that the new treatment gives lower total energies in comparison with both canonical LAPW and LAPW+LO, with the energy difference more pronounced for intermediate and poor LAPW basis sets.

Abstract: Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems—a particle solvated in a Lennard-Jones fluid, a Diels–Alder reaction, and crystallization of liquid sodium—to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.

Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.

Abstract: The equilibrium bond length, harmonic frequency, first and second order anharmonicity constants, rotational and centrifugal distortion constants, as well as the rotation-vibrational and centrifugal coupling constants for the ground X(2) Sigma(+) and first excited A(2) Pi states of the SiN radical have been calculated at the complete active space SCF (CASSCF), multireference CI (MRCI) and coupled cluster (CCSD(T)) levels using Dunning's correlation-consistent basis sets. The excitation energy of the A(2) Pi State has also been computed at these theoretical levels. Dipole moments of SiN in the X(2) Sigma(+) and A(2) Pi states are given. Our study shows that core correlation must be considered in order to obtain satisfactory accuracy for the spectroscopic constants.

Abstract: Multi-wall carbon nanotubes have been produced by the catalytic decomposition of acetylene. Go-Mo, Co-V and Co-Fe mixtures supported either on zeolite or corundum alumina were used as catalysts. When Fe or V is added to Co, the carbon deposit increases. The nanotubes were characterized by both low and high resolution TEM. From histograms representing the outer diameter distributions, it is clear that the outer diameter of the nanotubes can be controlled by choosing the appropriate catalyst. (C) 2000 Elsevier Science B.V. All rights reserved.