Source author record

Miguel A. Caro

Miguel A. Caro appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.mtrl-sci cond-mat.mes-hall Machine Learning physics.comp-ph

Catalog footprint

What is connected

7works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Density-dependent sodium-storage mechanisms in hard carbon materials

Understanding the sodium-storage mechanism in hard carbon (HC) anodes is crucial for advancing sodium-ion battery (SIB) technology. However, the intrinsic complexity of HC microstructures and their interactions with sodium remain not fully elucidated. We present a multiscale methodology that integrates grand-canonical Monte Carlo (GCMC) simulations with a machine-learning interatomic potential based on the Gaussian approximation potential (GAP) framework to investigate sodium insertion mechanisms in hard carbons with different levels of porosity, achieved by simulating structural models with densities ranging from 0.7 to 1.9 g cm$^{-3}$. Structural and thermodynamic analyses reveal the interplay between pore size and accessibility and the relative contributions of adsorption, intercalation, and pore filling to the overall storage capacity. Low-density carbons favor pore-filling, achieving extremely high capacities at near-zero voltages, whereas high-density carbons primarily store sodium through adsorption and intercalation, leading to lower but more stable capacities. Intermediate-density carbons ($1.3-1.6$ g cm$^{-3}$) provide the most balanced performance, combining moderate capacity (480 and 310 mAh g$^{-1}$), safe operating voltages, and minimal volume expansion ($<10$\%). These findings establish a direct correlation between carbon density and electrochemical behavior, providing atomic-scale insight into how hard carbon morphology governs sodium-storage. The proposed framework offers a rational design principle for optimizing HC-based SIB anodes toward high energy density and long-term cycling stability.

preprint2022arXiv

Accurate computational prediction of core-electron binding energies in carbon-based materials: A machine-learning model combining density-functional theory and $\boldsymbol{GW}$

We present a quantitatively accurate machine-learning (ML) model for the computational prediction of core-electron binding energies, from which x-ray photoelectron spectroscopy (XPS) spectra can be readily obtained. Our model combines density functional theory (DFT) with $GW$ and uses kernel ridge regression for the ML predictions. We apply the new approach to materials and molecules containing carbon, hydrogen and oxygen, and obtain qualitative and quantitative agreement with experiment, resolving spectral features within 0.1 eV of reference experimental spectra. The method only requires the user to provide a structural model for the material under study to obtain an XPS prediction within seconds. Our new tool is freely available online through the XPS Prediction Server.

preprint2022arXiv

Reassignment of magic numbers for icosahedral Au clusters: 310, 564, 928 and 1426

Icosahedral Au clusters with three and four shells of atoms are found to deviate significantly from the commonly assumed Mackay structures. By introducing additional atoms in the surface shell and creating a vacancy in the center of the cluster, the calculated energy per atom can be lowered significantly, according to several different descriptions of the interatomic interaction. Analogous icosahedral structures with five and six shells of atoms are generated using the same structural motifs and are similarly found to be more stable than Mackay icosahedra. The lowest energy per atom is obtained with clusters containing 310, 564, 928 and 1426 atoms, as compared with the commonly assumed magic numbers of 309, 561, 923 and 1415. Some of the vertices in the optimized clusters have a hexagonal ring of atoms, rather than a pentagon, with the vertex atom missing. An inner shell atom in some cases moves outwards by more than an Ångström into the surface shell at the vertex site. This feature, as well as the wide distribution of nearest-neighbor distances in the surface layer, can strongly influence the catalytic properties of icosahedral clusters. The structural optimization is initially carried out using the GOUST method with atomic forces estimated with the EMT empirical potential function, but the atomic coordinates are then refined by minimization using electron density functional theory (DFT) or Gaussian approximation potential (GAP). A single energy barrier is found to separate the Mackay icosahedron from a lower energy structure where a string of atoms moves outwards in a concerted manner from the center so as to create a central vacancy while placing an additional atom in the surface shell.

preprint2020arXiv

Particle Swarm Based Hyper-Parameter Optimization for Machine Learned Interatomic Potentials

