Electronic Structure of First and Second Row Atoms under Harmonic Confinement

Andres Robles-Navarro; Carlos Cardenas; Francisco Mu\~nos; Patricio Fuentealba

arxiv: 1907.10186 · v1 · pith:3UVWOKQ4new · submitted 2019-07-24 · ⚛️ physics.chem-ph · physics.atom-ph

Electronic Structure of First and Second Row Atoms under Harmonic Confinement

Andres Robles-Navarro , Patricio Fuentealba , Francisco Mu\~nos , Carlos Cardenas This is my paper

Pith reviewed 2026-05-24 17:03 UTC · model grok-4.3

classification ⚛️ physics.chem-ph physics.atom-ph

keywords harmonic confinementatomic ionization pressureelectronic structureshell contractionCCSDCASSCFelectron correlationhigh pressure chemistry

0 comments

The pith

Harmonic confinement calibrated to pressure shows first-row atoms ionize at pressures from 28 GPa for Li to 10.8 TPa for Ne, with ground-state reconfiguration in Li and Be.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper studies the electronic structure of atoms H through Ne when placed inside a harmonic potential that models high pressure. It calibrates confinement strength to physical pressure units and uses CCSD for dynamic correlation plus CASSCF for static correlation to track how shells contract, when atoms ionize, and when excited configurations cross to become the ground state. A reader would care because the results indicate both core and valence shells shrink enough that ordinary pseudopotentials may break down in dense matter, and that lithium and beryllium switch to new ground-state configurations at high pressure. The ionization pressures follow a clear periodic pattern across the row. Correlation energy rises for atoms with sparse valence shells but is overtaken by rising kinetic energy in atoms with fuller shells.

Core claim

Under harmonic confinement mapped to pressure, atoms from H to Ne lose an electron at critical pressures that increase periodically from 28 GPa for lithium to 10.8 TPa for neon. At sufficiently high pressure the ground state of lithium becomes the doublet 1s²2p configuration and that of beryllium becomes the triplet 1s²2s2p configuration. Both K and L shells contract, correlation energy changes in a manner consistent with the RPA picture, and the amount of correlation grows with fewer valence electrons while kinetic-energy growth dominates for atoms with more valence electrons.

What carries the argument

Calibration procedure that converts harmonic-confinement strength into equivalent physical pressure, used together with CCSD and CASSCF calculations on the confined atoms.

If this is right

Core and valence shell contraction is large enough that standard pseudopotentials become unreliable for solids at extreme pressures.
Pressure-induced configuration mixing in Li and Be alters their electronic ground states and may change their chemical behavior under compression.
Ionization pressure increases across the row in a periodic manner.
Correlation energy increases for atoms with few valence electrons but is dominated by kinetic-energy rise for atoms with more valence electrons.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same calibration approach could be applied to molecules or heavier atoms to predict high-pressure ionization thresholds.
The observed level crossings suggest possible changes in bonding preferences inside compressed materials containing Li or Be.
Comparison of the calculated ionization pressures against diamond-anvil-cell data on light elements would test the pressure mapping.
The periodic trend in ionization pressure may connect to trends in compressibility or metallization thresholds in the same elements.

Load-bearing premise

The chosen harmonic-potential strength can be mapped accurately onto real pressure, and the CCSD plus CASSCF calculations capture the dominant correlation and shell-contraction effects without large errors at strong confinement.

What would settle it

Direct experimental measurement of the pressure at which neon loses an electron, or spectroscopic confirmation that compressed beryllium adopts the 1s²2s2p triplet ground state above a few TPa.

read the original abstract

