Showing 1–30 of 30 results for author: Kaxiras, E

Plasmonic polarization sensing of electrostatic superlattice potentials

Authors: Shuai Zhang, Jordan Fonseca, Daniel Bennett, Zhiyuan Sun, Junhe Zhang, Ran Jing, Suheng Xu, Leo He, S. L. Moore, S. E. Rossi, Dmitry Ovchinnikov, David Cobden, Pablo. Jarillo-Herrero, M. M. Fogler, Philip Kim, Efthimios Kaxiras, Xiaodong Xu, D. N. Basov

Abstract: Plasmon polaritons are formed by coupling light with delocalized electrons. The half-light and half-matter nature of plasmon polaritons endows them with unparalleled tunability via a range of parameters, such as dielectric environments and carrier density. Therefore, plasmon polaritons are expected to be tuned when in proximity to polar materials since the carrier density is tuned by an electrosta… ▽ More Plasmon polaritons are formed by coupling light with delocalized electrons. The half-light and half-matter nature of plasmon polaritons endows them with unparalleled tunability via a range of parameters, such as dielectric environments and carrier density. Therefore, plasmon polaritons are expected to be tuned when in proximity to polar materials since the carrier density is tuned by an electrostatic potential; conversely, the plasmon polariton response might enable the sensing of polarization. Here, we use infrared nano-imaging and nano-photocurrent measurements to investigate heterostructures composed of graphene and twisted hexagonal boron nitride (t-BN), with alternating polarization in a triangular network of moiré stacking domains. We observe that the carrier density and the corresponding plasmonic response of graphene are modulated by polar domains in t-BN. In addition, we demonstrate that the nanometer-wide domain walls of graphene moirés superlattices, created by the polar domains of t-BN, provide momenta to assist the plasmonic excitations. Furthermore, our studies establish that the plasmon of graphene could function as a delicate sensor for polarization textures. The evolution of polarization textures in t-BN under uniform electric fields is tomographically examined via plasmonic imaging. Strikingly, no noticeable polarization switching is observed under applied electric fields up to 0.23 V/nm, at variance with transport reports. Our nano-images unambiguously reveal that t-BN with triangular domains acts like a ferrielectric, rather than ferroelectric claimed by many previous studies. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: 41 pages, 20 figures

Alloy Informatics through Ab Initio Charge Density Profiles: Case Study of Hydrogen Effects in Face-Centered Cubic Crystals

Authors: Dario Massa, Efthimios Kaxiras, Stefanos Papanikolaou

Abstract: Materials design has traditionally evolved through trial-error approaches, mainly due to the non-local relationship between microstructures and properties such as strength and toughness. We propose 'alloy informatics' as a machine learning based prototype predictive approach for alloys and compounds, using electron charge density profiles derived from first-principle calculations. We demonstrate t… ▽ More Materials design has traditionally evolved through trial-error approaches, mainly due to the non-local relationship between microstructures and properties such as strength and toughness. We propose 'alloy informatics' as a machine learning based prototype predictive approach for alloys and compounds, using electron charge density profiles derived from first-principle calculations. We demonstrate this framework in the case of hydrogen interstitials in face-centered cubic crystals, showing that their differential electron charge density profiles capture crystal properties and defect-crystal interaction properties. Radial Distribution Functions (RDFs) of defect-induced differential charge density perturbations highlight the resulting screening effect, and, together with hydrogen Bader charges, strongly correlate to a large set of atomic properties of the metal species forming the bulk crystal. We observe the spontaneous emergence of classes of charge responses while coarse-graining over crystal compositions. Nudge-Elastic-Band calculations show that RDFs and charge features also connect to hydrogen migration energy barriers between interstitial sites. Unsupervised machine-learning on RDFs supports classification, unveiling compositional and configurational non-localities in the similarities of the perturbed densities. Electron charge density perturbations may be considered as bias-free descriptors for a large variety of defects. △ Less

Submitted 24 September, 2023; originally announced September 2023.

Twisto-electrochemical activity volcanoes in Trilayer Graphene

Authors: Mohammad Babar, Ziyan Zhu, Rachel Kurchin, Efthimios Kaxiras, Venkatasubramanian Viswanathan

Abstract: In this work, we develop a twist-dependent electrochemical activity map, combining a tight-binding electronic structure model with modified Marcus-Hush-Chidsey kinetics in trilayer graphene. We identify a counterintuitive rate enhancement region spanning the magic angle curve and incommensurate twists of the system geometry. We find a broad activity peak with a ruthenium hexamine redox couple in r… ▽ More In this work, we develop a twist-dependent electrochemical activity map, combining a tight-binding electronic structure model with modified Marcus-Hush-Chidsey kinetics in trilayer graphene. We identify a counterintuitive rate enhancement region spanning the magic angle curve and incommensurate twists of the system geometry. We find a broad activity peak with a ruthenium hexamine redox couple in regions corresponding to both magic angles and incommensurate angles, a result qualitatively distinct from the twisted bilayer case. Flat bands and incommensurability offer new avenues for reaction rate enhancements in electrochemical transformations. △ Less

Submitted 31 May, 2023; originally announced June 2023.

Comments: 6 pages, 4 figures, Supporting Information

Journal ref: J. Am. Chem. Soc. 2024, 146, 23, 16105-16111

Highly-confined and tunable plasmonics based on two-dimensional solid-state defect lattices

Authors: Ali Ghorashi, Nicholas Rivera, Bowen Shi, Ravishankar Sundararaman, Efthimios Kaxiras, John Joannopoulos, Marin Soljacic