Modeling non-empirical and highly flexible interatomic potential energy surfaces (PES) using machine learning (ML) approaches is becoming popular in molecular and materials research. Training an ML-PES is typically performed in two stages: feature extraction and structure-property relationship modeling. The feature extraction stage transforms atomic positions into a symmetry-invariant mathematical representation. This representation can be fine-tuned by adjusting on a set of so-called "hyper-parameters" (HPs). Subsequently, an ML algorithm such as neural networks or Gaussian process regression (GPR) is used to model the structure-PES relationship based on another set of HPs. Choosing optimal values for the two sets of HPs is critical to ensure the high quality of the resulting ML-PES model. In this paper, we explore HP optimization strategies tailored for ML-PES generation using a custom-coded parallel particle swarm optimizer (available freely at https://github.com/suresh0807/PPSO.git). We employ the smooth overlap of atomic positions (SOAP) descriptor in combination with GPR-based Gaussian approximation potentials (GAP) and optimize HPs for four distinct systems: a toy C dimer, amorphous carbon, $α$-Fe, and small organic molecules (QM9 dataset). We propose a two-step optimization strategy in which the HPs related to the feature extraction stage are optimized first, followed by the optimization of the HPs in the training stage. This strategy is computationally more efficient than optimizing all HPs at the same time by means of significantly reducing the number of ML models needed to be trained to obtain the optimal HPs. This approach can be trivially extended to other combinations of descriptor and ML algorithm and brings us another step closer to fully automated ML-PES generation.

preprint2015arXiv

Atomistic analysis of the impact of alloy and well-width fluctuations on the electronic and optical properties of InGaN/GaN quantum wells

We present an atomistic description of the electronic and optical properties of $\text{In}_{0.25}\text{Ga}_{0.75}$N/GaN quantum wells. Our analysis accounts for fluctuations of well width, local alloy composition, strain and built-in field fluctuations as well as Coulomb effects. We find a strong hole and much weaker electron wave function localization in InGaN random alloy quantum wells. The presented calculations show that while the electron states are mainly localized by well-width fluctuations, the holes states are already localized by random alloy fluctuations. These localization effects affect significantly the quantum well optical properties,leading to strong inhomogeneous broadening of the lowest interband transition energy. Our results are compared with experimental literature data.

preprint2014arXiv

Extended scheme for the projection of material tensors of arbitrary symmetry onto a higher symmetry tensor

I propose a straightforward generalization of the projection scheme for elastic tensors introduced by Moakher and Norris [J. Elasticity 85, 215 (2006)] that takes into account also rotations. The "closest" tensor of any desired symmetry to the original tensor of lower symmetry is "closer" in this generalized scheme. The method has an important application in the context of the special quasirandom structure (SQS) method for the computational modeling of alloys, whereby the supercell's symmetry, and therefore that of the tensors representing its properties, is reduced with respect to the material's underlying symmetry. The approach allows to extract the tensor components most representative of the macroscopic symmetry of the material. Although the approach is general, in the present case I apply it to the elastic tensor and give numerical examples. Simple approximate analytical expressions for cubic materials are also provided.

preprint2013arXiv

Theory of local electric polarization and its relation to internal strain: impact on the polarization potential and electronic properties of group-III nitrides

We present a theory of local electric polarization in crystalline solids and apply it to study the case of wurtzite group-III nitrides. We show that a local value of the electric polarization, evaluated at the atomic sites, can be cast in terms of a summation over nearest-neighbor distances and Born effective charges. Within this model, the local polarization shows a direct relation to internal strain and can be expressed in terms of internal strain parameters. The predictions of the present theory show excellent agreement with a formal Berry phase calculation for random distortions of a test-case CuPt-like InGaN alloy and InGaN supercells with randomly placed cations. While the present level of theory is appropriate for highly ionic compounds, we show that a more complex model is needed for less ionic materials, in which the strain dependence of Born effective charges has to be taken into account. Moreover, we provide ab initio parameters for GaN, InN and AlN, including hybrid functional values for the piezoelectric coefficients and the spontaneous polarization, which we use to accurately implement the local theory expressions. In order to calculate the local polarization potential, we also present a point dipole method. This method overcomes several limitations related to discretization and resolution which arise when obtaining the local potential by solving Poisson's equation on an atomic grid. Finally, we perform tight-binding supercell calculations to assess the impact of the local polarization potential arising from alloy fluctuations on the electronic properties of InGaN alloys. In particular, we find that the large upward bowing with composition of the InGaN valence band edge is strongly influenced by local polarization effects. Furthermore, our analysis allows us to extract composition-dependent bowing parameters for the energy gap and valence and conduction band edges.