Abstract: The Coulomb interaction and the correlation of a remote electron with a single layer of graphene is investigated in the presence of a magnetic field applied perpendicular to the graphene layer. The remote electron polarizes the electron gas in the graphene layer, which we describe in terms of excitations of virtual plasmons in graphene. The composite quasiparticle formed by electron plus polarization is called a plasmon polaron. The ground-state energy of this quasiparticle is calculated within perturbation theory for remote electrons in different environments.

Abstract: We present a comprehensive study of the crystal structure, magnetic structure, and microscopic magnetic model of (CuBr)LaNb2O7, the Br analog of the spin-gap quantum magnet (CuCl) LaNb2O7. Despite similar crystal structures and spin lattices, the magnetic behavior and even peculiarities of the atomic arrangement in the Cl and Br compounds are very different. The high- resolution x-ray and neutron data reveal a split position of Br atoms in (CuBr) LaNb2O7. This splitting originates from two possible configurations developed by [CuBr] zigzag ribbons. While the Br atoms are locally ordered in the ab plane, their arrangement along the c direction remains partially disordered. The predominant and energetically more favorable configuration features an additional doubling of the c lattice parameter that was not observed in (CuCl) LaNb2O7. (CuBr) LaNb2O7 undergoes long-range antiferromagnetic ordering at T-N = 32 K, which is nearly 70% of the leading exchange coupling J4 similar or equal to 48 K. The Br compound does not show any experimental signatures of low-dimensional magnetism because the underlying spin lattice is three-dimensional. The coupling along the c direction is comparable to the couplings in the ab plane, even though the shortest Cu-Cu distance along c (11.69 angstrom) is three times larger than nearest-neighbor distances in the ab plane (3.55 angstrom). The stripe antiferromagnetic long-range order featuring columns of parallel spins in the ab plane and antiparallel spins along c is verified experimentally and confirmed by the microscopic analysis.

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 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: To date, several experimental groups reported measurements of the thickness dependence of T-c of atomically uniform single-crystalline Pb nanofilms. The reported amplitude of the T-c oscillations varies significantly from one experiment to another. Here we propose that the reason for this unresolved issue is an interplay of the quantum-size variations in the single-electron density of states with thickness-dependent oscillations in the phonon-mediated electron-electron coupling. Such oscillations in the coupling depend on the substrate material, the quality of the interface, the protection cover, and other details of the fabrication process, changing from one experiment to another. This explains why the available data do not exhibit one-voice consistency about the amplitude of the T-c oscillations. Our analyses are based on a numerical solution of the Bogoliubov-de Gennes equations for a superconducting slab.

Abstract: The lattice dynamics of low-frequency rigid-plane modes in metallic (graphene multilayers, GML) and in insulating (hexagonal boron-nitride multilayers, BNML) layered crystals is investigated. The frequencies of shearing and compression (stretching) modes depend on the layer number N and are presented in the form of fan diagrams. The results for GML and BNML are very similar. In both cases, only the interactions (van der Waals and Coulomb) between nearest-neighbor planes are effective, while the interactions between more distant planes are screened. A comparison with recent Raman scattering results on low-frequency shear modes in GML [Tan et al., Nat. Mater., in press, doi: 10.1038/nmat3245, (2012)] is made. Relations with the low-lying rigid-plane phonon dispersions in the bulk materials are established. Master curves, which connect the fan diagram frequencies for any given N, are derived. Static and dynamic thermal correlation functions for rigid-layer shear and compression modes are calculated. The results might be of use for the interpretation of friction force experiments on multilayer crystals.

Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.

Abstract: We consider two possible mechanisms for single-vortex fluctuative entry/exit through the surface barrier in ultrathin superconducting disk-shaped nanoislands made of Pb and consisting of just a few monoatomic layers, which can be fabricated using modern techniques. We estimate tunneling probabilities and establish criteria for the crossover between these two mechanisms depending on magnetic field and system sizes. For the case of vortex entry, quantum tunneling dominates on the major part of the temperature/flux phase diagram. For the case of vortex exit, thermal activation turns out to be more probable. This nontrivial result is due to the subtle balance between the barrier height and width, which determine rates of the thermal activation and quantum tunneling, respectively.

Abstract: The Southwest Coastal people of Bangladesh have introduced Tidal River Management (TRM) as an environmentally acceptable water resource management practice based on their indigenous knowledge of water logging of low lying coastal land. TRM helps to address problems resulting from different anthropogenic and structural development activities, and it has been successful in helping coastal communities to adapt to climate change and rising sea level vulnerability by forming new land in Tidal Basins. Hence, it is essential to measure sustainability impacts of TRM from the environmental, socio-economic and institutional perspectives. Therefore, firstly, the study identifies sustainability indicators of TRM considering ecosystem services and secondly, develops an inclusive conceptual framework to understand the important impacts of each indicator at various spatial and temporal scales. The conceptual framework is followed by the construction of a Sustainability Index of Tidal River Management (SITRM). It has advantages over the Ramsar Convention framework (2007) and the World Meteorological Organization (WMO) framework (2012) to measure water sustainability as it includes a sustainable model to project future vulnerability of the community, river and Tidal Basin, emphasizing on climate change issues. It also involves trade-offs analysis, livelihood analysis and SWOT (Strengths, Weaknesses, Opportunities and Threats) analysis for a complete impact assessment to enable decision-makers to focus on those services most likely to be of risks and weaknesses or opportunities and strengths for the sustainability of TRM. Moreover, the framework is a useful guide for policymakers in identifying the sustainability impacts of TRM so that they can choose best coping strategies for coastal people to effectively deal with adverse effects of water-logging and undesired climatic events as well as environmental and socio-economic changes in coastal areas.

