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: Chemisorption of hydrogen on graphene is studied using atomistic simulations with the second generation of reactive empirical bond order Brenner inter-atomic potential. The lowest energy adsorption sites and the most important metastable sites are determined. The H concentration is varied from a single H atom, to clusters of H atoms up to full coverage. We found that when two or more H atoms are present, the most stable configurations of H chemisorption on a single graphene layer are ortho hydrogen pairs adsorbed on one side or on both sides of the graphene sheet. The latter has the highest hydrogen binding energy. The next stable configuration is the orthopara pair combination, and then para hydrogen pairs. The structural changes of graphene caused by chemisorbed hydrogen are discussed and are compared with existing experimental data and other theoretical calculations. The obtained results will be useful for nanoengineering of graphene by hydrogenation and for hydrogen storage.

Abstract: We studied the structural and spectral properties of a classical system consisting of a finite number of particles, moving in two dimensions, and interacting through a repulsive logarithmic potential and held together by an anisotropic harmonic potential. Increasing the anisotropy of the confinement potential can drive the system from a two-dimensional (2D) to a one-dimensional (1D) configuration. This change occurs through a sequence of structural transitions of first and second order which are reflected in the normal mode frequencies. Our results of the ground state configurations are compared with recent experiments and we obtained a satisfactory agreement. The transition from the 1D line structure to the 2D structure occurs through a zigzag transition which is of second order. We found analytical expressions for the eigenfrequencies before the zigzag transition, which allowed us to obtain an analytical expression for the anisotropy parameter at which the zigzag transition occurs as a function of the number of particles in the system.

Abstract: We study the mixing of two different kinds of particles, having different charge and/or mass, interacting through a pure Coulomb potential, and confined in a parabolic trap. The structure of the cluster and its normal mode spectrum are analyzed as a function of the ratio of the charges (mass ratio) of the two types of particles. We show that particles are not always arranged in a shell structure. Mixing of the particles goes hand in hand with a large number of metastable states. The normal modes of the system are obtained, and we find that some of the special modes can be tuned by varying the ratio between the charges (masses) of the two species. The degree of mixing of the two type of particles is summarized in a phase diagram, and an order parameter that describes quantitatively the mixing between particles is defined.

Abstract: Animal models are being used extensively in pre-clinical and safety assessment studies to assess the effectiveness and safety of new chemical entities and delivery systems. Although never entirely replacing the need for animal testing, the use of computer simulations could eventually reduce the amount of animals needed for research purposes and refine the data acquired from the animal studies. Computational fluid dynamics is a powerful tool that makes it possible to simulate flow and particle behavior in animal or patient-specific respiratory models, for purposes of inhaled delivery. This tool requires an accurate representation of the respiratory system, respiration and dose delivery attributes. The aim of this study is to develop a representative airway model of the Sprague-Dawley rat using static and dynamic micro-CT scans. The entire respiratory tract was modeled, from the snout and nares down to the central airways at the point where no distinction could be made between intraluminal air and the surrounding tissue. For the selection of the representative model, variables such as upper airway movement, segmentation length, airway volume and size are taken into account. Dynamic scans of the nostril region were used to illustrate the characteristic morphology of this region in anaesthetized animals. It could be concluded from this study that it was possible to construct a highly detailed representative model of a Sprague-Dawley rat based on imaging modalities such as micro-CT scans

Abstract: Unidirectional chiral states are predicted in single layer graphene which originate from the breaking of the sublattice symmetry due to an asymmetric mass potential. The latter can be created experimentally using boron-nitride (BN) substrates with a line defect (B-B or N-N) that changes the induced mass potential in graphene. Solving the Dirac-Weyl equation, the obtained energy spectrum is compared with the one calculated using ab initio density functional calculations. We found that these one-dimensional chiral states are very robust and they can even exist in the presence of a small gap between the mass regions. In the latter case additional bound states are found that are topologically different from those chiral states.

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: We show that a superconducting quantum nanowire undergoes a new type of BCS-BEC crossover each time when an electron subband approaches the Fermi surface. In this case the longitudinal Cooper-pair size drops by two-three orders of magnitude down to a few nanometers. This unconventional BCS-BEC crossover is driven by quantum-size effects rather than by tuning the fermion-fermion interaction.

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 held exhibits quantum-size oscillations with pronounced resonant enhancements.

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: 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).