Abstract: Plasmons, collective excitations of electrons in solids, are associated with strongly confined electromagnetic fields, with wavelengths far below the wavelength of photons in free space. This strong confinement promises the realization of optoelectronic devices that could bridge the size difference between photonic and electronic devices. However, despite decades of research in plasmonics, many ap… ▽ More Plasmons, collective excitations of electrons in solids, are associated with strongly confined electromagnetic fields, with wavelengths far below the wavelength of photons in free space. This strong confinement promises the realization of optoelectronic devices that could bridge the size difference between photonic and electronic devices. However, despite decades of research in plasmonics, many applications remain limited by plasmonic losses, thus motivating a search for new engineered plasmonic materials with lower losses. A promising pathway for low-loss plasmonic materials is the engineering of materials with flat and energetically isolated metallic bands, which can strongly limit phonon-assisted optical losses, a major contributor to short plasmonic lifetimes. Such electronic band structures may be created by judiciously introducing an ordered lattice of defects in an insulating host material. Here, we explore this approach, presenting several low-loss, highly-confined, and tunable plasmonic materials based on arrays of carbon substitutions in hexagonal boron nitride (hBN) monolayers. From our first-principles calculations based on density functional theory (DFT), we find plasmonic structures with mid-infrared plasmons featuring very high confinements ($λ_{\text{vacuum}}/λ_{\text{plasmon}}$ exceeding 2000) and quality factors in excess of 1000. We provide a systematic explanation of how crystal structure, electronic bandwidth, and many-body effects affect the plasmonic dispersions and losses of these materials. The results are thus of relevance to low-loss plasmon engineering in other flat band systems. △ Less

Submitted 1 May, 2023; originally announced May 2023.

First principles physics-informed neural network for quantum wavefunctions and eigenvalue surfaces

Authors: Marios Mattheakis, Gabriel R. Schleder, Daniel T. Larson, Efthimios Kaxiras

Abstract: Physics-informed neural networks have been widely applied to learn general parametric solutions of differential equations. Here, we propose a neural network to discover parametric eigenvalue and eigenfunction surfaces of quantum systems. We apply our method to solve the hydrogen molecular ion. This is an ab-initio deep learning method that solves the Schrodinger equation with the Coulomb potential… ▽ More Physics-informed neural networks have been widely applied to learn general parametric solutions of differential equations. Here, we propose a neural network to discover parametric eigenvalue and eigenfunction surfaces of quantum systems. We apply our method to solve the hydrogen molecular ion. This is an ab-initio deep learning method that solves the Schrodinger equation with the Coulomb potential yielding realistic wavefunctions that include a cusp at the ion positions. The neural solutions are continuous and differentiable functions of the interatomic distance and their derivatives are analytically calculated by applying automatic differentiation. Such a parametric and analytical form of the solutions is useful for further calculations such as the determination of force fields. △ Less

Submitted 19 November, 2022; v1 submitted 8 November, 2022; originally announced November 2022.

Boosting the efficiency of ab initio electron-phonon coupling calculations through dual interpolation

Authors: Anderson S. Chaves, Alex Antonelli, Daniel T. Larson, Efthimios Kaxiras

Abstract: The coupling between electrons and phonons in solids plays a central role in describing many phenomena, including superconductivity and thermoelecric transport. Calculations of this coupling are exceedingly demanding as they necessitate integrations over both the electron and phonon momenta, both of which span the Brillouin zone of the crystal, independently. We present here an ab initio method fo… ▽ More The coupling between electrons and phonons in solids plays a central role in describing many phenomena, including superconductivity and thermoelecric transport. Calculations of this coupling are exceedingly demanding as they necessitate integrations over both the electron and phonon momenta, both of which span the Brillouin zone of the crystal, independently. We present here an ab initio method for efficiently calculating electron-phonon mediated transport properties by dramatically accelerating the computation of the double integrals with a dual interpolation technique that combines maximally localized Wannier functions with symmetry-adapted plane waves. The performance gain in relation to the current state-of-the-art Wannier-Fourier interpolation is approximately 2n_s \times M, where n_s is the number of crystal symmetry operations and M, a number in the range 5 - 60, governs the expansion in star functions. We demonstrate with several examples how our method performs some ab initio calculations involving electron-phonon interactions. △ Less

Submitted 10 July, 2020; v1 submitted 30 June, 2020; originally announced June 2020.

Comments: 9 pages, 5 figures

Journal ref: Phys. Rev. B 102, 125116 (2020)

Intrinsic 1T' phase induced in atomically thin 2H-MoTe$_2$ by a single terahertz pulse

Authors: Jiaojian Shi, Ya-Qing Bie, Alfred Zong, Shiang Fang, Wei Chen, Jinchi Han, Zhaolong Cao, Yong Zhang, Takashi Taniguchi, Kenji Watanabe, Vladimir Bulović, Efthimios Kaxiras, Edoardo Baldini, Pablo Jarillo-Herrero, Keith A. Nelson

Abstract: Polymorphic transitions in layered transition metal dichalcogenides provide an excellent platform for discovering exotic phenomena associated with metastable states, ranging from topological phase transitions to enhanced superconductivity. In particular, the transition from 2H to 1T'-MoTe$_2$, which was thought to be induced by high-energy photon irradiation among many other means, has been intens… ▽ More Polymorphic transitions in layered transition metal dichalcogenides provide an excellent platform for discovering exotic phenomena associated with metastable states, ranging from topological phase transitions to enhanced superconductivity. In particular, the transition from 2H to 1T'-MoTe$_2$, which was thought to be induced by high-energy photon irradiation among many other means, has been intensely studied for its technological relevance in nanoscale transistors. Despite the remarkable electrical performance arising from this 2H-to-1T' transition, it remains controversial whether a crystalline 1T' phase is produced because optical signatures of this putative transition are found to be associated with the formation of elemental Te clusters instead. Here, we demonstrate the creation of an intrinsic 1T' lattice after irradiating a mono- or few-layer 2H-MoTe$_2$ with a single field-enhanced terahertz pulse, whose low photon energy limits possible structural damage by optical pulses. To visualize the temporal evolution of this irreversible transition, we further develop a single-shot terahertz pump-second harmonic probe technique, and we find that the transition out of the 2H phase occurs within 10 ns after photoexcitation. Our results not only resolve the long-standing debate over the light-induced polymorphic transition in MoTe$_2$, they also highlight the unique capability of strong-field terahertz pulses in manipulating the structure of quantum materials. △ Less

