Abstract: The effect of a parallel magnetic field on the orbital motion of electrons in high-quality superconducting nanowires resulting in a superconductor-to-normal transition which occurs through a cascade of jumps in the order parameter as a function of the magnetic field. Such cascades originate from the transverse size quantization that splits the conduction band into a series of subbands. Here, based on a numerical solution of the Bogoliubov-de Gennes equations for a hollow nanocylinder, we investigate how the quantum-size cascades depend on the confining geometry, i.e., by changing the cylinder radius R and its thickness d we cover the range from the nanowire-like to the nanofilm-like regime. The cascades are shown to become much less pronounced when increasing R/d, i.e., when the nanofilm-like regime is approached. When the temperature is non-zero they are thermally smoothed. This includes the spin-magnetic-field interaction which reduces the critical (depairing) parallel magnetic field H-c,H-parallel to but does not have any qualitative effect on the quantum cascades. From our calculations it is seen that the paramagnetic limiting field H-par significantly exceeds H-c,H-parallel to even in extremely narrow nanocylinders, i.e., when R, d are down to a few nanometers, and H-c,H-parallel to is only about 10% larger when switching-off the spin-magnetic-field interaction in this case. Both characteristic fields, H-c,H-parallel to and H-par, exhibit pronounced quantum-size oscillations. We demonstrate that the quantum cascades and the quantum-size oscillations survive in the presence of surface roughness.

Abstract: Quantum confinement of the perpendicular motion of electrons in single-crystalline metallic superconducting nanofilms splits the conduction band into a series of single-electron subbands. A distinctive feature of such a nanoscale multi-band superconductor is that the energetic position of each subband can vary significantly with changing nanofilm thickness, substrate material, protective cover and other details of the fabrication process. It can occur that the bottom of one of the available subbands is situated in the vicinity of the Fermi level. We demonstrate that the character of the superconducting pairing in such a subband changes dramatically and exhibits a clear molecule-like trend, which is very similar to the well-known crossover from the Bardeen-Cooper-Schrieffer regime to Bose-Einstein condensation (BCS-BEC) observed in trapped ultracold fermions. For Pb nanofilms with thicknesses of 4 and 5 monolayers (MLs) this will lead to a spectacular scenario: up to half of all the Cooper pairs nearly collapse, shrinking in the lateral size (parallel to the nanofilm) down to a few nanometers. As a result, the superconducting condensate will be a coherent mixture of almost molecule-like fermionic pairs with ordinary, extended Cooper pairs.

Abstract: Superconducting correlations in an isolated metallic grain are governed by the interplay between two energy scales: the mean level spacing delta and the bulk pairing gap Delta(0), which are strongly influenced by the position of the chemical potential with respect to the closest single-electron level. In turn superconducting correlations affect the position of the chemical potential. Within the parity projected BCS model we investigate the probability distribution of the chemical potential in a superconducting grain with randomly distributed single-electron levels. Taking into account statistical fluctuations of the chemical potential due to the pairing interaction, we find that such fluctuations have a significant impact on the critical level spacing delta(c) at which the superconducting correlations cease: the critical ratio delta(c)/Delta(0) at which superconductivity disappears is found to be increased.

Abstract: Recent experiments have revealed that the macroscopic Kelvin equation remains surprisingly accurate even for nanoscale capillaries. This phenomenon was so far explained by the oscillatory behavior of the solid−liquid interfacial free energy. We here demonstrate thermodynamic and capillarity inconsistencies with this explanation. After revising the Kelvin equation, we ascribe its validity at nanoscale confinement to the effect of disjoining pressure.
To substantiate our hypothesis, we employed molecular dynamics simulations to evaluate interfacial heat transfer and wetting properties. Our assessments unveil a breakdown in a previously established proportionality between the work of adhesion and the Kapitza conductance at capillary heights below 1.3 nm, where the dominance of the work of adhesion shifts primarily from energy to entropy. Alternatively, the peak density of the initial water layer can effectively probe the work of adhesion. Unlike under bulk conditions, high confinement renders the work of adhesion entropically unfavorable.

Abstract: We investigate the stability of circular monolayer graphene subjected to a radial load using non-equilibrium molecular dynamics simulations. When monolayer graphene is radially stressed, after some small circular strain (~0.4%) it buckles and bends into a new bowl-like shape. Young's modulus is calculated from the linear relation between stress and strain before the buckling threshold, which is in agreement with experimental results. The prediction of elasticity theory for the buckling threshold of a radially stressed plate is presented and its results are compared to the one of our atomistic simulation. The Jarzynski equality is used to estimate the difference between the free energy of the non-compressed states and the buckled states. From a calculation of the free energy we obtain the optimum radius for which the system feels the minimum boundary stress.

Abstract: The electron states in graphene-based magnetic dot and magnetic ring structures and combinations of both are investigated. The corresponding spectra are studied as a function of the radii, the strengths of the inhomogeneous magnetic field and of a uniform background field, the strength of an electrostatic barrier and the angular momentum quantum number. In the absence of an external magnetic field we have only long-lived quasi-bound and scattering states and we assess their influence on the density of states. In addition, we consider elastic electron scattering by a magnetic dot, whose average B vanishes, and show that the Hall and longitudinal resistivities, as a function of the Fermi energy, exhibit a pronounced oscillatory structure due to the presence of quasi-bound states. Depending on the dot parameters this oscillatory structure differs substantially for energies below and above the first Landau level.

