BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
S. Gandolfi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
This white paper was submitted to the 2022 Fundamental Symmetries, Neutrons, and Neutrinos (FSNN) Town Hall Meeting in preparation for the next NSAC Long Range Plan. We advocate to support current and future theoretical and experimental searches for physics beyond the Standard Model using nuclear $β$ decay.
The neutron skin of atomic nuclei impacts the structure of neutron-rich nuclei, the equation of state of nucleonic matter, and the size of neutron stars. Here we predict the neutron skin of selected light- and medium-mass nuclei using coupled-cluster theory and the auxiliary field diffusion Monte Carlo method with two- and three-nucleon forces from chiral effective field theory. We find a linear correlation between the neutron skin and the isospin asymmetry in agreement with the liquid-drop model and compare with data. We also extract the linear relationship that describes the difference between neutron and proton radii of mirror nuclei and quantify the effect of charge symmetry breaking terms in the nuclear Hamiltonian. Our results for the mirror-difference charge radii and binding energies per nucleon agree with existing data.
We report ab initio calculations of the S wave pairing gap in neutron matter calculated using realistic nuclear Hamiltonians that include two- and three-body interactions. We use a trial state, properly optimized to capture the essential pairing correlations, from which we extract ground state properties by means of auxiliary field diffusion Monte Carlo simulations. We extrapolate our results to the thermodynamic limit by studying the finite-size effects in the symmetry-restored projected Bardeen-Cooper-Schrieffer (PBCS) theory and compare our results to other ab initio studies done in the past. Our quantum Monte Carlo results for the pairing gap show a modest suppression with respect to the mean-field BCS values. These results can be connected to cold atom experiments, via the unitarity regime where fermionic superfluidity assumes a unified description, and they are important in the prediction of thermal properties and the cooling of neutron stars.
High-energy scattering processes, such as deep inelastic scattering (DIS) and quasielastic (QE) scattering provide a wealth of information about the structure of atomic nuclei. The remarkable discovery of the empirical linear relationship between the slope of the European Muon Collaboration (EMC) effect in DIS and the short-range-correlation (SRC) scaling factors $a_2$ in QE kinematics is naturally explained in terms of scale separation in effective field theory. This explanation has powerful consequences, allowing us to calculate and predict SRC scaling factors from ab initio low-energy nuclear theory. We present ab initio calculations of SRC scaling factors for a nucleus $A$ relative to the deuteron $a_2(A/d)$ and relative to $^3\rm He$ $a_2(A/^3\rm He)$ in light and medium-mass nuclei. Our framework further predicts that the EMC effect and SRC scaling factors have minimal or negligible isovector corrections.
We report Quantum Monte Carlo calculations of weak transitions in $A\leq 10$ nuclei, based on the Norfolk two- and three-nucleon chiral interactions, and associated one- and two-body axial currents. We find that the contribution from two-body currents is at the $2$ - $3\%$ level, with the exception of matrix elements entering the rates of $^8$Li, $^8$B, and $^8$He beta decays. These matrix elements are suppressed in impulse approximation based on the (leading order) Gamow Teller transition operator alone; two-body currents provide a $20$ - $30\%$ correction, which is, however, insufficient to bring theory in agreement with experimental data. For the other transitions, the agreement with the data is satisfactory, and the results exhibit a negligible to mild model dependence when different combinations of Norfolk interactions are utilized to construct the nuclear wave functions. We report a complete study of two-body weak transition densities which reveals the expected universal behavior of two-body currents at short distances throughout the range of $A\,$=$\,3$ to $A\,$=$\,10$ systems considered here.
We report a quantum Monte Carlo calculation of the equation of state of symmetric nuclear matter using local interactions derived from chiral effective field theory up to next-to-next-to-leading order fit to few-body observables only. The empirical saturation density and energy are well reproduced within statistical and systematic uncertainties. We have also derived the symmetry energy as a function of the density, finding good agreement with available experimentally derived constraints at saturation and twice saturation density. We find that the corresponding pressure is also in excellent agreement with recent constraints extracted from gravitational waves of the neutron-star merger GW170817 by the LIGO-Virgo detection.
We present quantum Monte Carlo calculations of light nuclei, neutron-$α$ scattering, and neutron matter using local two- and three-nucleon ($3N$) interactions derived from chiral effective field theory up to next-to-next-to-leading order (N$^2$LO). The two undetermined $3N$ low-energy couplings are fit to the $^4$He binding energy and, for the first time, to the spin-orbit splitting in the neutron-$α$ $P$-wave phase shifts. Furthermore, we investigate different choices of local $3N$-operator structures and find that chiral interactions at N$^2$LO are able to simultaneously reproduce the properties of $A=3,4,5$ systems and of neutron matter, in contrast to commonly used phenomenological $3N$ interactions.
The longitudinal and transverse electromagnetic response functions of $^{12}$C are computed in a "first-principles" Green's function Monte Carlo calculation, based on realistic two- and three-nucleon interactions and associated one- and two-body currents. We find excellent agreement between theory and experiment and, in particular, no evidence for the quenching of measured versus calculated longitudinal response. This is further corroborated by a re-analysis of the Coulomb sum rule, in which the contributions from the low-lying $J^π\,$=$\, 2^+$, $0^+_2$ (Hoyle), and $4^+$ states in $^{12}$C are accounted for explicitly in evaluating the total inelastic strength.
Chiral effective field theory (EFT) enables a systematic description of low-energy hadronic interactions with controlled theoretical uncertainties. For strongly interacting systems, quantum Monte Carlo (QMC) methods provide some of the most accurate solutions, but they require as input local potentials. We have recently constructed local chiral nucleon-nucleon (NN) interactions up to next-to-next-to-leading order (N$^2$LO). Chiral EFT naturally predicts consistent many-body forces. In this paper, we consider the leading chiral three-nucleon (3N) interactions in local form. These are included in auxiliary field diffusion Monte Carlo (AFDMC) simulations. We present results for the equation of state of neutron matter and for the energies and radii of neutron drops. In particular, we study the regulator dependence at the Hartree-Fock level and in AFDMC and find that present local regulators lead to less repulsion from 3N forces compared to the usual nonlocal regulators.
Ab initio calculations provide direct access to the properties of pure neutron systems that are challenging to study experimentally. In addition to their importance for fundamental physics, their properties are required as input for effective field theories of the strong interaction. In this work, we perform auxiliary-field diffusion Monte Carlo calculations of the ground and first excited state of two neutrons in a finite box, considering a simple contact potential as well as chiral effective field theory interactions. We compare the results against exact diagonalizations and present a detailed analysis of the finite-volume effects, whose understanding is crucial for determining observables from the calculated energies. Using the Lüscher formula, we extract the low-energy S-wave scattering parameters from ground- and excited-state energies for different box sizes.
Multi-neutron systems are crucial to understanding the physics of neutron-rich nuclei and neutron stars. Neutron drops, neutrons confined in an external field, are investigated systematically in both non-relativistic and relativistic density functional theories and with ab initio calculations. We demonstrate a strong linear correlation, which is universal in the realm of mean-field models, between the rms radii of neutron drops and the neutron skin thickness of Pb-208 and Ca-48; i.e., the difference between the neutron and proton rms radii of a nucleus. Due to its high quality, this correlation can be used to deduce the radii of neutron drops from the measured neutron skin thickness in a model-independent way, and the radii obtained for neutron drops can provide a useful constraint for realistic three neutron forces. This correlation, together with high- precision measurements of the neutron skin thicknesses of Pb-208 and Ca-48, will have an enduring impact on the understanding of multi-neutron interactions, neutron-rich nuclei, neutron stars, etc.
Precise measurement of neutrino oscillations, and hence the determination of their masses demands a quantitative understanding of neutrino-nucleus interactions. To this aim, two-body meson-exchange currents have to be accounted for along within realistic models of nuclear dynamics. We summarize our progresses towards the construction of a consistent framework, based on quantum Monte Carlo methods and on the spectral function approach, that can be exploited to accurately describe neutrino interactions with atomic nuclei over the broad kinematical region covered by neutrino experiments.
An ultracold Fermi atomic gas at unitarity presents universal properties that in the diluted limit can be well described by a contact interaction. By employing a guide function with correct boundary conditions and making simple modifications to the sampling procedure we are able to handle for the first time a true contact interaction in a quantum Monte Carlo calculation. The results are obtained with small variances. Our calculations for the Bertsch and contact parameters are in excellent agreement with published experiments. The possibility of using a more faithfully description of ultracold atomic gases can help uncover features yet unknown of the ultracold atomic gases. In addition, this work paves the way to perform quantum Monte Carlo calculations for systems interacting with contact interactions, where in many cases the description using potentials with finite effective range might not be accurate.
Ab initio calculations of the quasi-elastic electromagnetic and neutral-weak response functions of 4He and 12C are carried out for the first time. They are based on a realistic approach to nuclear dynamics, in which the strong interactions are described by two- and three-nucleon potentials and the electroweak interactions with external fields include one- and two-body terms. The Green's function Monte Carlo method is used to calculate directly the Laplace transforms of the response functions, and maximum-entropy techniques are employed to invert the resulting imaginary-time correlation functions with associated statistical errors. The theoretical results, confirmed by experiment in the electromagnetic case, show that two-body currents generate excess transverse strength from threshold to the quasi-elastic to the dip region and beyond. These findings challenge the conventional picture of quasi-elastic inclusive scattering as being largely dominated by single-nucleon knockout processes.
Quantum Monte Carlo methods have proved very valuable to study the structure and reactions of light nuclei and nucleonic matter starting from realistic nuclear interactions and currents. These ab-initio calculations reproduce many low-lying states, moments and transitions in light nuclei, and simultaneously predict many properties of light nuclei and neutron matter over a rather wide range of energy and momenta. We review the nuclear interactions and currents, and describe the continuum Quantum Monte Carlo methods used in nuclear physics. These methods are similar to those used in condensed matter and electronic structure but naturally include spin-isospin, tensor, spin-orbit, and three-body interactions. We present a variety of results including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. We also describe low-energy scattering techniques, studies of the electroweak response of nuclei relevant in electron and neutrino scattering, and the properties of dense nucleonic matter as found in neutron stars. A coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.
We present microscopic calculations of light and medium mass nuclei and the equation of state of symmetric and asymmetric nuclear matter using different nucleon-nucleon forces, including a new Argonne version that has the same spin/isospin structure as local chiral forces at next-to-next-to-leading order (N2LO). The calculations are performed using Auxiliary Field Diffusion Monte Carlo (AFDMC) combined with an improved variational wave function. We show that the AFDMC method can now be used to successfully calculate the energies of very light to medium mass nuclei as well as the energy of isospin-asymmetric nuclear matter, demonstrating microscopically the quadratic dependence of the energy on the symmetry energy.
We present details of the derivation of local chiral effective field theory interactions to next-to-next-to-leading order, and show results for nucleon-nucleon phase shifts and deuteron properties for these potentials. We then perform systematic auxiliary-field diffusion Monte Carlo calculations for neutron matter based on the developed local chiral potentials at different orders. This includes studies of the effects of the spectral-function regularization and of the local regulators. For all orders, we compare the quantum Monte Carlo results with perturbative many-body calculations and find excellent agreement for low cutoffs.
An ab initio calculation of the sum rules of the neutral weak response functions in $^{12}$C is reported, based on a realistic Hamiltonian, including two- and three-nucleon potentials, and on realistic currents, consisting of one- and two-body terms. We find that the sum rules of the response functions associated with the longitudinal and transverse components of the (space-like) neutral current are largest and that a significant portion ($\simeq 30$ %) of the calculated strength is due to two-body terms. This fact may have implications for the MiniBooNE and other neutrino quasi-elastic scattering data on nuclei.
We investigate the inhomogeneous unitary Fermi gas and use the long-wavelength properties to predict the energies of small clusters of unitary fermions trapped in harmonic potentials. The large pairing gap and scale invariance place severe restrictions on the form of the density functional. We determine the relevant universal constants needed to constrain the functional from calculations of the bulk in oscillating external potentials. Comparing with exact Quantum Monte Carlo calculations, we find that the same functional correctly predicts the lack of shell closures for small clusters of fermions trapped in harmonic wells as well as their absolute energies. A rapid convergence to the bulk limit in three dimensions, where the surface to volume ratio is quite large, is demonstrated. The resulting functional can be tested experimentally, and is a key ingredient in predicting possible polarized superfluid phases and the properties of the unitary Fermi gas in optical lattices.
We present the first Green's function Monte Carlo calculations of light nuclei with nuclear interactions derived from chiral effective field theory up to next-to-next-to-leading order. Up to this order, the interactions can be constructed in a local form and are therefore amenable to quantum Monte Carlo calculations. We demonstrate a systematic improvement with each order for the binding energies of $A=3$ and $A=4$ systems. We also carry out the first few-body tests to study perturbative expansions of chiral potentials at different orders, finding that higher-order corrections are more perturbative for softer interactions. Our results confirm the necessity of a three-body force for correct reproduction of experimental binding energies and radii, and pave the way for studying few- and many-nucleon systems using quantum Monte Carlo methods with chiral interactions.
An "ab initio" calculation of the Carbon-12 elastic form factor, and sum rules of longitudinal and transverse response functions measured in inclusive (e,e') scattering, is reported, based on realistic nuclear potentials and electromagnetic currents. The longitudinal elastic form factor and sum rule are found to be in satisfactory agreement with available experimental data. A direct comparison between theory and experiment is difficult for the transverse sum rule. However, it is shown that the calculated one has large contributions from two-body currents, indicating that these mechanisms lead to a significant enhancement of the quasi-elastic transverse response. This fact may have implications for the anomaly observed in recent neutrino quasi-elastic charge-changing scattering data off Carbon-12.
Background: The equation of state (EoS) of nucleonic matter is central for the understanding of bulk nuclear properties, the physics of neutron star crusts, and the energy release in supernova explosions. Purpose: This work presents coupled-cluster calculations of infinite nucleonic matter using modern interactions from chiral effective field theory (EFT). It assesses the role of correlations beyond particle-particle and hole-hole ladders, and the role of three-nucleon-forces (3NFs) in nuclear matter calculations with chiral interactions. Methods: This work employs the optimized nucleon-nucleon NN potential NNLOopt at next-to-next-to leading-order, and presents coupled-cluster computations of the EoS for symmetric nuclear matter and neutron matter. The coupled-cluster method employs up to selected triples clusters and the single-particle space consists of a momentum-space lattice. We compare our results with benchmark calculations and control finite-size effects and shell oscillations via twist-averaged boundary conditions. Results: We provide several benchmarks to validate the formalism and show that our results exhibit a good convergence toward the thermodynamic limit. While neutron matter with interactions from chiral EFT is perturbative, symmetric nuclear matter requires nonperturbative approaches. Correlations beyond the standard particle-particle ladder approximation yield non-negligible contributions. The saturation point of symmetric nuclear matter is sensitive to the employed 3NFs and the employed regularization scheme. 3NFs with nonlocal cutoffs exhibit a considerably improved convergence than their local cousins. We are unable to find values for the parameters of the short-range part of the local 3NF that simultaneously yield acceptable values for the saturation point in symmetric nuclear matter and the binding energies of light nuclei.
Background: The calculation of the hyperon binding energy in hypernuclei is crucial to understanding the interaction between hyperons and nucleons. Purpose: We assess the relative importance of two- and three-body hyperon-nucleon force by studying the effect of the hyperon-nucleon-nucleon interaction in closed shell Λ-hypernuclei from A=5 to 91. Methods: The Λ-binding energy has been calculated using the auxiliary field diffusion Monte Carlo method for the first time, to study light and heavy hypernuclei within the same model. Results: Our results show that including a three-body component in the hyperon-nucleon interaction leads to a saturation of the Λ-binding energy remarkably close to the experimental data. In contrast, the two-body force alone gives an unphysical limit for the binding energy. Conclusions: The repulsive contribution of the three-body hyperon-nucleon-nucleon force is essential to reproduce, even qualitatively, the binding energy of the hypernuclei in the mass range considered.
We calculate properties of neutron drops in external potentials using both quantum Monte Carlo and no-core full configuration techniques. The properties of the external wells are varied to examine different density profiles. We compare neutron drop results given by a selection of nuclear Hamiltonians, including realistic two-body interactions as well as several three-body forces. We compute a range of properties for the neutron drops: ground-state energies, spin-orbit splittings, excitation energies, radial densities and rms radii. We compare the equations of state for neutron matter for several of these Hamiltonians. Our results can be used as benchmarks to test other many-body techniques, and to constrain properties of energy-density functionals.
We present the first quantum Monte Carlo (QMC) calculations with chiral effective field theory (EFT) interactions. To achieve this, we remove all sources of nonlocality, which hamper the inclusion in QMC calculations, in nuclear forces to next-to-next-to-leading order. We perform auxiliary-field diffusion Monte Carlo (AFDMC) calculations for the neutron matter energy up to saturation density based on local leading-order, next-to-leading order, and next-to-next-to-leading order nucleon-nucleon interactions. Our results exhibit a systematic order-by-order convergence in chiral EFT and provide nonperturbative benchmarks with theoretical uncertainties. For the softer interactions, perturbative calculations are in excellent agreement with the AFDMC results. This work paves the way for QMC calculations with systematic chiral EFT interactions for nuclei and nuclear matter, for testing the perturbativeness of different orders, and allows for matching to lattice QCD results by varying the pion mass.
We study the spin response of cold dense neutron matter in the limit of zero momentum transfer, and show that the frequency dependence of the long-wavelength spin response is well constrained by sum-rules and the asymptotic behavior of the two-particle response at high frequency. The sum-rules are calculated using Auxiliary Field Diffusion Monte Carlo technique and the high frequency two-particle response is calculated for several nucleon-nucleon potentials. At nuclear saturation density, the sum-rules suggest that the strength of the spin response peaks at $ω\simeq$ 40--60 MeV, decays rapidly for $ω\geq $100 MeV, and has a sizable strength below 40 MeV. This strength at relatively low energy may lead to enhanced neutrino production rates in dense neutron-rich matter at temperatures of relevance to core-collapse supernova.
We review the calculation of the equation of state of pure neutron matter using quantum Monte Carlo (QMC) methods. QMC algorithms permit the study of many-body nuclear systems using realistic two- and three-body forces in a nonperturbative framework. We present the results for the equation of state of neutron matter, and focus on the role of three-neutron forces at supranuclear density. We discuss the correlation between the symmetry energy, the neutron star radius and the symmetry energy. We also combine QMC and theoretical models of the three-nucleon interactions, and recent neutron star observations to constrain the value of the symmetry energy and its density dependence.
We investigate the role of two- and three-body $Λ$-nucleon forces by computing the ground state of a few $Λ$-hypernuclei with the Auxiliary Field Diffusion Monte Carlo algorithm. Calculations have been performed for masses up to A=41, including some open-shell hypernuclei. The results show that the use of a bare hyperon-nucleon force fitted on the available scattering data yields a consistent overestimate of the $Λ$-separation energy $B_Λ$. The inclusion of a hyperon-nucleon-nucleon interaction systematically reduces $B_Λ$, leading to a qualitatively good agreement with experimental data over the range of masses investigated.
Using a phenomenological form of the equation of state of neutron matter near the saturation density which has been previously demonstrated to be a good characterization of quantum Monte Carlo simulations, we show that currently available neutron star mass and radius measurements provide a significant constraint on the equation of state of neutron matter. At higher densities we model the equation of state using polytropes and a quark matter model, and we show that our results do not change strongly upon variation of the lower boundary density where these polytropes begin. Neutron star observations offer an important constraint on a coefficient which is directly connected to the strength of the three-body force in neutron matter, and thus some theoretical models of the three-body may be ruled out by currently available astrophysical data. In addition, we obtain an estimate of the symmetry energy of nuclear matter and its slope that can be directly compared to the experiment and other theoretical calculations.
The symmetry energy contribution to the nuclear Equation of State (EoS) impacts various phenomena in nuclear astrophysics, nuclear structure, and nuclear reactions. Its determination is a key objective of contemporary nuclear physics with consequences for the understanding of dense matter within neutron stars. We examine the results of laboratory experiments that have provided initial constraints on the nuclear symmetry energy and its density dependence at and somewhat below normal nuclear matter density. Some of these constraints have been derived from properties of nuclei. Others have been derived from the nuclear response to electroweak and hadronic probes. We also examine the most frequently used theoretical models that predict the symmetry energy and its slope. By comparing existing constraints on the symmetry pressure to theories, we demonstrate how the contribution of the three-body force, an essential ingredient in neutron matter models, can be determined.
Background: Neutrino-nucleus quasi-elastic scattering is crucial to interpret the neutrino oscillation results in long baseline neutrino experiments. There are rather large uncertainties in the cross section, due to insufficient knowledge on the role of two-body weak currents. Purpose: Determine the role of two-body weak currents in neutrino-deuteron quasi-elastic scattering up to GeV energies. Methods: Calculate cross sections for inclusive neutrino scattering off deuteron induced by neutral and charge-changing weak currents, from threshold up to GeV energies, using the Argonne $v_{18}$ potential and consistent nuclear electroweak currents with one- and two-body terms. Results: Two-body contributions are found to be small, and increase the cross sections obtained with one-body currents by less than 10% over the whole range of energies. Total cross sections obtained by describing the final two-nucleon states with plane waves differ negligibly, for neutrino energies $\gtrsim 500$ MeV, from those in which interaction effects in these states are fully accounted for. The sensitivity of the calculated cross sections to different models for the two-nucleon potential and/or two-body terms in the weak current is found to be weak. Comparing cross sections to those obtained in a naive model in which the deuteron is taken to consist of a free proton and neutron at rest, nuclear structure effects are illustrated to be non-negligible. Conclusion: Contributions of two-body currents in neutrino-deuteron quasi-elastic scattering up to GeV are found to be smaller than 10%. Finally, it should be stressed that the results reported in this work do not include pion production channels.
We calculate the equation of state of neutron matter with realistic two- and three-nucleon interactions using quantum Monte Carlo techniques, and illustrate that the short-range three-neutron interaction determines the correlation between neutron matter energy at nuclear saturation density and higher densities relevant to neutron stars. Our model also makes an experimentally testable prediction for the correlation between the nuclear symmetry energy and its density dependence -- determined solely by the strength of the short-range terms in the three neutron force. The same force provides a significant constraint on the maximum mass and radius of neutron stars.
The contact parameter in unitary Fermi Gases governs the short-distance, high-momentum, and high-energy properties of the system. We perform accurate quantum Monte Carlo calculations with highly optimized trial functions to precisely determine this parameter at T=0, demonstrate its universal application to a variety of observables, and determine the regions of momentum and energy over which the leading short-range behavior is dominant. We derive Tan's expressions for the contact parameter using just the short-range behavior of the ground-state many-body wave function, and use this behavior to calculate the two-body distribution function, one-body density matrix, and the momentum distribution of unitary Fermi gases; providing a precise value of the contact parameter that can be compared to experiments.
The properties of inhomogeneous neutron matter are crucial to the physics of neutron-rich nuclei and the crust of neutron stars. Advances in computational techniques now allow us to accurately determine the binding energies and densities of many neutrons interacting via realistic microscopic interactions and confined in external fields. We perform calculations for different external fields and across several shells to place important constraints on inhomogeneous neutron matter, and hence the large isospin limit of the nuclear energy density functionals that are used to predict properties of heavy nuclei and neutron star crusts. We find important differences between microscopic calculations and current density functionals; in particular the isovector gradient terms are significantly more repulsive than in traditional models, and the spin-orbit and pairing forces are comparatively weaker.
We examine the physics of two, three, and four heavy fermions interacting with a single light fermion via short-range interactions. Four-particle bosonic Efimov states have proven important experimentally and also been the subject of significant theoretical effort. Similar fermionic systems are just now being investigated. We find that with some simple interactions the four- and five-particle states collapse to the interaction range at smaller mass ratios than the three-body state, and also before larger clusters can collapse. These states and their excitations can be studied in cold atom experiments, providing unique insights into the role of few-body systems in many-body physics.
We present results for neutron star models constructed with a new equation of state for nuclear matter at zero temperature. The ground state is computed using the Auxiliary Field Diffusion Monte Carlo (AFDMC) technique, with nucleons interacting via a semi-phenomenological Hamiltonian including a realistic two-body interaction. The effect of many-body forces is included by means of additional density-dependent terms in the Hamiltonian. In this letter we compare the properties of the resulting neutron-star models with those obtained using other nuclear Hamiltonians, focusing on the relations between mass and radius, and between the gravitational mass and the baryon number.
We introduce a variant to the Diffusion Monte Carlo algorithm that can be employed to study the effects of the Rashba interaction in a many electron systems. Because of the spin--orbit nature of Rashba interaction a standard algorithm cannot be applied and therefore a specific imaginary time spin dependent propagator has been developed and implemented following previous work developed in the framework of nuclear physics. We computed the ground state energy of the 2D electron gas at different densities for several values of the Rashba interaction strength as a function of "Rashba spin states" polarization. Comparison is given with analytically known Hartree-Fock results and for the system in absence of Coulomb interaction.