Submitted 18 October, 2022; v1 submitted 29 October, 2019; originally announced October 2019.

Comments: 10 pages, 3 figures

Journal ref: Nat. Commun. 14, 5905 (2023)

Predicting outcomes of catalytic reactions using machine learning

Authors: Trevor David Rhone, Robert Hoyt, Christopher R. O'Connor, Matthew M. Montemore, Challa S. S. R. Kumar, Cynthia M. Friend, Efthimios Kaxiras

Abstract: Predicting the outcome of a chemical reaction using efficient computational models can be used to develop high-throughput screening techniques. This can significantly reduce the number of experiments needed to be performed in a huge search space, which saves time, effort and expense. Recently, machine learning methods have been bolstering conventional structure-activity relationships used to advan… ▽ More Predicting the outcome of a chemical reaction using efficient computational models can be used to develop high-throughput screening techniques. This can significantly reduce the number of experiments needed to be performed in a huge search space, which saves time, effort and expense. Recently, machine learning methods have been bolstering conventional structure-activity relationships used to advance understanding of chemical reactions. We have developed a model to predict the products of catalytic reactions on the surface of oxygen-covered and bare gold using machine learning. Using experimental data, we developed a machine learning model that maps reactants to products, using a chemical space representation. This involves predicting a chemical space value for the products, and then matching this value to a molecular structure chosen from a database. The database was developed by applying a set of possible reaction outcomes using known reaction mechanisms. Our machine learning approach complements chemical intuition in predicting the outcome of several types of chemical reactions. In some cases, machine learning makes correct predictions where chemical intuition fails. We achieve up to 93% prediction accuracy for a small data set of less than two hundred reactions. △ Less

Submitted 28 August, 2019; originally announced August 2019.

Graphene epsilon-near-zero plasmonic crystals

Authors: Marios Mattheakis, Matthias Maier, Wei Xi Boo, Efthimios Kaxiras

Abstract: Plasmonic crystals are a class of optical metamaterials that consist of engineered structures at the sub-wavelength scale. They exhibit optical properties that are not found under normal circumstances in nature, such as negative-refractive-index and epsilon-near-zero (ENZ) behavior. Graphene-based plasmonic crystals present linear, elliptical, or hyperbolic dispersion relations that exhibit ENZ be… ▽ More Plasmonic crystals are a class of optical metamaterials that consist of engineered structures at the sub-wavelength scale. They exhibit optical properties that are not found under normal circumstances in nature, such as negative-refractive-index and epsilon-near-zero (ENZ) behavior. Graphene-based plasmonic crystals present linear, elliptical, or hyperbolic dispersion relations that exhibit ENZ behavior, normal or negative-index diffraction. The optical properties can be dynamically tuned by controlling the operating frequency and the doping level of graphene. We propose a construction approach to expand the frequency range of the ENZ behavior. We demonstrate how the combination of a host material with an optical Lorentzian response in combination with a graphene conductivity that follows a Drude model leads to an ENZ condition spanning a large frequency range. △ Less

Submitted 31 May, 2019; originally announced June 2019.

Physical Symmetries Embedded in Neural Networks

Authors: M. Mattheakis, P. Protopapas, D. Sondak, M. Di Giovanni, E. Kaxiras

Abstract: Neural networks are a central technique in machine learning. Recent years have seen a wave of interest in applying neural networks to physical systems for which the governing dynamics are known and expressed through differential equations. Two fundamental challenges facing the development of neural networks in physics applications is their lack of interpretability and their physics-agnostic design… ▽ More Neural networks are a central technique in machine learning. Recent years have seen a wave of interest in applying neural networks to physical systems for which the governing dynamics are known and expressed through differential equations. Two fundamental challenges facing the development of neural networks in physics applications is their lack of interpretability and their physics-agnostic design. The focus of the present work is to embed physical constraints into the structure of the neural network to address the second fundamental challenge. By constraining tunable parameters (such as weights and biases) and adding special layers to the network, the desired constraints are guaranteed to be satisfied without the need for explicit regularization terms. This is demonstrated on upervised and unsupervised networks for two basic symmetries: even/odd symmetry of a function and energy conservation. In the supervised case, the network with embedded constraints is shown to perform well on regression problems while simultaneously obeying the desired constraints whereas a traditional network fits the data but violates the underlying constraints. Finally, a new unsupervised neural network is proposed that guarantees energy conservation through an embedded symplectic structure. The symplectic neural network is used to solve a system of energy-conserving differential equations and out-performs an unsupervised, non-symplectic neural network. △ Less

Submitted 29 January, 2020; v1 submitted 18 April, 2019; originally announced April 2019.

Comments: This is the same manuscript with version 1 (arXiv:1904.08991v1) which accidentally was replaced 16 pages, 8 figures

Imaginary-time time-dependent density functional theory and its application for robust convergence of electronic states

Authors: Cedric Flamant, Grigory Kolesov, Efstratios Manousakis, Efthimios Kaxiras

