The equation of state of neutron matter, symmetry energy, and neutron star structure

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.Comment: 11 pages, 11 figure

Uncertainty in the determination of soil hydraulic parameters and its influence on the performance of two hydrological models of different complexity

Data of soil hydraulic properties forms often a limiting factor in unsaturated zone modelling, especially at the larger scales. Investigations for the hydraulic characterization of soils are time-consuming and costly, and the accuracy of the results obtained by the different methodologies is still debated. However, we may wonder how the uncertainty in soil hydraulic parameters relates to the uncertainty of the selected modelling approach. We performed an intensive monitoring study during the cropping season of a 10 ha maize field in Northern Italy. The data were used to: i) compare different methods for determining soil hydraulic parameters and ii) evaluate the effect of the uncertainty in these parameters on different variables (i.e. evapotranspiration, average water content in the root zone, flux at the bottom boundary of the root zone) simulated by two hydrological models of different complexity: SWAP, a widely used model of soil moisture dynamics in unsaturated soils based on Richards equation, and ALHyMUS, a conceptual model of the same dynamics based on a reservoir cascade scheme. We employed five direct and indirect methods to determine soil hydraulic parameters for each horizon of the experimental profile. Two methods were based on a parameter optimization of: a) laboratory measured retention and hydraulic conductivity data and b) field measured retention and hydraulic conductivity data. The remaining three methods were based on the application of widely used Pedo-Transfer Functions: c) Rawls and Brakensiek, d) HYPRES, and e) ROSETTA. Simulations were performed using meteorological, irrigation and crop data measured at the experimental site during the period June – October 2006. Results showed a wide range of soil hydraulic parameter values generated with the different methods, especially for the saturated hydraulic conductivity Ksat and the shape parameter a of the van Genuchten curve. This is reflected in a variability of the modeling results which is, as expected, different for each model and each variable analysed. The variability of the simulated water content in the root zone and of the bottom flux for different soil hydraulic parameter sets is found to be often larger than the difference between modeling results of the two models using the same soil hydraulic parameter set. Also we found that a good agreement in simulated soil moisture patterns may occur even if evapotranspiration and percolation fluxes are significantly different. Therefore multiple output variables should be considered to test the performances of methods and model

Quantum Monte Carlo calculations of weak transitions in A A\,= \,6--10 nuclei

Ab initio calculations of the Gamow-Teller (GT) matrix elements in the $\beta$ decays of $^6$He and $^{10}$C and electron captures in $^7$Be are carried out using both variational and Green's function Monte Carlo wave functions obtained from the Argonne $v_{18}$ two-nucleon and Illinois-7 three-nucleon interactions, and axial many-body currents derived from either meson-exchange phenomenology or chiral effective field theory. The agreement with experimental data is excellent for the electron captures in $^7$Be, while theory overestimates the $^6$He and $^{10}$C data by $\sim 2\%$ and $\sim 10\%$, respectively. We show that for these systems correlations in the nuclear wave functions are crucial to explain the data, while many-body currents increase by $\sim 2$--$3\%$ the one-body GT contributions. These findings suggest that the longstanding $g_A$-problem, i.e., the systematic overprediction ($\sim 20 \%$ in $A\le 18$ nuclei) of GT matrix elements in shell-model calculations, may be resolved, at least partially, by correlation effects.Comment: 6 pages, 2 figure

Properties of nuclei up to A=16A=16 using local chiral interactions

We report accurate quantum Monte Carlo calculations of nuclei up to $A=16$ based on local chiral two- and three-nucleon interactions up to next-to-next-to-leading order. We examine the theoretical uncertainties associated with the chiral expansion and the cutoff in the theory, as well as the associated operator choices in the three-nucleon interactions. While in light nuclei the cutoff variation and systematic uncertainties are rather small, in $^{16}$O these can be significant for large coordinate-space cutoffs. Overall, we show that chiral interactions constructed to reproduce properties of very light systems and nucleon-nucleon scattering give an excellent description of binding energies, charge radii, and form factors for all these nuclei, including open-shell systems in $A=6$ and 12.Comment: 6 pages, 4 figure

Theory of Impurity Effects on the Spin Nematic State

The effect of magnetic bond disorder in otherwise antiferro nematic ordered system is investigated. We introduced triangular-shaped ferromagnetic bond disorder in the S=1 bilinear-biquadratic model on a triangular lattice. It is shown that the coupling between the impurity magnetic moment and nonmagnetic excitation in the bulk yields single-moment anisotropy and long-range anisotropic interaction between impurity magnetic moments. This interaction can induce unconventional spin-freezing phenomena observed in triangular magnet, NiGa2S4.Comment: 19 pages, 14 figure

Confirmation of an exoplanet using the transit color signature: Kepler-418b, a blended giant planet in a multiplanet system

We announce confirmation of Kepler-418b, one of two proposed planets in this system. This is the first confirmation of an exoplanet based primarily on the transit color signature technique. We used the Kepler public data archive combined with multicolor photometry from the Gran Telescopio de Canarias and radial velocity follow-up using FIES at the Nordic Optical Telescope for confirmation. We report a confident detection of a transit color signature that can only be explained by a compact occulting body, entirely ruling out a contaminating eclipsing binary, a hierarchical triple, or a grazing eclipsing binary. Those findings are corroborated by our radial velocity measurements, which put an upper limit of ~1 Mjup on the mass of Kepler-418b. We also report that the host star is significantly blended, confirming the ~10% light contamination suspected from the crowding metric in the Kepler light curve measured by the Kepler team. We report detection of an unresolved light source that contributes an additional ~40% to the target star, which would not have been detected without multicolor photometric analysis. The resulting planet-star radius ratio is 0.110 +/- 0.0025, more than 25% more than the 0.087 measured by Kepler, leading to a radius of 1.20 +/- 0.16 Rjup instead of the 0.94 Rjup measured by the Kepler team. This is the first confirmation of an exoplanet candidate based primarily on the transit color signature, demonstrating that this technique is viable from ground for giant planets. It is particularly useful for planets with long periods such as Kepler-418b, which tend to have long transit durations. Additionally, multicolor photometric analysis of transits can reveal unknown stellar neighbors and binary companions that do not affect the classification of the transiting object but can have a very significant effect on the perceived planetary radius.Comment: accepted by Astronomy & Astrophysic

Zero Temperature Dynamics of 2D and 3D Ising Ferromagnets

We consider zero-temperature, stochastic Ising models with nearest-neighbor interactions in two and three dimensions. Using both symmetric and asymmetric initial configurations, we study the evolution of the system with time. We examine the issue of convergence of the dynamics and discuss the nature of the final state of the system. By determining a relation between the median number of spin flips per site, the probability p that a spin in the initial spin configuration takes the value +1, and lattice size, we conclude that in two and three dimensions, the system converges to a frozen (but not necessarily uniform) state when p is not equal to 1/2. Results for p=1/2 in three dimensions are consistent with the conjecture that the system does not evolve towards a fully frozen limiting state. Our simulations also uncover `striped' and `blinker' states first discussed by Spirin et al., and their statistical properties are investigated.Comment: 17 pages, 12 figure