Researcher profile

Gediminas Gaigalas

Gediminas Gaigalas contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 published item(s)

preprint2022arXiv

Opacity of the highly ionized lanthanides and the effect on the early kilonova

We investigate the effect of the presence of lanthanides (Z = 57- 71) on the kilonova at t~hours after the neutron star merger for the first time. For this purpose, we calculate the atomic structures and the opacities for selected lanthanides: Nd (Z = 60), Sm (Z = 62), and Eu (Z = 63). We consider the ionization degree up to tenth (XI), applicable for the ejecta at t ~ a few hours after the merger, when the temperature is T ~ 10^5 K. We find that the opacities for the highly ionized lanthanides are exceptionally high, reaching k_exp~1000 cm^2/g for Eu, due to the highly dense energy levels. Using the new opacity, we perform radiative transfer simulations to show that the early light curves become fainter by a (maximum) factor of four, in comparison to lanthanide-free ejecta at t~0.1 day. However, the period at which the light curves are affected is relatively brief due to the rapid time evolution of the opacity in the outermost layer of the ejecta. We predict that for a source at a distance of ~100 Mpc, UV brightness for lanthanide-rich ejecta shows a drop to ~21-22 mag at t~0.1 day and the UV peaks around t~0.2 day with a magnitude of ~19 mag. Future detection of such a kilonova by the existing UV satellite like Swift or the upcoming UV satellite ULTRASAT will provide useful constraints on the abundance in the outer ejecta and the corresponding nucleosynthesis conditions in the neutron star mergers.

preprint2022arXiv

Tungsten vs Selenium as a potential source of kilonova nebular emission observed by Spitzer

Infrared emission lines arising from transitions between fine structure levels of heavy elements are expected to produce kilonova nebular emission. For the kilonova in GW170817, strong emission at 4.5 ${\rm μm}$ at late times was detected by the Spitzer Space Telescope but no source was detected at 3.6 ${\rm μm}$. This peculiar spectrum indicates that there exist strong line emitters around 4.5 ${\rm μm}$ and the absence of strong lines around 3.6 ${\rm μm}$. To model the spectrum we prepare a line list based on the selection rules in LS coupling from the experimentally calibrated energy levels in the NIST database. This method enables to generate the synthetic spectra with accurate line wavelengths. We find that the spectrum is sensitive to the abundance pattern whether or not the first r-process peak elements are included. In both cases, the synthetic spectra can match the observed data, leading to two possible interpretations. If the first peak elements are abundant a Se III line dominates the flux. If otherwise, W III with Os III, Rh III, and Ce IV can be the main sources. Observing nebular spectra for the future kilonovae in a wider wavelength range can provide more conclusive elemental identification.

preprint2020arXiv

$\textit{Ab initio}$ electronic factors of the $A$ and $B$ hyperfine structure constants for the $5s^25p6s \; ^{1,3}\! P^{\rm o}_{1}$ states in Sn I

Large-scale $\textit{ab initio}$ calculations of the electric field gradient, which constitutes the electronic contribution to the electric quadrupole hyperfine constant $B$, were performed for the $5s^25p6s$ $^{1,3}\!P^{\rm o}_1$ excited states of tin, using three independent computational strategies of the variational multiconfiguration Dirac-Hartree-Fock method and a fourth approach based on the configuration interaction Dirac-Fock-Sturm theory. For the $5s^25p6s$ $^{1}\!P^{\rm o}_1$ state, the final value of $B/Q =703(50)$ MHz/b differs by $0.4\%$ from the one recently used by Yordanov ${\it et~al.}$ [Communications Physics ${\bf 3}$, 107 (2020)] to extract the nuclear quadrupole moments, $Q$, for tin isotopes in the range $^{(117-131)}$Sn from collinear laser spectroscopy measurements. Efforts were made to provide a realistic theoretical uncertainty for the final $B/Q$ value of the $5s^25p6s\,^{1}\!P^{\rm o}_1$ state based on statistical principles and on correlation with the magnetic dipole hyperfine constant $A$.

preprint2020arXiv

Energy level structure and transition data of Er$^{2+}$

The main aim of this paper is to present accurate energy levels of the ground [Xe]$4f^{12}$ and first excited [Xe]$4f^{11}5d$ configurations of Er$^{2+}$. The energy level structure of the Er$^{2+}$ ion was computed using the multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction (RCI) methods, as implemented in the GRASP2018 program package. The Breit interaction, self-energy and vacuum polarization corrections were included in the RCI computations. The zero-first-order approach was used in the computations. Energy levels with the identification in $LS$ coupling for all (399) states belonging to the [Xe]$4f^{12}$ and [Xe]$4f^{11}5d$ configurations are presented. Electric dipole (E1) transition data between the levels of these two configurations are computed. The accuracy of the these data are evaluated by studying the behaviour of the transition rates as functions of the gauge parameter as well as by evaluating the cancellation factors. The core electron correlations were studied using different strategies. Root-mean-square deviations obtained in this study for states of the ground and excited configurations from the available experimental or semi-empirical data are 649 cm$^{-1}$, and 747 cm$^{-1}$, respectively.