Abstract: We study the effect of electron confinement on the superconducting-to-normal phase transition driven by a magnetic field and/or on the current-carrying state of the superconducting condensate in nanowires. Our investigation is based on a self-consistent numerical solution of the Bogoliubov-de Gennes equations. We show that in a parallel magnetic field and/or in the presence of supercurrent the transition from superconducting to normal phase occurs as a cascade of discontinuous jumps in the superconducting order parameter for diameters D < 10 divided by 15 nm at T = 0. The critical magnetic field exhibits quantum-size oscillations with pronounced resonant enhancements.

Abstract: We investigate the dynamics of a charged particle moving in a graphene layer and in a two-dimensional electron gas, where it obeys the Dirac and the Schrödinger equations, respectively. The charge carriers are described as Gaussian wavepackets. The dynamics of the wavepackets is studied numerically by solving both quantum-mechanical and relativistic equations of motion. The scattering of such wavepackets by step-like magnetic and potential barriers is analysed for different values of wavepacket energy and width. We find: (1) that the average position of the wavepacket does not coincide with the classical trajectory, and (2) that, for slanted incidence, the path of the centre of mass of the wavepacket does not have to penetrate the barrier during the scattering process. Trembling motion of the charged particle in graphene is observed in the absence of an external magnetic field and can be enhanced by a substrate-induced mass term.

Abstract: Two-dimensional materials composed of transition metal carbides and nitrides (MXenes) are poised to revolutionize energy conversion and storage. In this work, we used density functional theory (DFT) to investigate the adsorption of Mg and Na adatoms on five M2CS2 monolayers (where M = Mo, Nb, Ti, V, and Zr) for battery applications. We assessed the stability of the adatom (i.e. Na and Mg)-monolayer systems by calculating adsorption and formation energies, as well as voltages as a function of surface coverage. For instance, we found that Mo2CS2 cannot support a full layer of Na nor even a single Mg atom. Na and Mg exhibit the strongest binding on Zr2CS2, followed by Ti2CS2, Nb2CS2 and V2CS2. Using the nudged elastic band method (NEB), we computed promising diffusion barriers for both dilute and nearly full ion surface coverage cases. In the dilute ion adsorption case, a single Mg and Na atom on Ti2CS2 experience similar to 0.47 eV and similar to 0.10 eV diffusion barriers between the lowest energy sites, respectively. For a nearly full surface coverage, a Na ion moving on Ti2CS2 experiences a similar to 0.33 eV energy barrier, implying a concentration-dependent diffusion barrier. Our molecular dynamics results indicate that the three (one) layers (layer) of the Mg (Na) ion on both surfaces of Ti2CS2 remain stable at T = 300 K. While, according to voltage calculations, Zr2CS2 can store Na up to three atomic layers, our MD simulations predict that the outermost layers detach from the Zr2CS2 monolayer due to the weak interaction between Na ions and the monolayer. This suggests that MD simulations are essential to confirm the stability of an ion-electrode system – an insight that is mostly absent in previous studies.

Abstract: Topological properties of a double-stranded DNA (dsDNA) proximity-coupled by an s-wave superconductor are investigated, in which the energy spectra and the differential conductance are calculated within the framework of tight-binding approximation. Our results indicate that this dsDNA-superconductor system hosts Majorana zero modes (MZMs) when the Zeeman field is perpendicular to the helix axis, whereas no MZM could be observed when the Zeeman field is parallel to the helix axis, in sharp contrast to previous studies on nanowires including single-stranded DNA. In particular, two topological phase transitions could take place in the dsDNA-superconductor system by changing the Zeeman field, one from a topological trivial phase to a topological nontrivial phase with one pair of MZMs in small Zeeman field regime, and the other from a phase with one pair of MZMs to a phase with two pairs of MZMs by further increasing the Zeeman field. In the presence of a gate field normal to the helix axis, the topological nontrivial phase with two pairs of MZMs can transform into the phase with one pair of MZMs. The topological phase with one pair of MZMs is more stable and robust against Anderson disorder.

Abstract: The highly tunable band structure of the zero-energy Landau level (zLL) of bilayer graphene makes it an ideal platform for engineering novel quantum states. However, the zero-energy Landau level at high electric fields has remained largely unexplored. Here we present magnetotransport measurements of bilayer graphene in high transverse electric fields. We observe previously undetected Landau level crossings at filling factors nu = -2, 1, and 3 at high electric fields. These crossings provide constraints for theoretical models of the zero-energy Landau level and show that the orbital, valley, and spin character of the quantum Hall states at high electric fields is very different from low electric fields. At high E, new transitions between states at nu = -2 with different orbital and spin polarization can be controlled by the gate bias, while the transitions between nu = 0 -> 1 and nu = 2 -> 3 show anomalous behavior.

Abstract: We study vortex configurations in mesoscopic superconducting squares. Our theoretical approach is based on the analytical solution of the London equation using the Green's function method. The potential energy landscape found is then used in Langevin-type molecular-dynamics simulations to obtain stable vortex configurations. We show that the filling rules for vortices in squares with increasing applied magnetic field can be formulated, although in a different manner than in disks, in terms of the formation of vortex 'shells'. We discuss metastable states and the stability of the vortex configurations found with respect to variations of the material parameters and deformations of the shape of the sample.