Abstract: The spiral magnetic structure and underlying spin lattice of BiMnFe2O6 are investigated by low-temperature neutron powder diffraction and density functional theory band structure calculations. In spite of the random distribution of the Mn3+ and Fe3+ cations, this centrosymmetric compound undergoes a transition into an incommensurate antiferromagnetically ordered state below TN≃220 K. The magnetic structure is characterized by the propagation vector k=[0,β,0] with β≃0.14 and the P221211′(0β0)0s0s magnetic superspace symmetry. It comprises antiferromagnetic helixes propagating along the b axis. The magnetic moments lie in the ac plane and rotate about π(1+β)≃204.8-deg angle between the adjacent magnetic atoms along b. The spiral magnetic structure arises from the peculiar frustrated arrangement of exchange couplings in the ab plane. The antiferromagnetic coupling along the c axis cancels the possible electric polarization and prevents ferroelectricity in BiMnFe2O6.

Abstract: Iron carbides play a crucial role in steel manufacturing and processing and to a large extent determine the physical properties of steel products. The modified embedded atom method (MEAM) in combination with Lee's Fe-C potential is a good candidate for molecular dynamics simulations on larger Fe-C systems. Here, we investigate the stability and crystal structures of pure iron and binary iron carbides using MEAM and compare them with the experimental data and quantum-mechanical density functional theory calculations. The analysis shows that the Fe-C potential gives reasonable results for the relative stability of iron and iron carbides. The performance of MEAM for the prediction of the potential energy and the calculated lattice parameters at elevated temperature for pure iron phases and cementite are investigated as well. The conclusion is that Lee's MEAM Fe-C potential provides a promising basis for further molecular dynamics simulations of Fe-C alloys and steels at lower temperatures (up to 800 K).

Abstract: Using an effective one-dimensional cluster expansion in combination with first-principles electronic structure calculations we have studied the energetics and electronic properties of Sb-Te layered systems. For a Te concentration between 0 and 60 at. % an almost continuous series of metastable structures is obtained consisting of consecutive Sb bilayers next to consecutive Sb2Te3 units, with the general formula (Sb-2)(n)(Sb2Te3)(m) (n, m = 1,2, ... ). Between 60 and 100 at.% no stable structures are found. We account explicitly for the weak van derWaals bonding between Sb bilayers and Sb2Te3 units by using a recently developed functional, which strongly improves the interlayer bonding distances. At T = 0 K, no evidence is found for the existence of two separate single-phase regions delta and gamma and a two-phase region delta + gamma. Metastable compounds with a Te concentration between 0 and 40 at. % are semimetallic, whereas compounds with a Te concentration between 50 and 60 at. % are semiconducting. Compounds with an odd number of Sb layers are metallic and have a much higher formation energy than those with an even number of consecutive Sb layers, thereby favoring the formation of Sb bilayers.

Abstract: We have revealed thermodynamically stable vortex-antivortex configurations in mesoscopic type I superconducting equilateral triangles. The revealed effect is explained by two factors: (i) vortex confinement in mesoscopic triangles and (ii) appearance of the vortex-antivortex repulsion in mesoscopic type I superconductor triangles in the vicinity of the dual point. The stability of the obtained solutions is examined as a function of material and external parameters – the Ginzburg-Landau parameter K and temperature – as well as a function of a distortion of the sample's shape. The vortex-antivortex. “molecule” turns out to be stable in a wide range of the parameters that makes it possible its experimental visualization.

Abstract: The effect of different stacking order of graphene multilayers on the electric field induced band gap is investigated. We considered a positively charged top and a negatively charged back gate in order to independently tune the band gap and the Fermi energy of three and four layer graphene systems. A tight-binding approach within a self-consistent Hartree approximation is used to calculate the induced charges on the different graphene layers. We found that the gap for trilayer graphene with the ABC stacking is much larger than the corresponding gap for the ABA trilayer. Also we predict that for four layers of graphene the energy gap strongly depends on the choice of stacking, and we found that the gap for the different types of stacking is much larger as compared to the case of Bernal stacking. Trigonal warping changes the size of the induced electronic gap by approximately 30% for intermediate and large values of the induced electron density.

Abstract: During the synthesis of ultrathin materials with hexagonal lattice structure Stone-Wales (SW) type of defects are quite likely to be formed and the existence of such topological defects in the graphenelike structures results in dramatic changes of their electronic and mechanical properties. Here we investigate the formation and reactivity of such SW defects in silicene. We report the energy barrier for the formation of SW defects in freestanding (similar to 2.4 eV) and Ag(111)-supported (similar to 2.8 eV) silicene and found it to be significantly lower than in graphene (similar to 9.2 eV). Moreover, the buckled nature of silicene provides a large energy barrier for the healing of the SW defect and therefore defective silicene is stable even at high temperatures. Silicene with SW defects is semiconducting with a direct band gap of 0.02 eV and this value depends on the concentration of defects. Furthermore, nitrogen substitution in SW-defected silicene shows that the defect lattice sites are the least preferable substitution locations for the N atoms. Our findings show the easy formation of SW defects in silicene and also provide a guideline for band gap engineering in silicene-based materials through such defects.