preprint2020arXiv

Extended calculations of energy levels and transition rates for singly ionized lanthanide elements I: Pr - Gd

Lanthanide elements play important roles as an opacity source in the ejected material from neutron star mergers. Accurate and complete atomic data are necessary to evaluate the opacities and to analyze the observed data. In this paper, we perform extended, ab-initio atomic calculations from Pr II (Z=59) to Gd II (Z=64). By using multiconfiguration Dirac-Hartree-Fock and relativistic configuration-interaction methods, implemented in the general-purpose relativistic atomic structure package GRASP2K, we calculate the energy levels and transition data of electric dipole transitions. These computations are based on strategies (with small variations) of Nd II published by Gaigalas et al. (2019). Accuracy of data is evaluated by comparing computed energy levels with the NIST database or other works. For the energy levels, we obtain the average relative accuracy of 8%, 12%, 6%, 8%, and 7% for Pr II, Pm II, Sm II, Eu II, and Gd II ions, respectively as compared with the NIST data. Accuracy of energy transfer to the wavelength as 3%, 14% and 11% for Pr II, Eu II and Gd II. Our computed E1 type transition probabilities are in good agreement with experimental values presented by other authors especially for strong transitions.

preprint2020arXiv

Systematic Opacity Calculations for Kilonovae

Coalescence of neutron stars gives rise to kilonova, thermal emission powered by radioactive decays of freshly synthesized r-process nuclei. Although observational properties are largely affected by bound-bound opacities of r-process elements, available atomic data have been limited. In this paper, we study element-to-element variation of the opacities in the ejecta of neutron star mergers by performing systematic atomic structure calculations of r-process elements for the first time. We show that the distributions of energy levels tend to be higher as electron occupation increases for each electron shell due to the larger energy spacing caused by larger effects of spin-orbit and electron-electron interactions. As a result, elements with a fewer number of electrons in the outermost shells tend to give larger contributions to the bound-bound opacities. This implies that Fe is not representative for the opacities of light r-process elements. The average opacities for the mixture of r-process elements are found to be kappa ~ 20-30 cm^2 g^{-1} for the electron fraction of Ye < 0.20, kappa ~ 3-5 cm^2 g^{-1} for Ye = 0.25-0.35, and kappa ~ 1 cm^2 g^{-1} for Ye = 0.40 at T = 5,000-10,000 K, and they steeply decrease at lower temperature. We show that, even with the same abundance or Ye, the opacity in the ejecta changes with time by one order of magnitude from 1 to 10 days after the merger. Our radiative transfer simulations with the new opacity data confirm that ejecta with a high electron fraction (Ye >~ 0.25, with no lanthanide) are needed to explain the early, blue emission in GW170817/AT2017gfo while lanthanide-rich ejecta (with a mass fraction of lanthanides ~ 5 x 10^{-3}) reproduce the long-lasting near-infrared emission.

preprint2013arXiv

Cascade emission in electron beam ion trap plasma

We present investigation of the influence of cascade emission to the formation of spectra from plasma created by electron beam ion trap (EBIT) in electron trapping mode. It has been shown that cascade emission can play an important role in the formation of spectra from the EBIT plasma. Process of the cascade emission takes place when ion having cycloidal orbit leaves electron beam where coronal approximation is applicable. Thus both processes - excitation from ground or metastable levels and cascade emission - take part in the population of levels. Demonstration is based on the investigation of $W^{13+}$ spectra. The present investigation helps to resolve long-standing discrepancies; in particular, the present structure of $W^{13+}$ spectra is in good agreement with measurements on electron beam ion trap. Lines in the experimental spectra are identified as $4f^{13} 5s 5p \rightarrow 4f^{13} 5s^{2}$ and $4f^{12} 5s 5p^{2} \rightarrow 4f^{12} 5s^{2} 5p$ transitions from Dirac-Fock-Slater calculations.

preprint2013arXiv

Multiconfiguration Dirac-Hartree-Fock calculations of atomic electric dipole moments of 225^Ra, 199^Hg, and 171^Yb