Atoms under pressure undergo a series of processes and modification of its electronic structure. Examples are the spontaneous ionization, stabilization of excited-state configurations that result in a level-crossing with the ground state of the free atom, and contraction of atomic shells. In this work, we do a systematic study of the effects of confinement with harmonic potential on the electronic structure of atoms from H to Ne. Dynamic and static correlation is taken into account by performing CCSD and CASSCF calculations, respectively. Because the strength of harmonic confinement cannot be translated into pressure, we envisioned a "calibration" method to transform confinement into pressure. We focused on the effect of confinement on: i) changes of electron distribution and localization within the $K$ and $L$ atomic shells, ii) confinement-induced ionization pressure, iii) level crossing of electronic states, and iv) the electron correlation energy. We found that contraction of valence and core shells are not negligible and that the use of standard pseudopotentials might be not adequate to study solids under extreme pressures. The critical pressure at which and atom ionizes follows a periodic trend, and it ranges from $28$ GPa for Li to $10.8$ TPa for Ne. In Li and a Be, pressure induces mixing of the ground state configuration with excited states. At high pressure, the ground state of Li and Be becomes a doublet and a triplet with configurations $1s^22p$ and $1s^22s2p$ respectively. The potential consequences of these changes of configuration on the chemistry of Be are discussed. Finally, the changes in the amount of electron correlation are characterized and analyzed in terms of the RPA approximation. For atoms with fewer electrons in the valence shell correlation increases, but for atoms with more electron, the increasing of kinetic energy dominates over electron correlation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper supplies concrete ionization pressures and config switches for H-Ne under confinement, but the harmonic-to-pressure calibration is the shaky step.

read the letter

The main new pieces are the ionization pressures (28 GPa for Li up to 10.8 TPa for Ne) and the high-pressure ground-state switches for Li (to 1s²2p doublet) and Be (to 1s²2s2p triplet). They also track shell contraction and correlation energy changes across the series. CCSD plus CASSCF is a reasonable choice for including both dynamic and static correlation, and running the full H-to-Ne set gives a consistent picture rather than isolated cases. The calibration that converts the harmonic parameter into physical pressure is their attempt to make the numbers usable outside the model, and the note that standard pseudopotentials may fail under strong confinement is a practical observation. The correlation trends (increasing for light valence shells, kinetic energy dominating for heavier ones) follow from the calculations and are worth checking against other work. The calibration remains the soft spot. The abstract states there is no direct mapping, so they introduce one, but without details on the reference observable, atom dependence, or checks against independent high-pressure benchmarks, the absolute GPa/TPa values and the claimed periodic trend rest on that step. The config changes themselves come straight from the energy ordering and are less sensitive to the pressure label. Basis-set and convergence information is not visible in the abstract, which matters when shells tighten. This is for the niche group working on confined atoms or extreme-pressure modeling. The calculations are standard enough and the claims specific enough that it deserves a serious referee, mainly to examine the calibration procedure and any validation they provide.

Referee Report

2 major / 2 minor

Summary. The paper computationally examines the electronic structure of atoms H–Ne under harmonic confinement using CCSD (dynamic correlation) and CASSCF (static correlation). It introduces a calibration procedure to map confinement strength ω to physical pressure P, then reports shell contractions, ionization pressures that follow a periodic trend (28 GPa for Li to 10.8 TPa for Ne), pressure-induced configuration changes (Li doublet 1s²2p, Be triplet 1s²2s2p), and trends in electron correlation energy analyzed via RPA.

Significance. If the calibration mapping is reliable and the calculations converge, the results would usefully quantify how extreme pressure alters atomic shells and ground-state configurations, with direct implications for the validity of standard pseudopotentials in high-pressure solids and for the high-pressure chemistry of Be. The direct extraction of configuration crossings from the ab initio energies is a clear strength independent of the pressure label.

major comments (2)

[Calibration procedure] Calibration procedure (following the computational methods): the mapping from harmonic confinement strength to pressure is load-bearing for all reported absolute ionization pressures and the claimed periodic trend, yet the manuscript supplies no quantitative details on the reference observables chosen for calibration, whether the mapping is atom-dependent, or any validation against independent high-pressure data or known limits (e.g., experimental or other theoretical ionization pressures).
[Computational details] Computational details section: no basis-set specifications, active-space choices for CASSCF, or convergence tests with respect to basis size or confinement strength are provided, so it is impossible to judge whether the reported shell contractions, level crossings, and correlation energies are numerically stable under strong confinement.