Abstract: Using atomistic simulations we investigate the morphological properties of graphene deposited on top of a nanostructured substrate. Sinusoidally corrugated surfaces, steps, elongated trenches, one-dimensional and cubic barriers, spherical bubbles, Gaussian bumps, and Gaussian depressions are considered as support structures for graphene. The graphene-substrate interaction is governed by van der Waals forces and the profile of the graphene layer is determined by minimizing the energy using molecular dynamics simulations. Based on the obtained optimum configurations, we found that (i) for graphene placed over sinusoidally corrugated substrates with corrugation wavelengths longer than 2 nm, the graphene sheet follows the substrate pattern while for supported graphene it is always suspended across the peaks of the substrate, (ii) the conformation of graphene to the substrate topography is enhanced when increasing the energy parameter in the van der Waals model, (iii) the adhesion of graphene into the trenches depends on the width of the trench and on the graphene's orientation, i. e., in contrast to a small-width (3 nm) nanoribbon with armchair edges, the one with zigzag edges follows the substrate profile, (iv) atomic-scale graphene follows a Gaussian bump substrate but not the substrate with a Gaussian depression, and (v) the adhesion energy due to van der Waals interaction varies in the range [0.1-0.4] J/m(2).

Abstract: We present a tight-binding investigation of strained bilayer graphene within linear elasticity theory, focusing on the different environments experienced by the A and B carbon atoms of the different sublattices. We find that the inequivalence of the A and B atoms is enhanced by the application of perpendicular strain epsilon(zz), which provides a physical mechanism for opening a band gap, most effectively obtained when pulling the two graphene layers apart. In addition, perpendicular strain introduces electron-hole asymmetry and can result in linear electronic dispersion near the K point. Our findings suggest experimental means for strain-engineered band gaps in bilayer graphene.

Abstract: We study the influence of the inelastic relaxation time τ̃E of the quasiparticle distribution function f(E) on the phase slip process in quasi-one-dimensional superconducting rings at a temperature close to the critical temperature Tc. We find that the initial time of growth of the order parameter |Δ| in the phase slip core after the phase slip is a nonmonotonic function of τ̃E which has a maximum at τ̃E≃τ̃GL=πℏ/8kB(Tc−T) and has a tendency to saturate for large τ̃E⪢τ̃GL. The effective heating of the electron subsystem due to the increase in |Δ| in the phase slip center together with the above effect result in a nonmonotonic dependence of the number of subsequent phase slips on τ̃E in rings of relatively large radius (in which each phase slip reduces the current density to a small fraction of its initial value). During the phase slip process the order parameter distribution has two peaks near the phase slip core due to the diffusion of the nonequilibrium quasiparticles from that region.

Abstract: Structural phase transitions in two vertically or horizontally coupled channels of strongly interacting particles are investigated. The particles are free to move in the x direction but are confined by a parabolic potential in the y direction. They interact with each other through a screened power-law potential (r(-n)e(-r/lambda)). In vertically coupled systems, the channels are stacked above each other in the direction perpendicular to the (x, y) plane, while in horizontally coupled systems both channels are aligned in the confinement direction. Using Monte Carlo (MC) simulations we obtain the ground-state configurations and the structural transitions as a function of the linear particle density and the separation between the channels. At zero temperature, the vertically coupled system exhibits a rich phase diagram with continuous and discontinuous transitions. On the other hand, the horizontally coupled system exhibits only a very limited number of phase transitions due to its symmetry. Further, we calculated the normal modes for the Wigner crystals in both cases. From MC simulations, we found that in the case of vertically coupled systems, the zigzag transition is only possible for low densities. A Ginzburg-Landau theory for the zigzag transition is presented, which predicts correctly the behavior of this transition from which we interpret the structural phase transition of the Wigner crystal through the reduction of the Brillouin zone.

Abstract: In this paper, we report on the electronic transport properties of mesoscopic, ultrathin graphite flakes with a thickness corresponding to a stack of 150 graphene layers. The graphite flakes show an unexpectedly strong positive magnetoresistance (PMR) already at room temperature, which scales in good approximation with the square of the magnetic field. Furthermore, we show that the resistivity is unaffected by magnetic fields oriented in plane with the graphene layers. Hall effect measurements indicate that the charge carriers are p-type and their concentration increases with increasing temperature while the mobility is decreasing. The Hall voltage is non-linear in higher magnetic fields. Possible origins of the observed effects are discussed. Ball and stick model of the two topmost carbon layers of the hexagonal graphite structure.

Abstract: The interplay between quantum Hall states and Cooper pairs is usually hindered by the suppression of the superconducting state due to the strong magnetic fields needed to observe the quantum Hall effect. From this point of view, graphene is special since it allows the creation of strong pseudomagnetic fields due to strain. We show that in a Josephson junction made of strained graphene, Cooper pairs will diffuse into the strained region. The pair correlation function will be sublattice polarized due to the polarization of the local density of states in the zero pseudo-Landau level. We uncover two regimes: (1) one in which the cyclotron radius is larger than the junction length, in which case the supercurrent will be enhanced, and (2) the long junction regime where the supercurrent is strongly suppressed because the junction becomes an insulator. In the latter case quantized Hall states form and Andreev scattering at the normal/superconducting interface will induce edge states. Our numerical calculation has become possible due to an extension of the Chebyshev-Bogoliubovde Gennes method to computations on video cards (GPUs).