The multiconfiguration Dirac-Hartree-Fock (MCDHF) method has been employed to calculate atomic electric dipole moments (EDM) of 225^Ra, 199^Hg, and 171^Yb. For the calculations of the matrix elements we extended the relativistic atomic structure package GRASP2K. The extension includes programs to evaluate matrix elements of (P, T)-odd e-N tensor-pseudotensor and pseudoscalar-scalar interactions, the atomic electric dipole interaction, the nuclear Schiff moment, and the interaction of the electron electric dipole moment with nuclear magnetic moments. The interelectronic interactions were accounted for through valence and core-valence electron correlation effects. The electron shell relaxation was included with separately optimised wave functions of opposite parities.

preprint2013arXiv

Toward Calculations with Spectroscopic Accuracy: the 2s2.2p 2Po_{3/2} - 2s2p2 4P_{5/2} Excitation Energy in Boron

No lines have been observed for transitions between the doublet and quartet levels of B I. Consequently, energy levels based on observation for the latter are obtained through extrapolation of wavelengths along the iso-electronic sequence for the 2s2.2p 2Po_{3/2} - 2s2p2 4P_{5/2} transition. In this paper, accurate theoretical excitation energies from a newly developed partitioned correlation function interaction (PCFI) method are reported for B I that include both relativistic effects in the Breit-Pauli approximation and a finite mass correction. Results are compared with extrapolated values from observed data. For B I our estimate of the excitation energy, 28959 +/- 5 cm^{-1}, is in better agreement with the values obtained by Edlén et al. (1969) than those reported by Kramida and Ryabtsev (2007). Our method is validated by applying the same procedure to the separation of these levels in C II

preprint2012arXiv

Effects of the electron correlation and Breit and hyperfine interactions on the lifetime of the 2p$^5$3s states in neutral neon

In the framework of the multiconfiguration Dirac-Hartree-Fock method, we investigate the transition properties of four excited states in the $2p^53s$ configuration of neutral neon. The electron correlation effects are taken into account systematically by using the active space approach. The effect of higher-order correlation on fine structures is shown. We also study the influence of the Breit interaction and find that it reduces the oscillator strength of the $^3P^o_1 - ^1S_0$ transition by 17%. It turns out that the inclusion of the Breit interaction is essential even for such a light atomic system. Our ab initio calculated line strengths, oscillator strengths and transition rates are compared with other theoretical values and experimental measurements. Good agreement is found except for the $^3P^o_2 - ^1S_0$ M2 transition for which discrepancies of around 15% between theories and experiments remain. In addition, the impact of hyperfine interactions on the lifetimes of the $^3P^o_0$ and $^3P^o_2$ metastable states is investigated for the $^{21}$Ne isotope (I=3/2). We find that hyperfine interactions reduce the lifetimes drastically. For the $^3P^o_0$ state the lifetime is decreased by a factor of 630.

preprint2012arXiv

Mass- and field-shift isotope parameters for the $2s - 2p$ resonance doublet of lithium-like ions

It was recently shown that dielectronic recombination measurements can be used for accurately inferring changes in the nuclear mean-square charge radii of highly-charged lithium-like neodymium [Brandau et al., Phys. Rev. Lett. 100 073201 (2008)]. To make use of this method to derive information about the nuclear charge distribution for other elements and isotopes, accurate electronic isotope shift parameters are required. In this work, we calculate and discuss the relativistic mass- and field-shift factors for the two $2s ^{2}S_{1/2} - 2p ^{2}P^{o}_{1/2,3/2}$ transitions along the lithium isoelectronic sequence. Based on the multiconfiguration Dirac-Hartree-Fock method, the electron correlation and the Breit interaction are taken into account systematically. The analysis of the isotope shifts for these two transitions along the isoelectronic sequence demonstrates the importance and competition between the mass shifts and the field shifts.

preprint2010arXiv

$M_F$-dependent Hyperfine Induced Transition Rates in an External Magnetic Field for Be-like $^{47}$Ti$^{18+}$

Hyperfine induced $2s2p~^3P_0 \rightarrow 2s^2~^1S_0$ transition rates in an external magnetic field for Be-like $^{47}$Ti were calculated based on the multiconfiguration Dirac-Fock method. It was found that the transition probability is dependent on the magnetic quantum number $M_F$ of the excited state, even in the weak field. The present investigation clarified that the difference of the hyperfine induced transition rate of Be-like Ti ions between experiment [Schippers {\sl et al.}, Phys Rev Lett {\bf 98}, (2007) 033001(4)] and theory does not result from the influence of external magnetic field.

preprint2009arXiv

MCDHF calculations of the electric dipole moment of radium induced by the nuclear Schiff moment

The multiconfiguration Dirac-Hartree-Fock theory (MCDHF) has been employed to calculate the electric dipole moment of the 7s6d 3D2 state of radium induced by the nuclear Schiff moment. The results are dominated by valence and core-valence electron correlation effects. We show that the correlation effects can be evaluated in a converged series of multiconfiguration expansions.