Abstract: Reliable and robust convergence to the electronic ground state within density functional theory (DFT) Kohn-Sham (KS) calculations remains a thorny issue in many systems of interest. In such cases, charge sloshing can delay or completely hinder the convergence. Here, we use an approach based on transforming the time-dependent DFT equations to imaginary time, followed by imaginary-time evolution, as… ▽ More Reliable and robust convergence to the electronic ground state within density functional theory (DFT) Kohn-Sham (KS) calculations remains a thorny issue in many systems of interest. In such cases, charge sloshing can delay or completely hinder the convergence. Here, we use an approach based on transforming the time-dependent DFT equations to imaginary time, followed by imaginary-time evolution, as a reliable alternative to the self-consistent field (SCF) procedure for determining the KS ground state. We discuss the theoretical and technical aspects of this approach and show that the KS ground state should be expected to be the long-imaginary-time output of the evolution, independent of the exchange-correlation functional or the level of theory used to simulate the system. By maintaining self-consistency between the single-particle wavefunctions and the electronic density throughout the determination of the stationary state, our method avoids the typical difficulties encountered in SCF. To demonstrate dependability of our approach, we apply it to selected systems which struggle to converge with SCF schemes. In addition, through the van Leeuwen theorem, we affirm the physical meaningfulness of imaginary time TDDFT, justifying its use in certain topics of statistical mechanics such as in computing imaginary time path integrals. △ Less

Submitted 17 June, 2019; v1 submitted 2 March, 2019; originally announced March 2019.

Comments: 22 pages, 4 figures, 2 tables; changed order of theoretical section for clarity, text made more concise, tables/figures improved. Added additional references

Journal ref: J. Chem. Theory Comput. 2019, 15, 11, 6036-6045

Machine learning with observers predicts complex spatiotemporal behavior

Authors: G. Neofotistos, M. Mattheakis, G. D. Barmparis, J. Hizanidis, G. P. Tsironis, E. Kaxiras

Abstract: Chimeras and branching are two archetypical complex phenomena that appear in many physical systems; because of their different intrinsic dynamics, they delineate opposite non-trivial limits in the complexity of wave motion and present severe challenges in predicting chaotic and singular behavior in extended physical systems. We report on the long-term forecasting capability of Long Short-Term Memo… ▽ More Chimeras and branching are two archetypical complex phenomena that appear in many physical systems; because of their different intrinsic dynamics, they delineate opposite non-trivial limits in the complexity of wave motion and present severe challenges in predicting chaotic and singular behavior in extended physical systems. We report on the long-term forecasting capability of Long Short-Term Memory (LSTM) and reservoir computing (RC) recurrent neural networks, when they are applied to the spatiotemporal evolution of turbulent chimeras in simulated arrays of coupled superconducting quantum interference devices (SQUIDs) or lasers, and branching in the electronic flow of two-dimensional graphene with random potential. We propose a new method in which we assign one LSTM network to each system node except for "observer" nodes which provide continual "ground truth" measurements as input; we refer to this method as "Observer LSTM" (OLSTM). We demonstrate that even a small number of observers greatly improves the data-driven (model-free) long-term forecasting capability of the LSTM networks and provide the framework for a consistent comparison between the RC and LSTM methods. We find that RC requires smaller training datasets than OLSTMs, but the latter requires fewer observers. Both methods are benchmarked against Feed-Forward neural networks (FNNs), also trained to make predictions with observers (OFNNs). △ Less

Submitted 27 July, 2018; originally announced July 2018.

High performance Tunnel Field Effect Transistors based on in-plane transition metal dichalcogenide heterojunctions

Authors: Jean Choukroun, Marco Pala, Shiang Fang, Efthimios Kaxiras, Philippe Dollfus

Abstract: In-plane heterojunction tunnel field effect transistors based on monolayer transition metal dichalcogenides are studied by means of self-consistent non-equilibrium Green's functions simulations and an atomistic tight-binding Hamiltonian. We start by comparing several heterojunctions before focusing on the most promising ones, i.e WTe2-MoS2 and MoTe2-MoS2. The scalability of those devices as a func… ▽ More In-plane heterojunction tunnel field effect transistors based on monolayer transition metal dichalcogenides are studied by means of self-consistent non-equilibrium Green's functions simulations and an atomistic tight-binding Hamiltonian. We start by comparing several heterojunctions before focusing on the most promising ones, i.e WTe2-MoS2 and MoTe2-MoS2. The scalability of those devices as a function of channel length is studied, and the influence of backgate voltages on device performance is analysed. Our results indicate that, by fine-tuning the design parameters, those devices can yield extremely low sub-threshold swings (below 5mV/decade) and Ion/Ioff ratios higher than 1e8 at a supply voltage of 0.3V, making them ideal for ultra-low power consumption. △ Less

Submitted 18 July, 2018; originally announced July 2018.

Comments: 10 pages

Emergence and dynamical properties of stochastic branching in the electronic flows of disordered Dirac solids

Authors: Marios Mattheakis, G. P. Tsironis, Efthimios Kaxiras

Abstract: Graphene as well as more generally Dirac solids constitute two dimensional materials where the electronic flow is ultra relativistic. When a Dirac solid is deposited on a different substrate surface with roughness, a local random potential develops through an inhomogeneous charge impurity distribution. This external potential affects profoundly the charge flow and induces a chaotic pattern of curr… ▽ More Graphene as well as more generally Dirac solids constitute two dimensional materials where the electronic flow is ultra relativistic. When a Dirac solid is deposited on a different substrate surface with roughness, a local random potential develops through an inhomogeneous charge impurity distribution. This external potential affects profoundly the charge flow and induces a chaotic pattern of current branches that develops through focusing and defocusing effects produced by the randomness of the surface. An additional bias voltage may be used to tune the branching pattern of the charge carrier currents. We employ analytical and numerical techniques in order to investigate the onset and the statistical properties of carrier branches in Dirac solids. We find a specific scaling-type relationship that connects the physical scale for the occurrence of branches with the characteristic medium properties, such as disorder and bias field. We use numerics to test and verify the theoretical prediction as well as a perturbative approach that gives a clear indication of the regime of validity of the approach. This work is relevant to device applications and may be tested experimentally. △ Less

