Abstract: Starting from a model of microscopic interactions between C-70 molecules, we have developed a theory which describes the orientational dynamics and its coupling to lattice displacements in the rhombohedral phase of C-70 fullerite. The Landau free energy is calculated. We obtain a first-order phase transition to a monoclinic structure with the space group P2(1)/m. The transition is driven by the condensation of orientational quadrupoles at the F point of the Brillouin zone of the rhombohedral lattice. We find no evidence that the monoclinic structure is connected with the freezing in of orientations around the fivefold molecular axis. We calculate the lattice strains that are associated with the transition to the monoclinic structure. The theory is compared with a range of experimental data on the phase transition.

Abstract: Using first-principles calculations and reactive force field molecular dynamics simulations, we study the structural properties and dynamics of a fluorine (F) atom, either adsorbed on the surface of single layer graphene (F/GE) or between the layers of AB stacked bilayer graphene (F@ bilayer graphene). It is found that the diffusion of the F atom is very different in those cases, and that the mobility of the F atom increases by about an order of magnitude when inserted between two graphene layers. The obtained diffusion constant for F/GE is twice larger than that experimentally found for gold adatom and theoretically found for C-60 molecule on graphene. Our study provides important physical insights into the dynamics of fluorine atoms between and on graphene layers and explains the mechanism behind the separation of graphite layers due to intercalation of F atoms.

Abstract: Anharmonic effects in an atomic monolayer thin crystal with honeycomb lattice structure are investigated by analytical and numerical lattice dynamical methods. Starting from a semiempirical model for anharmonic couplings of third and fourth orders, we study the in-plane and out-of-plane (flexural) mode components of the generalized wave vector dependent Gruneisen parameters, the thermal tension and the thermal expansion coefficients as a function of temperature and crystal size. From the resonances of the displacement-displacement correlation functions, we obtain the renormalization and decay rate of in-plane and flexural phonons as a function of temperature, wave vector, and crystal size in the classical and in the quantum regime. Quantitative results are presented for graphene. There, we find that the transition temperature T-alpha from negative to positive thermal expansion is lowered with smaller system size. Renormalization of the flexural mode has the opposite effect and leads to values of T-alpha approximate to 300 K for systems of macroscopic size. Extensive numerical analysis throughout the Brillouin zone explores various decay and scattering channels. The relative importance of normal and umklapp processes is investigated. The work is complementary to crystalline membrane theory and computational studies of anharmonic effects in two-dimensional crystals.

Abstract: The Majorana zero modes and persistent spin current in mesoscopic d-wave-superconducting loops with spin-orbit (SO) interaction are investigated by numerically solving the spin-generalized Bogoliubov-de Gennes equations self-consistently. For some appropriate strength of the SO coupling, Majorana zero-energy states and sharp jumps of the spin-polarized currents can be observed when the highest energy levels cross the Fermi energy in the spectrum, leading to spin currents with opposite chirality flowing near the inner and outer edges of the sample. When the threaded magnetic flux turns on, four flux-dependent patterns of the persistent spin current with step-like features show up, accompanied by Majorana edge modes at flux values where the energy gap closes. Moreover, the Majorana zero mode is highly influenced by the direction of the Zeeman field. A finite in-plane field can lead to the gap opening since the inversion symmetry is broken. Remarkably, multiple Majorana zero-energy states occur in the presence of an out-of-plane field h(z), and the number of steps in the spin current evolution can be effectively tuned by the field strength due to the shift of Majorana zero modes. Finally, when the loop sample contains surface indentation defects, zero-energy modes can always show up in the presence of an appropriate h(z). Interestingly, multiple Majorana states may be present in the system with a corner defect even if h(z) = 0.

Abstract: In order to study quantum effects in a two-dimensional crystal lattice of a free-standing monolayer graphene, we have performed both path-integral Monte Carlo (PIMC) and classical Monte Carlo (MC) simulations for temperatures up to 2000 K. The REBO potential is used for the interatomic interaction. The total energy, interatomic distance, root-mean-square displacement of the atom vibrations, and the free energy of the graphene layer are calculated. The obtained lattice vibrational energy per atom from the classical MC simulation is very close to the energy of a three-dimensional harmonic oscillator 3k(B)T. The PIMC simulation shows that quantum effects due to zero-point vibrations are significant for temperatures T < 1000 K. The quantum contribution to the lattice vibrational energy becomes larger than that of the classical lattice for T < 400 K. The lattice expansion due to the zero-point motion causes an increase of 0.53% in the lattice parameter. A minimum in the lattice parameter appears at T similar or equal to 500 K. Quantum effects on the atomic vibration amplitude of the graphene lattice and its free energy are investigated.