Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fitting”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

Estimating the color lifetime of energetic quarks

Using a simple geometric framework with a realistic nuclear density distribution, we fit published HERMES data to determine fundamental properties of hadronization using the nuclear medium as a spatial analyzer. Our approach uses a fit to the transverse momentum broadening observable and the hadronic multiplicity ratio; the simultaneous fit to two different observables strongly constrains the outcome. Using the known sizes of the target nuclei, we extract the color lifetime, finding a z h -dependent range of values from 2 to 8 fm/c for these data. We also extract estimates for the qˆ transport coefficient characterizing the strength of the interaction between the quark and the cold nuclear medium, finding an average value of 0.035±0.011 GeV 2 /fm. With a three-parameter model we obtain satisfactory fits to the data with a goodness-of-fit parameter χ 2 /dof of 1.1 or less. In a secondary fit of the results from that model we independently find a value for the Lund String Model string tension of 1.00±0.05 GeV/fm. We evaluated the sensitivity for extracting quark energy loss and effective in-medium hadronic cross sections using four-parameter variants of the model, finding large uncertainties in both cases. Our results suggest that hadronic interaction of forming hadrons in the nuclear medium is the primary dynamical cause of meson attenuation in the HERMES data, with quark energy loss playing a more minor role.

Electron-ion collider↗

Interpreting the Operando XANES of Surface-Supported Subnanometer Clusters: When Fluxionality, Oxidation State, and Size Effect Fight.

X-ray absorption near edge structure (XANES) spectroscopy is widely used for operando catalyst characterization. We show that, for highly fluxional supported nanoclusters, the customary extraction of the oxidation state of the metal from the XANES data by fitting to the bulk standards is highly questionable. The XANES signatures as well as the apparent oxidation state for such clusters arise from a complex combination of many factors, and not only from the chemical composition in reaction conditions (e.g., oxygen content in oxidizing atmosphere). The thermally accessible isomerization and population of several structurally distinct cluster forms, cluster-support interaction, and intrinsic size effects all impact the metal oxidation state and XANES signal. We demonstrate this on copper oxide clusters with different compositions, Cu4Ox(x = 2-5) and Cu5Oy (y = 3, 5), deposited on amorphous alumina and ultrananocrystalline diamond, for which we computed the XANES spectra and compare the results to the experiment. We show in addition that fitting the experimental spectrum to calculated spectra of supported clusters can, in contrast, provide good agreement and insight into the spectrum-composition-structure relation. Experimental XANES interpreted using the proposed fitting scheme shows the partial reduction of Cu oxide clusters at rising temperatures, and pinpoints the specific stoichiometries that dominate in the ensemble of cluster states as the temperature changes.

Zandkarimi, Borna↗

WUS256: An Adjoint Waveform Tomography Model of the Crust and Upper Mantle of the Western United States for Improved Waveform Simulations