Submitted 24 March, 2018; v1 submitted 24 January, 2018; originally announced January 2018.

Universal behavior of dispersive Dirac cone in gradient-index plasmonic metamaterials

Authors: Matthias Maier, Marios Mattheakis, Efthimios Kaxiras, Mitchell Luskin, Dionisios Margetis

Abstract: We demonstrate analytically and numerically that the dispersive Dirac cone emulating an epsilon-near-zero (ENZ) behavior is a universal property within a family of plasmonic crystals consisting of two-dimensional (2D) metals. Our starting point is a periodic array of 2D metallic sheets embedded in an inhomogeneous and anisotropic dielectric host that allows for propagation of transverse-magnetic (… ▽ More We demonstrate analytically and numerically that the dispersive Dirac cone emulating an epsilon-near-zero (ENZ) behavior is a universal property within a family of plasmonic crystals consisting of two-dimensional (2D) metals. Our starting point is a periodic array of 2D metallic sheets embedded in an inhomogeneous and anisotropic dielectric host that allows for propagation of transverse-magnetic (TM) polarized waves. By invoking a systematic bifurcation argument for arbitrary dielectric profiles in one spatial dimension, we show how TM Bloch waves experience an effective dielectric function that averages out microscopic details of the host medium. The corresponding effective dispersion relation reduces to a Dirac cone when the conductivity of the metallic sheet and the period of the array satisfy a critical condition for ENZ behavior. Our analytical findings are in excellent agreement with numerical simulations. △ Less

Submitted 12 January, 2018; v1 submitted 6 November, 2017; originally announced November 2017.

Journal ref: Phys. Rev. B 97, 035307 (2018)

Graphene and active metamaterials: theoretical methods and physical properties

Authors: Marios Mattheakis, Giorgos P. Tsironis, Efthimios Kaxiras

Abstract: The interaction of light with matter has triggered the interest of scientists for long time. The area of plasmonics emerges in this context through the interaction of light with valence electrons in metals. The random phase approximation in the long wavelength limit is used for analytical investigation of plasmons in three-dimensional metals, in a two-dimensional electron gas and finally in the mo… ▽ More The interaction of light with matter has triggered the interest of scientists for long time. The area of plasmonics emerges in this context through the interaction of light with valence electrons in metals. The random phase approximation in the long wavelength limit is used for analytical investigation of plasmons in three-dimensional metals, in a two-dimensional electron gas and finally in the most famous two-dimensional semi-metal, namely graphene. We show that plasmons in bulk metals as well as in a two-dimensional electron gas originate from classical laws, whereas, quantum effects appear as non-local corrections. On the other hand, graphene plasmons are purely quantum modes and, thus, they would not exist in a "classical world". Furthermore, under certain circumstances, light is able to couple with plasmons on metallic surfaces, forming a surface plasmon polariton, which is very important in nanoplasmonics due to its subwavelength nature. In addition, we outline two applications that complete our theoretical investigation. Firstly, we examine how the presence of gain (active) dielectrics affects surface plasmon polariton properties and we find that there is a gain value for which the metallic losses are completely eliminated resulting to lossless plasmon propagation. Secondly, we combine monolayers of graphene in a periodic order and construct a plasmonic metamaterial that provides tunable wave propagation properties, such as epsilon-near-zero behavior, normal and negative refraction. △ Less

Submitted 6 April, 2017; originally announced April 2017.

Comments: This is a Chapter of the book: Nanoplasmonics - Fundamentals and Applications. ISBN 978-953-51-5231-6

Manipulating polarized light with a planar slab of Black Phosphorus

Authors: Constantinos A. Valagiannopoulos, Marios Mattheakis, Sharmila N. Shirodkar, Efthimios Kaxiras

Abstract: Wave polarization contains valuable information for electromagnetic signal processing and the ability to manipulate it can be extremely useful in photonic devices. In this work, we propose designs comprised of one of the emerging and interesting two-dimensional media: Black Phosphorus. Due to substantial in-plane anisotropy, a single slab of Black Phosphorus can be very efficient for manipulating… ▽ More Wave polarization contains valuable information for electromagnetic signal processing and the ability to manipulate it can be extremely useful in photonic devices. In this work, we propose designs comprised of one of the emerging and interesting two-dimensional media: Black Phosphorus. Due to substantial in-plane anisotropy, a single slab of Black Phosphorus can be very efficient for manipulating the polarization state of electromagnetic waves. We investigate Black Phosphorus slabs that filter the fields along one direction, or polarization axis rotation, or convert linear polarization to circular. These slabs can be employed as components in numerous mid-IR integrated devices. △ Less

Submitted 17 September, 2017; v1 submitted 18 March, 2017; originally announced March 2017.

Epsilon-Near-Zero behavior from plasmonic Dirac point: theory and realization using two-dimensional materials

Authors: Marios Mattheakis, Constantinos A. Valagiannopoulos, Efthimios Kaxiras

Abstract: The electromagnetic response of a two-dimensional metal embedded in a periodic array of a dielectric host can give rise to a plasmonic Dirac point that emulates Epsilon-Near-Zero (ENZ) behavior. This theoretical result is extremely sensitive to tructural features like periodicity of the dielectric medium and thickness imperfections. We propose that such a device can actually be realized by using g… ▽ More The electromagnetic response of a two-dimensional metal embedded in a periodic array of a dielectric host can give rise to a plasmonic Dirac point that emulates Epsilon-Near-Zero (ENZ) behavior. This theoretical result is extremely sensitive to tructural features like periodicity of the dielectric medium and thickness imperfections. We propose that such a device can actually be realized by using graphene as the 2D metal and materials like the layered semiconducting transition-metal dichalcogenides or hexagonal boron nitride as the dielectric host. We propose a systematic approach, in terms of design characteristics, for constructing metamaterials with linear, elliptical and hyperbolic dispersion relations which produce ENZ behavior, normal or negative diffraction. △ Less

Submitted 11 October, 2016; originally announced October 2016.

Journal ref: Phys. Rev. B 94, 201404 (2016)

Mapping Reactive Flow Patterns in Monolithic Nanoporous Catalysts

Authors: Giacomo Falcucci, Sauro Succi, Andrea Montessori, Simone Melchionna, Pietro Prestininzi, Cedric Barroo, David C. Bell, Monika M. Biener, Juergen Biener, Branko Zugic, Efthimios Kaxiras

Abstract: The development of high-efficiency porous catalyst membranes critically depends on our understanding of where the majority of the chemical conversions occur within the porous structure. This requires mapping of chemical reactions and mass transport inside the complex nano-scale architecture of porous catalyst membranes which is a multiscale problem in both the temporal and spatial domain. To addre… ▽ More The development of high-efficiency porous catalyst membranes critically depends on our understanding of where the majority of the chemical conversions occur within the porous structure. This requires mapping of chemical reactions and mass transport inside the complex nano-scale architecture of porous catalyst membranes which is a multiscale problem in both the temporal and spatial domain. To address this problem, we developed a multi-scale mass transport computational framework based on the Lattice Boltzmann Method (LBM) that allows us to account for catalytic reactions at the gas-solid interface by introducing a new boundary condition. In good agreement with experiments, the simulations reveal that most catalytic reactions occur near the gas-flow facing side of the catalyst membrane if chemical reactions are fast compared to mass transport within the porous catalyst membrane. △ Less

Submitted 11 July, 2016; originally announced July 2016.

Journal ref: Microfluidics and Nanofluidics, 20(7), 1-13, 2016

Quantized biopolymer translocation through nanopores: departure from simple scaling

Authors: Simone Melchionna, Massimo Bernaschi, Maria Fyta, Efthimios Kaxiras, Sauro Succi

Abstract: We discuss multiscale simulations of long biopolymer translocation through wide nanopores that can accommodate multiple polymer strands. The simulations provide clear evidence of folding quantization, namely, the translocation proceeds through multi-folded configurations characterized by a well-defined integer number of folds. As a consequence, the translocation time acquires a dependence on the… ▽ More We discuss multiscale simulations of long biopolymer translocation through wide nanopores that can accommodate multiple polymer strands. The simulations provide clear evidence of folding quantization, namely, the translocation proceeds through multi-folded configurations characterized by a well-defined integer number of folds. As a consequence, the translocation time acquires a dependence on the average folding number, which results in a deviation from the single-exponent power-law characterizing single-file translocation through narrow pores. The mechanism of folding quantization allows polymers above a threshold length (approximately $1,000$ persistence lengths for double-stranded DNA) to exhibit cooperative behavior and as a result to translocate noticeably faster. △ Less

Submitted 27 February, 2009; originally announced February 2009.

Comments: 6 pages, 4 figures, to appear in Phys. Rev. E (R)

Journal ref: Phys. Rev. E, 79 030901 (R) (2009)

Numerical simulation of conformational variability in biopolymer translocation through wide nanopores

Authors: Maria Fyta, Simone Melchionna, Massimo Bernaschi, Efthimios Kaxiras, Sauro Succi

Abstract: Numerical results on the translocation of long biopolymers through mid-sized and wide pores are presented. The simulations are based on a novel methodology which couples molecular motion to a mesoscopic fluid solvent. Thousands of events of long polymers (up to 8000 monomers) are monitored as they pass through nanopores. Comparison between the different pore sizes shows that wide pores can host… ▽ More Numerical results on the translocation of long biopolymers through mid-sized and wide pores are presented. The simulations are based on a novel methodology which couples molecular motion to a mesoscopic fluid solvent. Thousands of events of long polymers (up to 8000 monomers) are monitored as they pass through nanopores. Comparison between the different pore sizes shows that wide pores can host a larger number of multiple biopolymer segments, as compared to smaller pores. The simulations provide clear evidence of folding quantization in the translocation process as the biopolymers undertake multi-folded configurations, characterized by a well-defined integer number of folds. Accordingly, the translocation time is no longer represented by a single-exponent power law dependence on the length, as it is the case for single-file translocation through narrow pores. The folding quantization increases with the biopolymer length, while the rate of translocated beads at each time step is linearly correlated to the number of resident beads in the pore. Finally, analysis of the statistics over the translocation work unravels the importance of the hydrodynamic interactions in the process. △ Less

Submitted 16 January, 2009; originally announced January 2009.

Comments: 10 pages, 6 figures, to appear in J. Stat. (2009)

Hydrodynamic correlations in the translocation of biopolymer through a nanopore: theory and multiscale simulations

Authors: Maria Fyta, Simone Melchionna, Sauro Succi, Efthimios Kaxiras

Abstract: We investigate the process of biopolymer translocation through a narrow pore using a multiscale approach which explicitly accounts for the hydrodynamic interactions of the molecule with the surrounding solvent. The simulations confirm that the coupling of the correlated molecular motion to hydrodynamics results in significant acceleration of the translocation process. Based on these results, we… ▽ More We investigate the process of biopolymer translocation through a narrow pore using a multiscale approach which explicitly accounts for the hydrodynamic interactions of the molecule with the surrounding solvent. The simulations confirm that the coupling of the correlated molecular motion to hydrodynamics results in significant acceleration of the translocation process. Based on these results, we construct a phenomenological model which incorporates the statistical and dynamical features of the translocation process and predicts a power law dependence of the translocation time on the polymer length with an exponent $α$ $\approx 1.2$. The actual value of the exponent from the simulations is $α= 1.28 \pm 0.01$, which is in excellent agreement with experimental measurements of DNA translocation through a nanopore, and is not sensitive to the choice of parameters in the simulation. The mechanism behind the emergence of such a robust exponent is related to the interplay between the longitudinal and transversal dynamics of both translocated and untranslocated segments. The connection to the macroscopic picture involves separating the contributions from the blob shrinking and shifting processes, which are both essential to the translocation dynamics. △ Less

Submitted 5 September, 2008; originally announced September 2008.

Comments: 7 pages, 5 figures. to appear in Phys. Rev. E

Journal ref: Phys. Rev. E, 78 036704 (2008)

Quantized current blockade and hydrodynamic correlations in biopolymer translocation through nanopores: evidence from multiscale simulations

Authors: Massimo Bernaschi, Simone Melchionna, Sauro Succi, Maria Fyta, Efthimios Kaxiras

Abstract: We present a detailed description of biopolymer translocation through a nanopore in the presence of a solvent, using an innovative multi-scale methodology which treats the biopolymer at the microscopic scale as combined with a self-consistent mesoscopic description for the solvent fluid dynamics. We report evidence for quantized current blockade depending on the folding configuration and offer d… ▽ More We present a detailed description of biopolymer translocation through a nanopore in the presence of a solvent, using an innovative multi-scale methodology which treats the biopolymer at the microscopic scale as combined with a self-consistent mesoscopic description for the solvent fluid dynamics. We report evidence for quantized current blockade depending on the folding configuration and offer detailed information on the role of hydrodynamic correlations in speeding-up the translocation process. △ Less

Submitted 8 February, 2008; originally announced February 2008.

Comments: 6 pages, 5 figures, to appear in Nano Letters

Parallel multiscale modeling of biopolymer dynamics with hydrodynamic correlations

Authors: Maria Fyta, Jayanta Sircar, Efthimios Kaxiras, Simone Melchionna, Massimo Bernaschi, Sauro Succi

Abstract: We employ a multiscale approach to model the translocation of biopolymers through nanometer size pores. Our computational scheme combines microscopic Molecular Dynamics (MD) with a mesoscopic Lattice Boltzmann (LB) method for the solvent dynamics, explicitly taking into account the interactions of the molecule with the surrounding fluid. We describe an efficient parallel implementation of the me… ▽ More We employ a multiscale approach to model the translocation of biopolymers through nanometer size pores. Our computational scheme combines microscopic Molecular Dynamics (MD) with a mesoscopic Lattice Boltzmann (LB) method for the solvent dynamics, explicitly taking into account the interactions of the molecule with the surrounding fluid. We describe an efficient parallel implementation of the method which exhibits excellent scalability on the Blue Gene platform. We investigate both dynamical and statistical aspects of the translocation process by simulating polymers of various initial configurations and lengths. For a representative molecule size, we explore the effects of important parameters that enter in the simulation, paying particular attention to the strength of the molecule-solvent coupling and of the external electric field which drives the translocation process. Finally, we explore the connection between the generic polymers modeled in the simulation and DNA, for which interesting recent experimental results are available. △ Less

Submitted 16 November, 2007; originally announced November 2007.

Comments: 10 pages, 5 figures, 3 tables. to appear in International Journal for Multiscale Computational Engineering

Journal ref: vol. 6, issue 1 (2008)

Coupling Lattice Boltzmann with Atomistic Dynamics for the multiscale simulation of nano-biological flows

Authors: Maria Fyta, Simone Melchionna, Efthimios Kaxiras, Sauro Succi

Abstract: We describe a recent multiscale approach based on the concurrent coupling of constrained molecular dynamics for long biomolecules with a mesoscopic lattice Boltzmann treatment of solvent hydrodynamics. The multiscale approach is based on a simple scheme of exchange of space-time information between the atomistic and mesoscopic scales and is capable of describing self-consistent hydrodynamic effe… ▽ More We describe a recent multiscale approach based on the concurrent coupling of constrained molecular dynamics for long biomolecules with a mesoscopic lattice Boltzmann treatment of solvent hydrodynamics. The multiscale approach is based on a simple scheme of exchange of space-time information between the atomistic and mesoscopic scales and is capable of describing self-consistent hydrodynamic effects on molecular motion at a computational cost which scales linearly with both solute size and solvent volume. For an application of our multiscale method, we consider the much studied problem of biopolymer translocation through nanopores: we find that the method reproduces with remarkable accuracy the statistical scaling behavior of the translocation process and provides valuable insight into the cooperative aspects of biopolymer and hydrodynamic motion. △ Less

Submitted 5 October, 2007; originally announced October 2007.

Comments: 10 pages, 9 figures, to appear in Computing in Science and Engineering

Journal ref: March/April issue, p.21 (2008)

Multiscale modeling of biopolymer translocation through a nanopore

Authors: Maria Fyta, Simone Melchionna, Efthimios Kaxiras, Sauro Succi

Abstract: We employ a multiscale approach to model the translocation of biopolymers through nanometer size pores. Our computational scheme combines microscopic Langevin molecular dynamics (MD) with a mesoscopic lattice Boltzmann (LB) method for the solvent dynamics, explicitly taking into account the interactions of the molecule with the surrounding fluid. Both dynamical and statistical aspects of the tra… ▽ More We employ a multiscale approach to model the translocation of biopolymers through nanometer size pores. Our computational scheme combines microscopic Langevin molecular dynamics (MD) with a mesoscopic lattice Boltzmann (LB) method for the solvent dynamics, explicitly taking into account the interactions of the molecule with the surrounding fluid. Both dynamical and statistical aspects of the translocation process were investigated, by simulating polymers of various initial configurations and lengths. For a representative molecule size, we explore the effects of important parameters that enter in the simulation, paying particular attention to the strength of the molecule-solvent coupling and of the external electric field which drives the translocation process. Finally, we explore the connection between the generic polymers modeled in the simulation and DNA, for which interesting recent experimental results are available. △ Less

Submitted 7 March, 2007; originally announced March 2007.

Comments: 8 pages, 5 figures, to appear in Lecture Notes in Computer Science

Journal ref: Lecture Notes in Computer Science, 4487 786 (2007)

Exploring DNA translocation through a nanopore via a multiscale Lattice-Boltzmann Molecular-Dynamics methodology

Authors: Simone Melchionna, Maria Fyta, Efthimios Kaxiras, Sauro Succi

Abstract: A multiscale approach is used to simulate the translocation of DNA through a nanopore. Within this scheme, the interactions of the molecule with the surrounding fluid (solvent) are explicitly taken into account. By generating polymers of various initial configurations and lengths we map the probability distibutions of the passage times of the DNA through the nanopore. A scaling law behavior for… ▽ More A multiscale approach is used to simulate the translocation of DNA through a nanopore. Within this scheme, the interactions of the molecule with the surrounding fluid (solvent) are explicitly taken into account. By generating polymers of various initial configurations and lengths we map the probability distibutions of the passage times of the DNA through the nanopore. A scaling law behavior for the most probable of these times with respect to length is derived, and shown to exhibit an exponent that is in a good agreement with the experimental findings. The essential features of the DNA dynamics as it passes through the pore are explored. △ Less

Submitted 21 February, 2007; originally announced February 2007.

Comments: 5 pages, 4 figures. to appear in International Journal of Modern Physics C

Journal ref: Int. J. Mod. Phys. C, 18 685 (2007)

Multiscale coupling of molecular dynamics and hydrodynamics: application to DNA translocation through a nanopore

Authors: Maria G. Fyta, Simone Melchionna, Efthimios Kaxiras, Sauro Succi

Abstract: We present a multiscale approach to the modeling of polymer dynamics in the presence of a fluid solvent. The approach combines Langevin Molecular Dynamics (MD) techniques with a mesoscopic Lattice-Boltzmann (LB) method for the solvent dynamics. A unique feature of the present approach is that hydrodynamic interactions between the solute macromolecule and the aqueous solvent are handled explicitl… ▽ More We present a multiscale approach to the modeling of polymer dynamics in the presence of a fluid solvent. The approach combines Langevin Molecular Dynamics (MD) techniques with a mesoscopic Lattice-Boltzmann (LB) method for the solvent dynamics. A unique feature of the present approach is that hydrodynamic interactions between the solute macromolecule and the aqueous solvent are handled explicitly, and yet in a computationally tractable way due to the dual particle-field nature of the LB solver. The suitability of the present LB-MD multiscale approach is demonstrated for the problem of polymer fast translocation through a nanopore. We also provide an interpretation of our results in the context of DNA translocation through a nanopore, a problem that has attracted much theoretical and experimental attention recently. △ Less

Submitted 2 January, 2007; originally announced January 2007.

Comments: 18 pages, 8 figures

Journal ref: Multiscale Modeling and Simulation, 5 (2006) 1156

A lattice Boltzmann study of reactive microflows

Authors: A. Gabriellii, S. Succi, E. Kaxiras

Abstract: The role of geometrical micro-barriers on the conversion efficiency of reactive flows in narrow three-dimensional channels of millimetric size is investigated. Using a Lattice-Boltzmann-Lax-Wendroff code, we show that micro-barriers have an appreciable effect on the effective reaction efficiency of the device. If extrapolated to macroscopic scales, these effects can result in a sizeable increase… ▽ More The role of geometrical micro-barriers on the conversion efficiency of reactive flows in narrow three-dimensional channels of millimetric size is investigated. Using a Lattice-Boltzmann-Lax-Wendroff code, we show that micro-barriers have an appreciable effect on the effective reaction efficiency of the device. If extrapolated to macroscopic scales, these effects can result in a sizeable increase of the overall reaction efficiency. △ Less

Submitted 21 August, 2001; originally announced August 2001.

Comments: 5 pages, 7 figures

Chemical efficiency of reactive microflows with heterogeneus catalysis: a lattice Boltzmann study

Authors: S. Succi, A. Gabrielli, G. Smith, E. Kaxiras

Abstract: We investigate the effects of geometrical micro-irregularities on the conversion efficiency of reactive flows in narrow channels of millimetric size. Three-dimensional simulations, based upon a Lattice-Boltzmann-Lax-Wendroff code, indicate that periodic micro-barriers may have an appreciable effect on the effective reaction efficiency of the device. Once extrapolated to macroscopic scales, these… ▽ More We investigate the effects of geometrical micro-irregularities on the conversion efficiency of reactive flows in narrow channels of millimetric size. Three-dimensional simulations, based upon a Lattice-Boltzmann-Lax-Wendroff code, indicate that periodic micro-barriers may have an appreciable effect on the effective reaction efficiency of the device. Once extrapolated to macroscopic scales, these effects can result in a sizeable increase of the overall reaction efficiency. △ Less

Submitted 15 March, 2001; originally announced March 2001.

Comments: 12 pages, 12 figures