-
A variational method for the simulation of hydrogen diffusion in metals
Authors:
Eva M. Andrés,
Ignacio Romero
Abstract:
We present a new method for the approximate solution of the strongly coupled, nonlinear stress-diffusion problem that appears when modeling hydrogen transport in metals. The most salient feature of the proposed approximation is that it is fully variational, meaning that all the discrete equations are obtained from the optimality conditions of an incremental potential, even for inelastic mechanical…
▽ More
We present a new method for the approximate solution of the strongly coupled, nonlinear stress-diffusion problem that appears when modeling hydrogen transport in metals. The most salient feature of the proposed approximation is that it is fully variational, meaning that all the discrete equations are obtained from the optimality conditions of an incremental potential, even for inelastic mechanical behavior. Like other variational methods, the proposed algorithm has remarkable properties, including the symmetry of the tangent operator, making its solution extremely efficient compared to other similar methods available in the literature.
△ Less
Submitted 5 June, 2024;
originally announced June 2024.
-
Spectroscopy of two-dimensional interacting lattice electrons using symmetry-aware neural backflow transformations
Authors:
Imelda Romero,
Jannes Nys,
Giuseppe Carleo
Abstract:
Neural networks have shown to be a powerful tool to represent ground state of quantum many-body systems, including for fermionic systems. In this work, we introduce a framework for embedding lattice symmetries in Neural Slater-Backflow-Jastrow wavefunction ansatzes, and demonstrate how our model allows us to target the ground state and low-lying excited states. To capture the Hamiltonian symmetrie…
▽ More
Neural networks have shown to be a powerful tool to represent ground state of quantum many-body systems, including for fermionic systems. In this work, we introduce a framework for embedding lattice symmetries in Neural Slater-Backflow-Jastrow wavefunction ansatzes, and demonstrate how our model allows us to target the ground state and low-lying excited states. To capture the Hamiltonian symmetries, we introduce group-equivariant backflow transformations. We study the low-energy excitation spectrum of the t-V model on a square lattice away from half-filling, and find that our symmetry-aware backflow significantly improves the ground-state energies, and yields accurate low-lying excited states for up to 10x10 lattices. We additionally compute the two-point density-correlation function and the structure factor to detect the phase transition and determine the critical point. Finally, we quantify the variational accuracy of our model using the V-score.
△ Less
Submitted 13 June, 2024;
originally announced June 2024.
-
Variational Benchmarks for Quantum Many-Body Problems
Authors:
Dian Wu,
Riccardo Rossi,
Filippo Vicentini,
Nikita Astrakhantsev,
Federico Becca,
Xiaodong Cao,
Juan Carrasquilla,
Francesco Ferrari,
Antoine Georges,
Mohamed Hibat-Allah,
Masatoshi Imada,
Andreas M. Läuchli,
Guglielmo Mazzola,
Antonio Mezzacapo,
Andrew Millis,
Javier Robledo Moreno,
Titus Neupert,
Yusuke Nomura,
Jannes Nys,
Olivier Parcollet,
Rico Pohle,
Imelda Romero,
Michael Schmid,
J. Maxwell Silvester,
Sandro Sorella
, et al. (8 additional authors not shown)
Abstract:
The continued development of novel many-body approaches to ground-state problems in physics and chemistry calls for a consistent way to assess its overall progress. Here we introduce a metric of variational accuracy, the V-score, obtained from the variational energy and its variance. We provide the most extensive curated dataset of variational calculations of many-body quantum systems to date, ide…
▽ More
The continued development of novel many-body approaches to ground-state problems in physics and chemistry calls for a consistent way to assess its overall progress. Here we introduce a metric of variational accuracy, the V-score, obtained from the variational energy and its variance. We provide the most extensive curated dataset of variational calculations of many-body quantum systems to date, identifying cases where state-of-the-art numerical approaches show limited accuracy, and novel algorithms or computational platforms, such as quantum computing, could provide improved accuracy. The V-score can be used as a metric to assess the progress of quantum variational methods towards quantum advantage for ground-state problems, especially in regimes where classical verifiability is impossible.
△ Less
Submitted 9 February, 2023;
originally announced February 2023.
-
Grain growth competition during melt pool solidification -- Comparing phase-field and cellular automaton models
Authors:
S. M. Elahi,
R. Tavakoli,
I. Romero,
D. Tourret
Abstract:
A broad range of computational models have been proposed to predict microstructure development during solidification processing but they have seldom been compared to each other on a quantitative and systematic basis. In this paper, we compare phase-field (PF) and cellular automaton (CA) simulations of polycrystalline growth in a two-dimensional melt pool under conditions relevant to additive manuf…
▽ More
A broad range of computational models have been proposed to predict microstructure development during solidification processing but they have seldom been compared to each other on a quantitative and systematic basis. In this paper, we compare phase-field (PF) and cellular automaton (CA) simulations of polycrystalline growth in a two-dimensional melt pool under conditions relevant to additive manufacturing (powder-bed fusion). We compare the resulting grain structures using local (point-by-point) measurements, as well as averaged grain orientation distributions over several simulations. We explore the effect of the CA spatial discretization level and that of the melt pool aspect ratio upon the selected grain texture. Our simulations show that detailed microscopic features related to transient growth conditions and solid-liquid interface stability (e.g. the initial planar growth stage prior to its cellular/dendritic destabilization, or the early elimination of unfavorably oriented grains due to neighbor grain sidebranching) can only be captured by PF simulations. The resulting disagreement between PF and CA predictions can only be addressed partially by a refinement of the CA grid. However, overall grain distributions averaged over the entire melt pools of several simulations seem to lead to a notably better agreement between PF and CA, with some variability with the melt pool shape and CA grid. While further research remains required, in particular to identify the appropriate selection of CA spatial discretization and its link to characteristic microstructural length scales, this research provides a useful step forward in this direction by comparing both methods quantitatively at process-relevant length and time scales.
△ Less
Submitted 30 October, 2022;
originally announced October 2022.
-
Analysis and design of bistable and thermally reversible metamaterials inspired by shape-memory alloys
Authors:
Aditya Vasudevan,
José A. Rodríguez-Martínez,
Ignacio Romero
Abstract:
In this work, we study lattice structures that exhibit a bistable behavior, i. e., they can snap from one stable state to another, and are also completely reversible, capable of reverting back to its original state through a heat treatment. We design this behavior by constructing lattice structures using networks of nonlinear springs that display tension-compression asymmetry and have different th…
▽ More
In this work, we study lattice structures that exhibit a bistable behavior, i. e., they can snap from one stable state to another, and are also completely reversible, capable of reverting back to its original state through a heat treatment. We design this behavior by constructing lattice structures using networks of nonlinear springs that display tension-compression asymmetry and have different thermal expansion coefficients. The mismatch in the thermal expansion coefficients induces residual stresses in the springs which results in the lattice structure exhibiting bistability at low temperatures and monostability at high temperatures. This behavior mimics the crystallographic phase transformations of shape memory alloys, but here artificially introduced in a structural lattice. By analyzing a representative unit cell, we quantify the effect that the stiffness and the thermal expansion coefficient of the springs have on the stability of the structural lattice. In addition, for simple 2D lattices, using the concept of universal unfoldings of singularity theory, we perform a perturbation analysis to identify the key variables of the structure where controlling defects is important, as they lead to drastic changes in the bifurcation behavior of the lattice. Finally, we verify numerically our analytical predictions in both 2D and 3D simulations using continuation techniques. The examples proposed confirm that the bistable and reversible features of the unit cell carry on to the macroscale, opening the route for the design of lattice structures for energy absorption applications that can hea} with a heat treatment.
△ Less
Submitted 11 July, 2022;
originally announced July 2022.
-
Multiscale simulation of powder-bed fusion processing of metallic alloys
Authors:
S. M. Elahi,
R. Tavakoli,
A. K. Boukellal,
T. Isensee,
I. Romero,
D. Tourret
Abstract:
We present a computational framework for the simulations of powder-bed fusion of metallic alloys, which combines: (1) CalPhaD calculations of temperature-dependent alloy properties and phase diagrams, (2) macroscale finite element (FE) thermal simulations of the material addition and fusion, and (3) microscopic phase-field (PF) simulations of solidification in the melt pool. The methodology is app…
▽ More
We present a computational framework for the simulations of powder-bed fusion of metallic alloys, which combines: (1) CalPhaD calculations of temperature-dependent alloy properties and phase diagrams, (2) macroscale finite element (FE) thermal simulations of the material addition and fusion, and (3) microscopic phase-field (PF) simulations of solidification in the melt pool. The methodology is applied to simulate the selective laser melting (SLM) of an Inconel 718 alloy using realistic processing parameters. We discuss the effect of temperature-dependent properties and the importance of accounting for different properties between the powder bed and the dense material in the macroscale thermal simulations. Using a two-dimensional longitudinal slice of the thermal field calculated via FE simulations, we perform an appropriately-converged PF solidification simulation at the scale of the entire melt pool, resulting in a calculation with over one billion grid points, yet performed on a single cluster node with eight graphics processing units (GPUs). These microscale simulations provide new insight into the grain texture selection via polycrystalline growth competition under realistic SLM conditions, with a level of detail down to individual dendrites.
△ Less
Submitted 19 March, 2022;
originally announced March 2022.
-
Development of a thermo-mechanically coupled crystal plasticity modeling framework: application to polycrystalline homogenization
Authors:
Jifeng Li,
Ignacio Romero,
Javier Segurado
Abstract:
Accurate predictions of thermo-mechanically coupled process in metals can lead to a reduction of cost and an increase of productivity in manufacturing processes such as forming. For modeling these coupled processes with the finite element method, accurate descriptions of both the mechanical and the thermal responses of the material, as well as their interaction, are needed. Conventional material m…
▽ More
Accurate predictions of thermo-mechanically coupled process in metals can lead to a reduction of cost and an increase of productivity in manufacturing processes such as forming. For modeling these coupled processes with the finite element method, accurate descriptions of both the mechanical and the thermal responses of the material, as well as their interaction, are needed. Conventional material modeling employs empirical macroscopic constitutive relations but does not account for the actual thermo-mechanical mechanisms occurring at the microscopic level. However, the consideration of the latter might be crucial to obtain accurate predictions and a complete understanding of the underlying physics. In this work we describe a fully coupled implicit thermo-mechanical framework for crystal plasticity simulations. This framework includes thermal strains, temperature dependency of the crystal behavior and heat generation by dissipation due to plastic slip and allows the use of large deformation steps thanks to the implicit integration of the governing equations. Its use within computational homogenization simulations allows to bridge the plastic deformation and temperature gradients at the macroscopic scale with the microscopic slip at the grain scale. A series of numerical examples are presented to validate the approach.
△ Less
Submitted 8 March, 2019;
originally announced March 2019.
-
Non-linear charge and energy dynamics of an adiabatically driven interacting quantum dot
Authors:
Javier I. Romero,
Pablo Roura-Bas,
Armando A. Aligia,
Liliana Arrachea
Abstract:
We formulate a general theory to study the time-dependent charge and energy transport of an adiabatically driven interacting quantum dot in contact to a reservoir for arbitrary amplitudes of the driving potential. We study within this framework the Anderson impurity model with a local ac gate voltage. We show that the exact adiabatic quantum dynamics of this system is fully determined by the behav…
▽ More
We formulate a general theory to study the time-dependent charge and energy transport of an adiabatically driven interacting quantum dot in contact to a reservoir for arbitrary amplitudes of the driving potential. We study within this framework the Anderson impurity model with a local ac gate voltage. We show that the exact adiabatic quantum dynamics of this system is fully determined by the behavior of the charge susceptibility of the frozen problem. At $T=0$, we evaluate the dynamic response functions with the numerical renormalization group (NRG). The time-resolved heat production exhibits a pronounced feature described by an instantaneous Joule law characterized by an universal resistance quantum $R_0=h/(2 e^2)$ for each spin channel. We show that this law holds in non-interacting as well as in the interacting system and also when the system is spin-polarized. In addition, in the presence of a static magnetic field, the interplay between many-body interactions and spin polarization leads to a non-trivial energy exchange between electrons with different spin components.
△ Less
Submitted 11 May, 2017; v1 submitted 2 October, 2016;
originally announced October 2016.
-
Magnetic field modulated Kondo effect in a single-magnetic-ion molecule
Authors:
Javier I. Romero,
E. Vernek,
G. B. Martins,
E. R. Mucciolo
Abstract:
We study numerically the low-temperature electronic transport properties of a single-ion magnet with uniaxial and transverse spin anisotropies. We find clear signatures of a Kondo effect caused by the presence of a transverse (zero-field) anisotropy in the molecule. Upon applying a transverse magnetic field to the single-ion magnet, we observe oscillations of the Kondo effect due to the presence o…
▽ More
We study numerically the low-temperature electronic transport properties of a single-ion magnet with uniaxial and transverse spin anisotropies. We find clear signatures of a Kondo effect caused by the presence of a transverse (zero-field) anisotropy in the molecule. Upon applying a transverse magnetic field to the single-ion magnet, we observe oscillations of the Kondo effect due to the presence of diabolical (degeneracy) points of the energy spectrum of the molecule caused by a geometrical phase interference effects similar to those observed in the quantum tunneling of multi-ion molecular nanomagnets. The field-induced lifting of the ground state degeneracy competes with the interference modulation, resulting in some cases in a suppression of the Kondo peak.
△ Less
Submitted 17 November, 2014; v1 submitted 14 April, 2014;
originally announced April 2014.
-
The single-electron transport in a three-ion magnetic molecule modulated by a transverse field
Authors:
Javier I. Romero,
Eduardo R. Mucciolo
Abstract:
We study single-electron transport in a three-ion molecule with strong uniaxial anisotropy and in the presence of a transverse magnetic field. Two magnetic ions are connected to each other through a third, nonmagnetic ion. The magnetic ions are coupled to ideal metallic leads and a back gate voltage is applied to the molecule, forming a field-effect transistor. The microscopic Hamiltonian describi…
▽ More
We study single-electron transport in a three-ion molecule with strong uniaxial anisotropy and in the presence of a transverse magnetic field. Two magnetic ions are connected to each other through a third, nonmagnetic ion. The magnetic ions are coupled to ideal metallic leads and a back gate voltage is applied to the molecule, forming a field-effect transistor. The microscopic Hamiltonian describing this system includes inter-ion hopping, on-site repulsions, and magnetic anisotropies. For a range of values of the parameters of the Hamiltonian, we obtain an energy spectrum similar to that of single-molecule magnets in the giant-spin approximation where the two states with maximum spin projection along the uniaxial anisotropy axis are well separated from other states. In addition, upon applying an external in-plane magnetic field, the energy gap between the ground and first excited states of the molecule oscillates, going to zero at certain special values of the field, in analogy to the diabolical points resulting from Berry phase interference in the giant spin model. Thus, our microscopic model provides the same phenomenological behavior expected from the giant spin model of a single-molecule magnet but with direct access to the internal structure of the molecule, thus making it more appropriate for realistic electronic transport studies. To illustrate this point, the nonlinear electronic transport in the sequential tunneling regime is evaluated for values of the field near these degeneracy points. We show that the existence of these points has a clear signature in the I-V characteristics of the molecule, most notably the modulation of excitation lines in the differential conductance.
△ Less
Submitted 30 September, 2013; v1 submitted 19 February, 2013;
originally announced February 2013.
-
Tunable Plasmon Molecules in Overlapping Nanovoids
Authors:
I. Romero,
T. V. Teperik,
F. J. Garcia de Abajo
Abstract:
Coupled and shape-tailored metallic nanoparticles are known to exhibit hybridized plasmon resonances. This Letter discuss the optical properties of a complementary system formed by overlapped nanovoid dimers buried in gold and filled with silica. This is an alternative route for plasmon engineering that benefits from vanishing radiation losses. Our analysis demonstrates the possibility of design…
▽ More
Coupled and shape-tailored metallic nanoparticles are known to exhibit hybridized plasmon resonances. This Letter discuss the optical properties of a complementary system formed by overlapped nanovoid dimers buried in gold and filled with silica. This is an alternative route for plasmon engineering that benefits from vanishing radiation losses. Our analysis demonstrates the possibility of designing artificial plasmon molecules on the basis of void plasmon hybridization, which allows fine mode tuning by varying the overlap between voids. The proposed structures could find application to both signal processing through buried optical elements and tunable-plasmon biosensing.
△ Less
Submitted 16 November, 2007;
originally announced November 2007.
-
Plasmon tunability in metallodielectric metamaterials
Authors:
S. Riikonen,
I. Romero,
F. J. Garcia de Abajo
Abstract:
The dielectric properties of metamaterials consisting of periodically arranged metallic nanoparticles of spherical shape are calculated by rigorously solving Maxwell's equations. Effective dielectric functions are obtained by comparing the reflectivity of planar surfaces limiting these materials with Fresnel's formulas for equivalent homogeneous media, showing mixing and splitting of individual-…
▽ More
The dielectric properties of metamaterials consisting of periodically arranged metallic nanoparticles of spherical shape are calculated by rigorously solving Maxwell's equations. Effective dielectric functions are obtained by comparing the reflectivity of planar surfaces limiting these materials with Fresnel's formulas for equivalent homogeneous media, showing mixing and splitting of individual-particle modes due to inter-particle interaction. Detailed results for simple cubic and fcc crystals of aluminum spheres in vacuum, silver spheres in vacuum, and silver spheres in a silicon matrix are presented. The filling fraction of the metal f is shown to determine the position of the plasmon modes of these metamaterials. Significant deviations are observed with respect to Maxwell-Garnett effective medium theory for large f, and multiple plasmons are predicted to exist in contrast to Maxwell-Garnett theory.
△ Less
Submitted 7 August, 2007;
originally announced August 2007.
-
Plasmons in nearly touching metallic nanoparticles: singular response in the limit of touching dimers
Authors:
I. Romero,
J. Aizpurua,
G. W. Bryant,
F. J. Garcia de Abajo
Abstract:
The response of gold nanoparticle dimers is studied theoretically near and beyond the limit where the particles are touching. As the particles approach each other, a dominant dipole feature is observed that is pushed into the infrared due to interparticle coupling and that is associated with a large pileup of induced charge in the interparticle gap. The redshift becomes singular as the particle…
▽ More
The response of gold nanoparticle dimers is studied theoretically near and beyond the limit where the particles are touching. As the particles approach each other, a dominant dipole feature is observed that is pushed into the infrared due to interparticle coupling and that is associated with a large pileup of induced charge in the interparticle gap. The redshift becomes singular as the particle separation decreases. The response weakens for very small separation when the coupling across the interparticle gap becomes so strong that dipolar oscillations across the pair are inhibited. Lower-wavelength, higher-order modes show a similar separation dependence in nearly touching dimers. After touching, singular behavior is observed through the emergence of a new infrared absorption peak, also accompanied by huge charge pileup at the interparticle junction, if initial interparticle-contact is made at a single point. This new mode is distinctly different from the lowest mode of the separated dimer. When the junction is made by contact between flat surfaces, charge at the junction is neutralized and mode evolution is continuous through contact. The calculated singular response explains recent experiments on metallic nanoparticle dimers and is relevant in the design of nanoparticle-based sensors and plasmon circuits.
△ Less
Submitted 7 August, 2007;
originally announced August 2007.