Abstract We report a new model (WUS256) of radially anisotropic seismic wavespeeds of the crust and upper mantle of the western United States (WUS) obtained from adjoint waveform tomography for the purpose of improving synthetic waveform fits to observed data. WUS256 is based on inversion of over 94,000 waveforms from 72 earthquakes recorded by nearly 3,400 stations. We started with the SPiRaL global model (Simmons et al., 2021, https://doi.org/10.1093/gji/ggab277 ) and waveforms in the period band of 50–120 s. We followed a conservative multiscale inversion approach with eight stages and 256 total inversion iterations which enabled monotonic misfit reduction to 20‐s minimum‐period waves. WUS256 relied on time‐frequency (TF) phase misfits and a trust region limited memory Broyden–Fletcher–Goldfarb–Shanno (L‐BFGS) optimization. Hessian‐vector products were used to qualitatively assess model resolution. Results indicate that WUS256 has good coverage of the continental regions to depths of about 150 km and is able to resolve features on lateral scales of about 200 km. We quantify waveform fits by the reduction in TF and normalized amplitude difference misfits between WUS256 and the SPiRaL starting model. WUS256 significantly improves waveform fits with misfit reduction 64% for both inversion and validation data sets compared to the SPiRaL starting model and shows even better fits compared to other models. Waveform fits illustrate that WUS256 reproduces body‐waves, fundamental mode surface waves as well as late arriving dispersed and/or scattered short period surface waves. The improvement in waveform fit indicates that WUS256 can be used to reproduce path effects on regional complete waveforms and moment tensor inversions.

58 GEOSCIENCES↗

Low energy neutron light output characterization of EJ301D and deuterated stilbene with a comparison of light output characterization methods

The neutron-induced light yield of a 2.54 cm diameter by 2.54 cm long right circular cylinder of EJ301D and a (5.08 cm)3 custom made cube of deuterated trans-stilbene-d12 (d-stilbene) were measured over incident neutron energies from 300 keV to 2.2 MeV and 200 keV to 2.4 MeV, respectively. The measurements were performed using a time-of-flight experiment with a Cf-252 source and an approximately 1.5 m flight path. We compare three light output spectrum full energy deposition edge estimation methods: (1) simulating the neutron energy spectrum edge and fitting it to the light output spectrum, (2) using the inflection point of the light output spectrum edge (derivative method, a.k.a. Kornilov’s method), and (3) using an empirical model fit to the edge of the light output spectrum. Both the derivative and equation fit methods do not account for physical processes such as multiple neutron scattering in the detectors. They instead rely on assumptions about the linear shape continuum shape of the light output spectrum and the direct correlation between the location of the spectrum’s inflection point and maximum energy deposition. These assumptions were found to introduce bias into those methods when tested against simulated spectra with known edge locations. When tested against measured spectra the derivative method was found to differ from the simulation fit by greater than 30% at low energies with large discontinuities for adjacent data points above 800 keV incident neutron energy. The empirical equation fitting method was found to also exhibit bias of a similar magnitude, but with significantly more continuous behavior, especially with the lower count data of the smaller volumed EJ301D scintillator. Experimental light output yield for this neutron energy range is reported using the simulated spectrum fitting method because it includes physics neglected by the other methods, and did not exhibit the bias observed in the other methods

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Enabling complete multichannel nonadiabatic dynamics: A global representation of the two-channel coupled, 1,2 1 A and 1 3 A states of NH 3 using neural networks

Global coupled three-state two-channel potential energy and property/interaction (dipole and spin–orbit coupling) surfaces for the dissociation of NH 3 (Ã) into NH + H 2 and NH 2 + H are reported. The permutational invariant polynomial-neural network approach is used to simultaneously fit and diabatize the electronic Hamiltonian by fitting the energies, energy gradients, and derivative couplings of the two coupled lowest-lying singlet states as well as fitting the energy and energy gradients of the lowest-lying triplet state. The key issue in fitting property matrix elements in the diabatic basis is that the diabatic surfaces must be smooth, that is, the diabatization must remove spikes in the original adiabatic property surfaces attributable to the switch of electronic wavefunctions at the conical intersection seam. Here, we employ the fit potential energy matrix to transform properties in the adiabatic representation to a quasi-diabatic representation and remove the discontinuity near the conical intersection seam. The property matrix elements can then be fit with smooth neural network functions. The coupled potential energy surfaces along with the dipole and spin–orbit coupling surfaces will enable more accurate and complete treatment of optical transitions, as well as nonadiabatic internal conversion and intersystem crossing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Potential energy surfaces for high-energy N + O 2 collisions

Potential energy surfaces for high-energy collisions between an oxygen molecule and a nitrogen atom are useful for modeling chemical dynamics in shock waves. In the present work, we present doublet, quartet, and sextet potential energy surfaces that are suitable for studying collisions of O 2 ( 3 Σ$^{–}_{g}$) with N( 4 S) in the electronically adiabatic approximation. Two sets of surfaces are developed, one using neural networks (NNs) with permutationally invariant polynomials (PIPs) and one with the least-squares many-body (MB) method, where a two-body part is an accurate diatomic potential and the three-body part is expressed with connected PIPs in mixed-exponential-Gaussian bond order variables (MEGs). We find, using the same dataset for both fits, that the fitting performance of the PIP-NN method is significantly better than that of the MB-PIP-MEG method, even though the MB-PIP-MEG fit uses a higher-order PIP than those used in previous MB-PIP-MEG fits of related systems (such as N 4 and N 2 O 2 ). However, the evaluation of the PIP-NN fit in trajectory calculations requires about 5 times more computer time than is required for the MB-PIP-MEG fit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Atacama Cosmology Telescope: delensed power spectra and parameters

We present ΛCDM cosmological parameter constraints obtained from delensed microwave background power spectra. Lensing maps from a subset of DR4 data from the Atacama Cosmology Telescope (ACT) are used to undo the lensing effect in ACT spectra observed at 150 and 98 GHz. At 150 GHz, we remove the lensing distortion with an effective efficiency of 30% (TT), 30% (EE), 26% (TE) and 20% (BB); this results in detections of the delensing effect at 8.7σ (TT), 5.1σ (EE), 2.6σ (TE), and 2.4σ (BB) significance. The combination of 150 and 98 GHz TT, EE, and TE delensed spectra is well fit by a standard ΛCDM model. We also measure the shift in best-fit parameters when fitting delensed versus lensed spectra; while this shift does not inform our ability to measure cosmological parameters, it does provide a three-way consistency check among the lensing inferred from the best-fit parameters, the lensing in the CMB power spectrum, and the reconstructed lensing map. This shift is predicted to be zero when fitting with the correct model since both lensed and delensed spectra originate from the same region of sky. Furthermore, fitting with a ΛCDM model and marginalizing over foregrounds, we find that the shift in cosmological parameters is consistent with zero. Our results show that gravitational lensing of the microwave background is internally consistent within the framework of the standard cosmological model.

79 ASTRONOMY AND ASTROPHYSICS↗

Inverse modeling of circular lattices via orbit response measurements in the presence of degeneracy

The number and location of beam position monitors (BPMs) and steerers with respect to the quadrupoles in a circular lattice can lead to degeneracy in the context of fitting linear optics and extracting lattice information from measured closed orbits. Furthermore, the measurement uncertainties due to the imperfection of BPMs and steerers can be propagated by the fitting process in ways that prohibit the successful extraction of discrepancies between lattice elements in the real machine and their description in the corresponding model. We systematically studied the influence of the placement of BPMs and steerers on the reconstruction of linear optics and corresponding lattice information. The derivative of orbit response coefficients with respect to the quadrupole strengths, the Jacobian, is derived as an analytical formula. This analytical version of the Jacobian is used to further derive the theoretical limitations of fitting linear optics from closed orbits in terms of the placement of BPMs and steerers. It is further demonstrated that when evaluating the Jacobian during the fitting procedure, the analytical version can be used in place of the conventional finite-difference computation. This allows for greatly improved efficiency when computing the Jacobian during each iteration of the fitting procedure. The approach is tested with large-scale simulations and the findings are verified by measurement data taken on SIS18 synchrotron at GSI Helmholtz Centre for Heavy Ion Research. The presented methods are of general nature and can be applied to other accelerator lattices as well. The fitting procedure by using the analytical Jacobian is tested in conjunction with various methods for mitigating quasidegeneracy and the results agree with those obtained by using the conventional Jacobian via finite-difference approximation.

47 OTHER INSTRUMENTATION↗

Hidden conformal symmetry from the lattice

We analyze newly expanded and refined data from lattice studies of an SU(3) gauge theory with eight Dirac fermions in the fundamental representation. We focus on the light composite states emerging from these studies, consisting of a set of pseudoscalars and a single light scalar. We first consider the view that this theory is just outside the conformal window. In this case, the pseudoscalars arise from spontaneous breaking of chiral symmetry. Identifying the scalar in this case as an approximate dilaton, we fit the lattice data to a dilaton effective field theory, finding that it yields a good fit even at lowest order. For comparison, we then consider the possibility that the theory is inside the conformal window. The fermion mass provides a deformation, triggering confinement. We employ simple scaling laws to fit the lattice data, and find that it is of lesser quality.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Symmetry-mode analysis for local structure investigations using pair distribution function data

Symmetry-adapted distortion modes provide a natural way of describing distorted structures derived from higher-symmetry parent phases. Structural refinements using symmetry-mode amplitudes as fit variables have been used for at least ten years in Rietveld refinements of the average crystal structure from diffraction data; more recently, this approach has also been used for investigations of the local structure using real-space pair distribution function (PDF) data. Here, the value of performing symmetry-mode fits to PDF data is further demonstrated through the successful application of this method to two topical materials: TiSe2, where a subtle but long-range structural distortion driven by the formation of a charge-density wave is detected, and MnTe, where a large but highly localized structural distortion is characterized in terms of symmetry-lowering displacements of the Te atoms. Here, the analysis is performed using fully open-source code within the DiffPy framework via two packages developed for this work: isopydistort, which provides a scriptable interface to the ISODISTORT web application for group theoretical calculations, and isopytools, which converts the ISODISTORT output into a DiffPy-compatible format for subsequent fitting and analysis. These developments expand the potential impact of symmetry-adapted PDF analysis by enabling high-throughput analysis and removing the need for any commercial software.

36 MATERIALS SCIENCE↗

Advanced EXAFS analysis techniques applied to the L -edges of the lanthanide oxides

The unique properties of the lanthanide (Ln) elements make them critical components of modern technologies, such as lasers, anti-corrosive films and catalysts. Thus, there is significant interest in establishing structure–property relationships for Ln-containing materials to advance these technologies. Extended X-ray absorption fine structure (EXAFS) is an excellent technique for this task considering its ability to determine the average local structure around the Ln atoms for both crystalline and amorphous materials. However, the limited availability of EXAFS reference spectra of the Ln oxides and challenges in the EXAFS analysis have hindered the application of this technique to these elements. The challenges include the limited k-range available for the analysis due to the superposition of L-edges on the EXAFS, multielectron excitations (MEEs) creating erroneous peaks in the EXAFS and the presence of inequivalent absorption sites. Herein, we removed MEEs to model the local atomic environment more accurately for light Ln oxides. Further, we investigated the use of cubic and non-cubic lattice expansion to minimize the fitting parameters needed and connect the fitting parameters to physically meaningful crystal parameters. The cubic expansion reduced the number of fitting parameters but resulted in a statistically worse fit. The non-cubic expansion resulted in a similar quality fit and showed non-isotropic expansion in the crystal lattice of Nd 2 O 3 . In total, the EXAFS spectra and the fits for the entire set of Ln oxides (excluding promethium) are included. The knowledge developed here can assist in the structural determination of a wide variety of Ln compounds and can further studies on their structure–property relationships.

36 MATERIALS SCIENCE↗

Estimating uncertainty: A Bayesian approach to modelling photosynthesis in C3 leaves

The Farquhar-von Caemmerer-Berry (FvCB) model is extensively used to model pho-tosynthesis from gas exchange measurements. Since its publication, many methods have been developed to measure, or more accurately estimate, parameters of this model. Here, we have created a tool that uses Bayesian statistics to fit photosyn-thetic parameters using concurrent gas exchange and chlorophyll fluorescence mea-surements whilst evaluating the reliability of the parameter estimation. We have tested this tool on synthetic data and experimental data from rice leaves. Our results indicate that reliable parameter estimation can be achieved whilst only keeping one parameter, Km, that is, Michaelis constant for CO2 by Rubisco, prefixed. Additionally, we show that including detailed low CO2 measurements at low light levels increases reliability and suggests this as a new standard measurement protocol. By providing an estimated distribution of parameter values, the tool can be used to evaluate the quality of data from gas exchange and chlorophyll fluorescence measurement proto-cols. Compared to earlier model fitting methods, the use of a Bayesian statistics-based tool minimizes human interaction during fitting, reducing the subjectivity which is essential to most existing tools. A user friendly, interactive Bayesian tool script is provided.

Bayesian statistics, leaf photosynthesis, mesophyl↗

TRIZ-Based Design Improvement for Facilitating Transmission Control Unit Remanufacturing

Design for Remanufacturing (DfRem) centers on enhancing product design and operations to facilitate remanufacturing while increasing both economic and environmental sustainability. DfRem requires a comprehensive understanding of product characteristics, production conditions, and operational constraints. Thus, DfRem requires a systematic methodology to identify and implement alternative designs. Among many critical steps in remanufacturing, disassembly is essential because it directly supports remanufacturing by separation of product components for cleaning, inspection, and other subsequent process steps. This study proposes a TRIZ-based framework to improve product design for disassembly in support of remanufacturing. The framework is applied to a transmission control unit (TCU) case study. The current TCU design prevents remanufacturing due to the sealant-based component joinery, which complicates disassembly and risks damaging the printed circuit board (PCB). After defining technical contradictions for the TCU product design and reviewing the suggested TRIZ principles to solve the conflicts, a cantilever snap-fit design alternative is developed. The economic feasibility of the snap-fit design is assessed by comparing the current and snap-fit design costs for three life cycles. The cost analysis demonstrates that the cost of the snap-fit design remains the same as the current design. Additionally, the snap-fit design offers substantial cost savings for three life cycles compared to the current design. We demonstrate how snap-fit design supports both environmental and economic sustainability.

42 ENGINEERING↗

Data for KETCHUP: Parameterizing of Large-Scale Kinetic Models Using Multiple Datasets with Different Reference States

Repository for Kinetic Estimation Tool Capturing Heterogeneous Datasets Using Pyomo (KETCHUP), a flexible parameter estimation tool that leverages a primal-dual interior-point algorithm to solve a nonlinear programming (NLP) problem that identifies a set of parameters capable of recapitulating the steady-state fluxes and concentrations in wild-type and perturbed metabolic networks. KETCHUP can use K-FIT [2] input files. Example K-FIT input files are located in the K-FIT repository at https://github.com/maranasgroup/K-FIT.

Metabolomics↗

Analysis of picosecond coherent anti-Stokes Raman spectra for gas-phase diagnostics

We present a hybrid frequency- and time-domain solution, applicable to the case of picosecond coherent anti-Stokes Raman scattering (CARS), for gas-phase diagnostics. A solution has been derived based on both physical arguments and four-wave mixing equations for picosecond CARS, with pulse durations that are comparable to the dephasing time scale for gas-phase Raman coherence—a regime where commonly employed solutions for impulsive (femtosecond) or cw (nanosecond) pump/Stokes forcing are not strictly valid. We present the ps-CARS spectrum in the form of incoherent sums of CARS intensity spectra, calculated from the fundamental solution for impulsive pump/Stokes Raman preparation. The solution was examined for temperatures from 1000–3000 K, for four plausible experimental configurations, with laser pulse durations of 50–150 ps, and probe pulse delays from −20 to 240 ps. Approximations based on cw and impulsive pump/Stokes preparation to fit picosecond CARS spectra at atmospheric pressure were examined and the relative thermometric accuracy and computational cost of these approximations were quantified for the case of a zero nonresonant CARS contribution, and a nonresonant susceptibility equal to 10% of the Raman-resonant value at the N 2 bandhead. The nanosecond CARS approximation can result in large fitting errors when the probe pulse time delay is less than the probe pulse duration. Errors as large as 10–20% are observed in the fit temperatures for a zero picosecond probe pulse delay, when the nonresonant background is neglected, largely due to an inability of the time-independent cw model to capture transient frequency spread dephasing effects at the Q -branch bandhead. The inclusion of a nonresonant background results in 40–60% thermometry errors with a nanosecond model at a zero-probe delay. Time-dependent impulsive calculations used for femtosecond CARS better approximate the structure of the N 2 bandhead, reducing temperature fitting errors to 5–10% at a short probe pulse delay. The impulsive approximation results in errors up to 10% at intermediate probe pulse delays, where the coherence of the pump and probe pulses leads to multiple terms in the picosecond CARS solution. Both approximations improve as the probe pulse delay exceeds the probe duration. The nanosecond approximation results in a 2–3% error, while the impulsive model results in differences of less than 1% in some cases. Fits to experimental data obtained using short, ∼60ps pulses at a zero probe time delay and longer 100 ps pulses at a substantial 200 ps delay are presented with accuracies of 1–3% in the fit temperature.

Kearney, Sean P.↗

PBX9502 PAGOSA/SURF Calibration for 100 µm Grid Size

There are two methods to calibrate the SURF reactive burn explosives model.The first is to calibrate using Pop plot data. Because there are only two main parameters in SURF, called A and B, there is just enough flexibility in SURF to determine the Pop plot line, i.e., equivalent information to a slope and a single point on the line. However, due to the linearity of the Pop plot being based on logarithmic scales, a slight deviation off the Pop plot (though the fit to the Pop plot might be very good overall) can make considerable difference when comparing with the velocity profile data. The second calibration method thusly focuses strongly on matching entire velocity profile data. This is to say that individual velocity profiles extracted from embedded gages (in the case of PBX 9502, we refer to R. Gustavsen’s experimental data) may fit poorly, while the overall fit to the Pop plot may seem reasonable. In addition, there are other factors that could produce slightly different calibrations for SURF, while leading to an overall acceptable fit to data. One of those may be the physics code used and here we use Pagosa. Another is mesh size dependence. A third one stems from fitting a reactive burn model being done in conjunction with an equation of state (EOS) of the reactants and one for the products – different EOS can result in slightly different parameterizations. Finally, there are several secondary material parameters in SURF, and the overall fit will depend on the entire parameterization. The goal of this report is to clarify some of our previously made statements about SURF calibration, as well as to redo a systematic calibration using 100 µm grid resolution.

42 ENGINEERING↗

Non-conformal interface-cohesive modeling with the shifted boundary method

The accurate simulation of boundary- and interface-dominated problems on complex geometries remains challenging when boundary- or interface-fitted meshes are difficult to generate, particularly for curved boundaries, polycrystalline microstructures, and dense interface networks. The Shifted Boundary Method (SBM) alleviates this meshing burden by shifting the enforcement of boundary conditions from the true boundary to a nearby surrogate boundary and recovering the effect of the true boundary through geometric correction terms, thereby enabling standard finite element spaces on non-boundary-fitted meshes. In this report, we develop a general shiftedboundary and shifted-interface framework within the open-source MOOSE framework. We first present a general SBM implementation for complex geometries on non-boundary-fitted meshes. We then adopt the Shifted Interface Method (SIM) for internal interfaces and develop a unified shifted-interface treatment in which the interface law is enforced on a surrogate interface and the effect of the true interface is recovered through shifted jumps, fluxes, and tractions. This perspective brings scalar thermal-contact and vector-valued cohesive-zone mechanics into a single framework, the latter realized as the Shifted Cohesive Zone Method (SCZM) and coupled with history-dependent constitutive models from NEML2. We further extend the MOOSE mesh infrastructure to support cohesive-zone calculations on distributed meshes. The framework is verified and demonstrated through three progressive studies: Poisson’s equation on a smoothed starshaped domain, a manufactured thermal-contact problem on a non-interface-fitted mesh, and a two-dimensional polycrystalline representative volume element combining crystal plasticity with cohesive grain-boundary interfaces. Across these studies, the shifted formulations reproduce boundary- and interface-fitted reference solutions with high fidelity, indicating that the proposed framework provides an accurate and efficient route to boundary- and interface-dominated simulations on arbitrary geometries without requiring fitted meshes.

Yang, Cheng-Hau↗

The Gemini Planet Imager Exoplanet Survey: Dynamical Mass of the Exoplanet β Pictoris b from Combined Direct Imaging and Astrometry

Herein, we present new observations of the planet β Pictoris b from 2018 with the Gemini Planet Imager (GPI), the first GPI observations following conjunction. Based on these new measurements, we perform a joint orbit fit to the available relative astrometry from ground-based imaging, the Hipparcos Intermediate Astrometric Data (IAD), and the Gaia DR2 position, and demonstrate how to incorporate the IAD into direct imaging orbit fits. We find a mass consistent with predictions of hot-start evolutionary models and previous works following similar methods, though with larger uncertainties: $12.8^{+5.3}_{–3.2} M_{Jup}$. Our eccentricity determination of ${0.12}_{-0.03}^{+0.04}$ disfavors circular orbits. We consider orbit fits to several different imaging data sets, and find generally similar posteriors on the mass for each combination of imaging data. Our analysis underscores the importance of performing joint fits to the absolute and relative astrometry simultaneously, given the strong covariance between orbital elements. Time of conjunction is well-constrained within 2.8 days of 2017 September 13, with the star behind the planet's Hill sphere between 2017 April 11 and 2018 February 16 (±18 days). Following the recent radial velocity detection of a second planet in the system, β Pic c, we perform additional two-planet fits combining relative astrometry, absolute astrometry, and stellar radial velocities. These joint fits find a significantly smaller mass (8.0 ± 2.6 M Jup) for the imaged planet β Pic b, in a somewhat more circular orbit. We expect future ground-based observations to further constrain the visual orbit and mass of the planet in advance of the release of Gaia DR4.

79 ASTRONOMY AND ASTROPHYSICS↗