minor comments (2)

[Abstract] Abstract and introduction: the phrase “we envisioned a ‘calibration’ method” is informal; a concise statement of the actual procedure and its assumptions would improve clarity.
[Figures] Figure captions and text: several instances of inconsistent notation for the harmonic frequency (ω vs. ω₀) and pressure units (GPa vs. TPa) should be standardized.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed report and constructive suggestions. We address the two major comments below and will revise the manuscript accordingly to improve clarity and reproducibility.

read point-by-point responses

Referee: [Calibration procedure] Calibration procedure (following the computational methods): the mapping from harmonic confinement strength to pressure is load-bearing for all reported absolute ionization pressures and the claimed periodic trend, yet the manuscript supplies no quantitative details on the reference observables chosen for calibration, whether the mapping is atom-dependent, or any validation against independent high-pressure data or known limits (e.g., experimental or other theoretical ionization pressures).

Authors: We agree that the calibration procedure requires expanded quantitative documentation. In the revised manuscript we will specify the reference observables used to map ω to P (including the precise definition of the reference density or energy matching), state whether the mapping is performed separately for each atom or uses a universal function, and add direct comparisons to independent theoretical ionization pressures or known high-pressure limits for at least a subset of the atoms. These additions will allow readers to assess the reliability of the absolute pressures and the reported periodic trend. revision: yes
Referee: [Computational details] Computational details section: no basis-set specifications, active-space choices for CASSCF, or convergence tests with respect to basis size or confinement strength are provided, so it is impossible to judge whether the reported shell contractions, level crossings, and correlation energies are numerically stable under strong confinement.

Authors: We acknowledge that the computational details section is incomplete. The revised manuscript will report the specific basis sets (including cardinal number and augmentation), the active-space selections and orbital restrictions used in the CASSCF calculations, and the results of explicit convergence tests with respect to both basis-set size and increasing confinement strength ω. These data will confirm the numerical stability of the shell contractions, configuration crossings, and correlation-energy trends. revision: yes

Circularity Check

0 steps flagged

No significant circularity; central results are direct outputs of CCSD/CASSCF calculations under harmonic confinement.

full rationale

The derivation chain consists of standard electronic-structure computations (CCSD for dynamic correlation, CASSCF for static) on atoms confined by a harmonic potential. Ionization pressures and configuration changes are obtained by evaluating energies as a function of the confinement parameter ω and identifying crossings or thresholds. The calibration step that converts ω to physical pressure P is an auxiliary mapping introduced because no direct translation exists; it does not redefine the computed energies or make the reported pressures tautological with the input data. No self-citation is invoked as a uniqueness theorem, no ansatz is smuggled, and no fitted parameter is relabeled as an independent prediction. The periodic trend and specific values (28 GPa for Li, 10.8 TPa for Ne) therefore emerge from the quantum-chemical results rather than by construction from the inputs.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Abstract alone does not enumerate explicit free parameters or axioms; the mentioned calibration step is the most likely source of fitted quantities, and the adequacy of CCSD/CASSCF under confinement is an unstated domain assumption.

free parameters (1)

calibration mapping from harmonic confinement strength to pressure
Paper states a calibration method is used to convert confinement into GPa/TPa but provides no functional form or fitting procedure in the abstract.

axioms (1)

domain assumption CCSD and CASSCF sufficiently capture dynamic and static correlation for confined atoms
Explicitly invoked to account for correlation effects in the electronic-structure calculations.

pith-pipeline@v0.9.0 · 5878 in / 1356 out tokens · 25105 ms · 2026-05-24T17:03:20.163707+00:00 · methodology

Electronic Structure of First and Second Row Atoms under Harmonic Confinement

Core claim

What carries the argument

If this is right

Where Pith is reading between the lines

Load-bearing premise

What would settle it

discussion (0)