Abstract: A shift from batch processing towards continuous processing is of interest in the pharmaceutical industry. However, this transition requires detailed knowledge and process understanding of all consecutive unit operations in a continuous manufacturing line to design adequate control strategies. This can be facilitated by developing mechanistic models of the multi-phase systems in the process. Since modelling efforts only started recently in this field, uncertainties about the model predictions are generally neglected. However, model predictions have an inherent uncertainty (i.e. prediction uncertainty) originating from uncertainty in input data, model parameters, model structure, boundary conditions and software. In this paper, the model prediction uncertainty is evaluated for a model describing the continuous drying of single pharmaceutical wet granules in a six-segmented fluidized bed drying unit, which is part of the full continuous from-powder-to-tablet manufacturing line (Consigma (TM), GEA Pharma Systems). A validated model describing the drying behaviour of a single pharmaceutical granule in two consecutive phases is used. First of all, the effect of the assumptions at the particle level on the prediction uncertainty is assessed. Secondly, the paper focuses on the influence of the most sensitive parameters in the model. Finally, a combined analysis (particle level plus most sensitive parameters) is performed and discussed. To propagate the uncertainty originating from the parameter uncertainty to the model output, the Generalized Likelihood Uncertainty Estimation (GLUE) method is used. This method enables a modeller to incorporate the information obtained from the experimental data in the assessment of the uncertain model predictions and to find a balance between model performance and data precision. A detailed evaluation of the obtained uncertainty analysis results is made with respect to the model structure, interactions between parameters and uncertainty boundaries. (C) 2013 Elsevier B.V. All rights reserved.

Abstract: Plasma-catalytic dry reforming of CH4 (DRM) is promising to convert the greenhouse gasses CH4 and CO2 into value-added chemicals, thus simultaneously providing an alternative to fossil resources as feedstock for the chemical industry. However, while many experiments have been dedicated to plasma-catalytic DRM, there is no consensus yet in literature on the optimal choice of catalyst for targeted products, because the underlying mechanisms are far from understood. Indeed, plasma catalysis is very complex, as it encompasses various chemical and physical interactions between plasma and catalyst, which depend on many parameters. This complexity hampers the comparison of experimental results from different studies, which, in our opinion, is an important bottleneck in the further development of this promising research field. Hence, in this perspective paper, we describe the important physical and chemical effects that should be accounted for when designing plasma-catalytic experiments in general, highlighting the need for standardized experimental setups, as well as careful documentation of packing properties and reaction conditions, to further advance this research field. On the other hand, many parameters also create many windows of opportunity for further optimizing plasma-catalytic systems. Finally, various experiments also reveal the lack of improvement in plasma catalysis compared to plasma-only, specifically for DRM, but the underlying mechanisms are unclear. Therefore, we present our newly developed coupled plasma-surface kinetics model for DRM, to provide more insight in the underlying reasons. Our model illustrates that transition metal catalysts can adversely affect plasmacatalytic DRM, if radicals dominate the plasma-catalyst interactions. Thus, we demonstrate that a good understanding of the plasma-catalyst interactions is crucial to avoiding conditions at which these interactions negatively affect the results, and we provide some recommendations for improvement. For instance, we believe that plasma-catalytic DRM may benefit more from higher reaction temperatures, at which vibrational excitation can enhance the surface reactions.

Abstract: At the interface of electrostatic potential kink profiles, one-dimensional chiral states are found in bilayer graphene (BLG). Such structures can be created by applying an asymmetric potential to the upper and the lower layers of BLG. We found the following: (i) due to the strong confinement by the single kink profile, the unidirectional states are only weakly affected by a magnetic field; (ii) increasing the smoothness of the kink potential results in additional bound states, which are topologically different from those chiral states; and (iii) in the presence of a kink-antikink potential, the overlap between the oppositely moving chiral states results in the appearance of crossing and anticrossing points in the energy spectrum. This leads to the opening of tunable minigaps in the spectrum of the unidirectional topological states.

Abstract: Using the Ginzburg-Landau theory extended to the next-to-leading order, we determine numerically the healing lengths of the two order parameters at the two-gap superconductor/normal metal interface. We demonstrate on several examples that those can be different even in the strict domain of applicability of the Ginzburg-Landau theory. This justifies the use of this theory to describe relevant physics of two-gap superconductors, distinguishing them from their single-gap counterparts. The calculational degree of complexity increases only slightly with respect to the conventional Ginzburg-Landau expansion, thus the extended Ginzburg-Landau model remains numerically far less demanding compared to the full microscopic approaches.