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: Uniform-acceptance force-bias Monte Carlo (fbMC) methods have been shown to be a powerful technique to access longer timescales in atomistic simulations allowing, for example, phase transitions and growth. Recently, a new fbMC method, the time-stamped force-bias Monte Carlo (tfMC) method, was derived with inclusion of an estimated effective timescale; this timescale, however, does not seem able to explain some of the successes the method. In this contribution, we therefore explicitly quantify the effective timescale tfMC is able to access for a variety of systems, namely a simple single-particle, one-dimensional model system, the Lennard-Jones liquid, an adatom on the Cu(100) surface, a silicon crystal with point defects and a highly defected graphene sheet, in order to gain new insights into the mechanisms by which tfMC operates. It is found that considerable boosts, up to three orders of magnitude compared to molecular dynamics, can be achieved for solid state systems by lowering of the apparent activation barrier of occurring processes, while not requiring any system-specific input or modifications of the method. We furthermore address the pitfalls of using the method as a replacement or complement of molecular dynamics simulations, its ability to explicitly describe correct dynamics and reaction mechanisms, and the association of timescales to MC simulations in general.

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: Single-chain elasticity of polyethylene at theta point up to 90% of stretching with respect to its contour length is computed by Monte Carlo simulation of an atomistic model in continuous space. The elasticity law together with the free-energy and the internal energy variations with stretching are found to be very well represented by the wormlike chain model up to 65% of the chain elongation, provided the persistence length is treated as a temperature-dependent parameter. Beyond this value of elongation simple ideal chain models are not able to describe the Monte Carlo data in a thermodynamic consistent way. This study reinforces the use of the wormlike chain model to interpret experimental data on the elasticity of synthetic polymers in the finite extensibility regime, provided the chain is not yet in its fully stretched regime. Specific solvent effects on the elasticity law and the partition between energetic and entropic contributions to single chain elasticity are investigated. (C) 2002 American Institute of Physics.

Abstract: A microscopic theory, which describes the orientational dynamics of C60 molecules in the face-centered cubic phase of C60-fullerite, is formulated or the case of a complex molecular structure. Interaction centers which comprise atoms, double bonds, and single bonds as molecular constituents contribute to the intermolecular potential. Orientation dependent physical properties are described in terms of symmetry-adapted rotator functions. It is found that a same set of rotator functions is sufficient even in the case of a complex molecular structure. Phase transition temperatures are discussed for various models of molecular structure.