-
Kernelization Complexity of Solution Discovery Problems
Authors:
Mario Grobler,
Stephanie Maaz,
Amer E. Mouawad,
Naomi Nishimura,
Vijayaragunathan Ramamoorthi,
Sebastian Siebertz
Abstract:
In the solution discovery variant of a vertex (edge) subset problem $Π$ on graphs, we are given an initial configuration of tokens on the vertices (edges) of an input graph $G$ together with a budget $b$. The question is whether we can transform this configuration into a feasible solution of $Π$ on $G$ with at most $b$ modification steps. We consider the token sliding variant of the solution disco…
▽ More
In the solution discovery variant of a vertex (edge) subset problem $Π$ on graphs, we are given an initial configuration of tokens on the vertices (edges) of an input graph $G$ together with a budget $b$. The question is whether we can transform this configuration into a feasible solution of $Π$ on $G$ with at most $b$ modification steps. We consider the token sliding variant of the solution discovery framework, where each modification step consists of sliding a token to an adjacent vertex (edge). The framework of solution discovery was recently introduced by Fellows et al. [Fellows et al., ECAI 2023] and for many solution discovery problems the classical as well as the parameterized complexity has been established. In this work, we study the kernelization complexity of the solution discovery variants of Vertex Cover, Independent Set, Dominating Set, Shortest Path, Matching, and Vertex Cut with respect to the parameters number of tokens $k$, discovery budget $b$, as well as structural parameters such as pathwidth.
△ Less
Submitted 25 September, 2024;
originally announced September 2024.
-
Robust personalized pricing under uncertainty of purchase probabilities
Authors:
Shunnosuke Ikeda,
Naoki Nishimura,
Noriyoshi Sukegawa,
Yuichi Takano
Abstract:
This paper is concerned with personalized pricing models aimed at maximizing the expected revenues or profits for a single item. While it is essential for personalized pricing to predict the purchase probabilities for each consumer, these predicted values are inherently subject to unavoidable errors that can negatively impact the realized revenues and profits. To address this issue, we focus on ro…
▽ More
This paper is concerned with personalized pricing models aimed at maximizing the expected revenues or profits for a single item. While it is essential for personalized pricing to predict the purchase probabilities for each consumer, these predicted values are inherently subject to unavoidable errors that can negatively impact the realized revenues and profits. To address this issue, we focus on robust optimization techniques that yield reliable solutions to optimization problems under uncertainty. Specifically, we propose a robust optimization model for personalized pricing that accounts for the uncertainty of predicted purchase probabilities. This model can be formulated as a mixed-integer linear optimization problem, which can be solved exactly using mathematical optimization solvers. We also develop a Lagrangian decomposition algorithm combined with line search to efficiently find high-quality solutions for large-scale optimization problems. Experimental results demonstrate the effectiveness of our robust optimization model and highlight the utility of our Lagrangian decomposition algorithm in terms of both computational efficiency and solution quality.
△ Less
Submitted 21 July, 2024;
originally announced July 2024.
-
Strategic Coupon Allocation for Increasing Providers' Sales Experiences in Two-sided Marketplaces
Authors:
Koya Ohashi,
Sho Sekine,
Deddy Jobson,
Jie Yang,
Naoki Nishimura,
Noriyoshi Sukegawa,
Yuichi Takano
Abstract:
In a two-sided marketplace, network effects are crucial for competitiveness, and platforms need to retain users through advanced customer relationship management as much as possible. Maintaining numerous providers' stable and active presence on the platform is highly important to enhance the marketplace's scale and diversity. The strongest motivation for providers to continue using the platform is…
▽ More
In a two-sided marketplace, network effects are crucial for competitiveness, and platforms need to retain users through advanced customer relationship management as much as possible. Maintaining numerous providers' stable and active presence on the platform is highly important to enhance the marketplace's scale and diversity. The strongest motivation for providers to continue using the platform is to realize actual profits through sales. Then, we propose a personalized promotion to increase the number of successful providers with sales experiences on the platform. The main contributions of our research are twofold. First, we introduce a new perspective in provider management with the distribution of successful sales experiences. Second, we propose a personalized promotion optimization method to maximize the number of providers' sales experiences. By utilizing this approach, we ensure equal opportunities for providers to experience sales without being monopolized by a few providers. Through experiments using actual data on coupon distribution, we confirm that our method enables the implementation of coupon allocation strategies that significantly increase the total number of providers having sales experiences.
△ Less
Submitted 20 July, 2024;
originally announced July 2024.
-
Balancing Immediate Revenue and Future Off-Policy Evaluation in Coupon Allocation
Authors:
Naoki Nishimura,
Ken Kobayashi,
Kazuhide Nakata
Abstract:
Coupon allocation drives customer purchases and boosts revenue. However, it presents a fundamental trade-off between exploiting the current optimal policy to maximize immediate revenue and exploring alternative policies to collect data for future policy improvement via off-policy evaluation (OPE). To balance this trade-off, we propose a novel approach that combines a model-based revenue maximizati…
▽ More
Coupon allocation drives customer purchases and boosts revenue. However, it presents a fundamental trade-off between exploiting the current optimal policy to maximize immediate revenue and exploring alternative policies to collect data for future policy improvement via off-policy evaluation (OPE). To balance this trade-off, we propose a novel approach that combines a model-based revenue maximization policy and a randomized exploration policy for data collection. Our framework enables flexible adjustment of the mixture ratio between these two policies to optimize the balance between short-term revenue and future policy improvement. We formulate the problem of determining the optimal mixture ratio as multi-objective optimization, enabling quantitative evaluation of this trade-off. We empirically verified the effectiveness of the proposed mixed policy using synthetic data. Our main contributions are: (1) Demonstrating a mixed policy combining deterministic and probabilistic policies, flexibly adjusting the data collection vs. revenue trade-off. (2) Formulating the optimal mixture ratio problem as multi-objective optimization, enabling quantitative evaluation of this trade-off.
△ Less
Submitted 8 September, 2024; v1 submitted 5 July, 2024;
originally announced July 2024.
-
Closed-loop measurement in atom interferometer gyroscope with velocity-dependent phase dispersion compensation
Authors:
Tomoya Sato,
Naoki Nishimura,
Naoki Kaku,
Sotatsu Otabe,
Takuya Kawasaki,
Toshiyuki Hosoya,
Mikio Kozuma
Abstract:
Atom interferometer-based gyroscopes are expected to have a wide range of applications due to their high sensitivity. However, their dynamic range is limited by dephasing caused by velocity-dependent Sagnac phase shift in combination with the velocity distribution of the atoms, restricting measurements of large angular velocities. In this study, we present a method for restoring the contrast deter…
▽ More
Atom interferometer-based gyroscopes are expected to have a wide range of applications due to their high sensitivity. However, their dynamic range is limited by dephasing caused by velocity-dependent Sagnac phase shift in combination with the velocity distribution of the atoms, restricting measurements of large angular velocities. In this study, we present a method for restoring the contrast deterioration in rotation rate measurements with interferometer gyroscopes using atomic beams. Our findings confirm that by introducing the pseudo-rotation effect with appropriate two-photon detunings for Raman lights in the interferometer, it is possible to effectively cancel the rotation of all atoms in the velocity distribution of the beam. Consequently, the contrast is unaffected by the rotation. Furthermore, we applied this method to an interferometer gyroscope with counterpropagating atomic beams sharing the same Raman lights. We also found that the rotation rate of the system can be estimated through the detunings points where the phase difference between the two interferometers is zero. This approach ensures that the scale factor of the atom interferometer gyroscope is independent of the velocity fluctuation of the atomic beam. We demonstrate our technique using the interferometer gyroscope of thermal atomic beams of rubidium-87, achieving a measurement of rotation rate of $\mathrm{{1.0}^{\circ}/s}$ even with an acceleration of 0.68$\mathrm{m/s^2}$ on a three-axis rotation table. This simple and robust dispersion compensation method with Raman light detuning benefits dynamic rotation-rate measurements in field applications such as the inertial navigation of vehicles.
△ Less
Submitted 8 July, 2024;
originally announced July 2024.
-
Fast solution to the fair ranking problem using the Sinkhorn algorithm
Authors:
Yuki Uehara,
Shunnosuke Ikeda,
Naoki Nishimura,
Koya Ohashi,
Yilin Li,
Jie Yang,
Deddy Jobson,
Xingxia Zha,
Takeshi Matsumoto,
Noriyoshi Sukegawa,
Yuichi Takano
Abstract:
In two-sided marketplaces such as online flea markets, recommender systems for providing consumers with personalized item rankings play a key role in promoting transactions between providers and consumers. Meanwhile, two-sided marketplaces face the problem of balancing consumer satisfaction and fairness among items to stimulate activity of item providers. Saito and Joachims (2022) devised an impac…
▽ More
In two-sided marketplaces such as online flea markets, recommender systems for providing consumers with personalized item rankings play a key role in promoting transactions between providers and consumers. Meanwhile, two-sided marketplaces face the problem of balancing consumer satisfaction and fairness among items to stimulate activity of item providers. Saito and Joachims (2022) devised an impact-based fair ranking method for maximizing the Nash social welfare based on fair division; however, this method, which requires solving a large-scale constrained nonlinear optimization problem, is very difficult to apply to practical-scale recommender systems. We thus propose a fast solution to the impact-based fair ranking problem. We first transform the fair ranking problem into an unconstrained optimization problem and then design a gradient ascent method that repeatedly executes the Sinkhorn algorithm. Experimental results demonstrate that our algorithm provides fair rankings of high quality and is about 1000 times faster than application of commercial optimization software.
△ Less
Submitted 10 June, 2024;
originally announced June 2024.
-
Interpretable Price Bounds Estimation with Shape Constraints in Price Optimization
Authors:
Shunnosuke Ikeda,
Naoki Nishimura,
Shunji Umetani
Abstract:
This study addresses the interpretable estimation of price bounds in the context of price optimization. In recent years, price-optimization methods have become indispensable for maximizing revenue and profits. However, effective application of these methods to real-world pricing operations remains a significant challenge. It is crucial for operators responsible for setting prices to utilize reason…
▽ More
This study addresses the interpretable estimation of price bounds in the context of price optimization. In recent years, price-optimization methods have become indispensable for maximizing revenue and profits. However, effective application of these methods to real-world pricing operations remains a significant challenge. It is crucial for operators responsible for setting prices to utilize reasonable price bounds that are not only interpretable but also acceptable. Despite this necessity, most studies assume that price bounds are given constant values, and few have explored reasonable determinations of these bounds. Therefore, we propose a comprehensive framework for determining price bounds that includes both the estimation and adjustment of these bounds. Specifically, we first estimate price bounds using three distinct approaches based on historical pricing data. Then, we adjust the estimated price bounds by solving an optimization problem that incorporates shape constraints. This method allows the implementation of price optimization under practical and reasonable price bounds suitable for real-world applications. We report the effectiveness of our proposed method through numerical experiments using historical pricing data from actual services.
△ Less
Submitted 29 September, 2024; v1 submitted 23 May, 2024;
originally announced May 2024.
-
Robust portfolio optimization model for electronic coupon allocation
Authors:
Yuki Uehara,
Naoki Nishimura,
Yilin Li,
Jie Yang,
Deddy Jobson,
Koya Ohashi,
Takeshi Matsumoto,
Noriyoshi Sukegawa,
Yuichi Takano
Abstract:
Currently, many e-commerce websites issue online/electronic coupons as an effective tool for promoting sales of various products and services. We focus on the problem of optimally allocating coupons to customers subject to a budget constraint on an e-commerce website. We apply a robust portfolio optimization model based on customer segmentation to the coupon allocation problem. We also validate th…
▽ More
Currently, many e-commerce websites issue online/electronic coupons as an effective tool for promoting sales of various products and services. We focus on the problem of optimally allocating coupons to customers subject to a budget constraint on an e-commerce website. We apply a robust portfolio optimization model based on customer segmentation to the coupon allocation problem. We also validate the efficacy of our method through numerical experiments using actual data from randomly distributed coupons. Main contributions of our research are twofold. First, we handle six types of coupons, thereby making it extremely difficult to accurately estimate the difference in the effects of various coupons. Second, we demonstrate from detailed numerical results that the robust optimization model achieved larger uplifts of sales than did the commonly-used multiple-choice knapsack model and the conventional mean-variance optimization model. Our results open up great potential for robust portfolio optimization as an effective tool for practical coupon allocation.
△ Less
Submitted 21 May, 2024;
originally announced May 2024.
-
Reconfiguration of Multisets with Applications to Bin Packing
Authors:
Jeffrey Kam,
Shahin Kamali,
Avery Miller,
Naomi Nishimura
Abstract:
We use the reconfiguration framework to analyze problems that involve the rearrangement of items among groups. In various applications, a group of items could correspond to the files or jobs assigned to a particular machine, and the goal of rearrangement could be improving efficiency or increasing locality.
To cover problems arising in a wide range of application areas, we define the general Rep…
▽ More
We use the reconfiguration framework to analyze problems that involve the rearrangement of items among groups. In various applications, a group of items could correspond to the files or jobs assigned to a particular machine, and the goal of rearrangement could be improving efficiency or increasing locality.
To cover problems arising in a wide range of application areas, we define the general Repacking problem as the rearrangement of multisets of multisets. We present hardness results for the general case and algorithms for various classes of instances that arise in real-life scenarios. By limiting the total size of items in each multiset, our results can be viewed as an offline approach to Bin Packing, in which each bin is represented as a multiset.
In addition to providing the first results on reconfiguration of multisets, our contributions open up several research avenues: the interplay between reconfiguration and online algorithms and parallel algorithms; the use of the tools of linear programming in reconfiguration; and, in the longer term, a focus on resources in reconfiguration.
△ Less
Submitted 9 May, 2024;
originally announced May 2024.
-
Covariate Ordered Systematic Sampling as an Improvement to Randomized Controlled Trials
Authors:
Deddy Jobson,
Li Yilin,
Naoki Nishimura,
Yang Jie,
Koya Ohashi,
Takeshi Matsumoto
Abstract:
The Randomized Controlled Trial (RCT) or A/B testing is considered the gold standard method for estimating causal effects. Fisher famously advocated randomly allocating experiment units into treatment and control groups to preclude systematic biases. We propose a variant of systematic sampling called Covariate Ordered Systematic Sampling (COSS). In COSS, we order experimental units using a pre-exp…
▽ More
The Randomized Controlled Trial (RCT) or A/B testing is considered the gold standard method for estimating causal effects. Fisher famously advocated randomly allocating experiment units into treatment and control groups to preclude systematic biases. We propose a variant of systematic sampling called Covariate Ordered Systematic Sampling (COSS). In COSS, we order experimental units using a pre-experiment covariate and allocate them alternately into treatment and control groups. Using theoretical proofs, experiments on simulated data, and hundreds of A/B tests conducted within 3 real-world marketing campaigns, we show how our method achieves better sensitivity gains than commonly used variance reduction techniques like CUPED while retaining the simplicity of RCTs.
△ Less
Submitted 16 April, 2024;
originally announced April 2024.
-
Constraints on the Neutron-Star Structure from the Clocked X-Ray Burster 1RXS J180408.9$-$342058
Authors:
Akira Dohi,
Wataru Iwakiri,
Nobuya Nishimura,
Tsuneo Noda,
Shigehiro Nagataki,
Masa-aki Hashimoto
Abstract:
Type-I X-ray bursts are rapid-brightening transient phenomena on the surfaces of accreting neutron stars (NSs). Some X-ray bursts, called {\it clocked bursters}, exhibit regular behavior with similar light curve profiles in their burst sequences. The periodic nature of clocked bursters has the advantage of constraining X-ray binary parameters and physics inside the NS. In the present study, we com…
▽ More
Type-I X-ray bursts are rapid-brightening transient phenomena on the surfaces of accreting neutron stars (NSs). Some X-ray bursts, called {\it clocked bursters}, exhibit regular behavior with similar light curve profiles in their burst sequences. The periodic nature of clocked bursters has the advantage of constraining X-ray binary parameters and physics inside the NS. In the present study, we compute numerical models, based on different equations of state and NS masses, which are compared with the observation of a recently identified clocked burster, 1RXS J180408.9$-$342058. We find that the relation between accretion rate and recurrence time is highly sensitive to the NS mass and radius. We determine, in particular, that 1RXS J180408.9$-$342058 appears to possess a mass less than $1.7M_{\odot}$ and favors a stiffer nuclear equation of state (with an NS radius $\gtrsim12.7{\rm km}$). Consequently, the observations of this new clocked burster may provide additional constraints for probing the structure of NSs.
△ Less
Submitted 5 November, 2023;
originally announced November 2023.
-
An IPW-based Unbiased Ranking Metric in Two-sided Markets
Authors:
Keisho Oh,
Naoki Nishimura,
Minje Sung,
Ken Kobayashi,
Kazuhide Nakata
Abstract:
In modern recommendation systems, unbiased learning-to-rank (LTR) is crucial for prioritizing items from biased implicit user feedback, such as click data. Several techniques, such as Inverse Propensity Weighting (IPW), have been proposed for single-sided markets. However, less attention has been paid to two-sided markets, such as job platforms or dating services, where successful conversions requ…
▽ More
In modern recommendation systems, unbiased learning-to-rank (LTR) is crucial for prioritizing items from biased implicit user feedback, such as click data. Several techniques, such as Inverse Propensity Weighting (IPW), have been proposed for single-sided markets. However, less attention has been paid to two-sided markets, such as job platforms or dating services, where successful conversions require matching preferences from both users. This paper addresses the complex interaction of biases between users in two-sided markets and proposes a tailored LTR approach. We first present a formulation of feedback mechanisms in two-sided matching platforms and point out that their implicit feedback may include position bias from both user groups. On the basis of this observation, we extend the IPW estimator and propose a new estimator, named two-sided IPW, to address the position bases in two-sided markets. We prove that the proposed estimator satisfies the unbiasedness for the ground-truth ranking metric. We conducted numerical experiments on real-world two-sided platforms and demonstrated the effectiveness of our proposed method in terms of both precision and robustness. Our experiments showed that our method outperformed baselines especially when handling rare items, which are less frequently observed in the training data.
△ Less
Submitted 13 July, 2023;
originally announced July 2023.
-
Postfission properties of uranium isotopes: A hybrid method with Langevin dynamics and the Hauser-Feshbach statistical model
Authors:
Shyoya Tanaka,
Nobuya Nishimura,
Futoshi Minato,
Yoshihiro Aritomo
Abstract:
Background: Precise understanding of nuclear fission is crucial for experimental and theoretical nuclear physics, astrophysics, and industrial applications; however, the complete physical mechanics is unresolved due to the complexities.
Purpose: In this study, we present a new method to describe the dynamical-fission process and following prompt-neutron emission, where we combine the dynamical f…
▽ More
Background: Precise understanding of nuclear fission is crucial for experimental and theoretical nuclear physics, astrophysics, and industrial applications; however, the complete physical mechanics is unresolved due to the complexities.
Purpose: In this study, we present a new method to describe the dynamical-fission process and following prompt-neutron emission, where we combine the dynamical fission calculation based on the Langevin method and the Hauser-Feshbach statistical model.
Methods: Two methods are connected smoothly within the universal charge distribution and the energy conservation, allowing us to calculate a sequence of fission dynamics and post-fission phase, including prompt neutron emission.
Results: Using a certain set of model parameters, we successfully reproduce the experimental primary-fission yields, total kinetic energy, independent-fission yields, and prompt neutron emissions for the neutron induced fission of ${}^{236}$U, a compound nucleus of ${\rm n} + {}^{235}{\rm U}$. We elucidate the physical mechanism of the characteristic features observed in previous experiments, such as shell properties. Additionally, we apply our calculation to two very neutron-rich uranium isotopes, i.e., ${}^{250}$U and ${}^{255}$U, which are not experimentally confirmed but are important for r-process nucleosynthesis. Theoretical results indicate that ${}^{250}$U exhibits an asymmetric multiple-peak fission yield distribution, while the neutron-rich ${}^{255}$U has a single peak due to symmetric fission. Our method predicts post-neutron emission fragments, where ${}^{250}$U shows a stronger neutron emissivity than ${}^{255}$U.
Conclusions: Our framework is highly reproducible in the experiments and shows that the number of emitted neutrons after fission differs significantly in neutron-rich uranium fission depending on distributions of fission variables.
△ Less
Submitted 23 November, 2023; v1 submitted 18 July, 2023;
originally announced July 2023.
-
The Impacts of Neutron-Star Structure and Base Heating on Type I X-Ray Bursts and Code Comparison
Authors:
Guoqing Zhen,
Guoliang Lv,
Helei Liu,
Akira Dohi,
Nobuya Nishimura,
Chunhua Zhu,
Liyu Song,
Weiyang Wang,
Renxin Xu
Abstract:
Type I X-ray bursts are rapidly brightening phenomena triggered by thermonuclear burning on accreting layer of a neutron star (NS). The light curves represent the physical properties of NSs and the nuclear reactions on the proton-rich nuclei. The numerical treatments of the accreting NS and physics of the NS interior are not established, which shows uncertainty in modelling for observed X-ray ligh…
▽ More
Type I X-ray bursts are rapidly brightening phenomena triggered by thermonuclear burning on accreting layer of a neutron star (NS). The light curves represent the physical properties of NSs and the nuclear reactions on the proton-rich nuclei. The numerical treatments of the accreting NS and physics of the NS interior are not established, which shows uncertainty in modelling for observed X-ray light curves. In this study, we investigate theoretical X-ray-burst models, compared with burst light curves with GS~1826-24 observations. We focus on the impacts of the NS mass, the NS radius, and base-heating on the NS surface using the MESA code. We find a monotonic correlation between the NS mass and the parameters of the light curve. The higher the mass, the longer the recurrence time and the greater the peak luminosity. While the larger the radius, the longer the recurrence time, the peak luminosity remains nearly constant. In the case of increasing base heating, both the recurrence time and peak luminosity decrease. We also examine the above results using with a different numerical code, HERES, based on general relativity and consider the central NS. We find that the burst rate, burst energy and burst strength are almost same in two X-ray burst codes by adjusting the base-heat parameter in MESA (the relative errors $\lesssim5\%$), while the duration time and the rise time are significantly different between (the relative error is possibly $\sim50\%$). The peak luminosity and the e-folding time are ragged between two codes for different accretion rates.
△ Less
Submitted 14 April, 2023;
originally announced April 2023.
-
Cross Sections of the $^{83}$Rb(p,$γ)^{84}$Sr and $^{84}$Kr(p,$γ)^{85}$Rb Reactions at Energies Characteristic of the Astrophysical $γ$ Process
Authors:
M. Williams,
B. Davids,
G. Lotay,
N. Nishimura,
T. Rauscher,
S. A. Gillespie,
M. Alcorta,
A. M. Amthor,
G. C. Ball,
S. S. Bhattacharjee,
V. Bildstein,
W. N. Catford,
D. T. Doherty,
N. E. Esker,
A. B. Garnsworthy,
G. Hackman,
K. Hudson,
A. Lennarz,
C. Natzke,
B. Olaizola,
A. Psaltis,
C. E. Svensson,
J. Williams,
D. Walter,
D. Yates
Abstract:
We have measured the cross section of the $^{83}$Rb(p,$γ)^{84}$Sr radiative capture reaction in inverse kinematics using a radioactive beam of $^{83}$Rb at incident energies of 2.4 and $2.7 A$ MeV. Prior to the radioactive beam measurement, the $^{84}$Kr(p,$γ)^{85}$Rb radiative capture reaction was measured in inverse kinematics using a stable beam of $^{84}$Kr at an incident energy of $2.7 A$ MeV…
▽ More
We have measured the cross section of the $^{83}$Rb(p,$γ)^{84}$Sr radiative capture reaction in inverse kinematics using a radioactive beam of $^{83}$Rb at incident energies of 2.4 and $2.7 A$ MeV. Prior to the radioactive beam measurement, the $^{84}$Kr(p,$γ)^{85}$Rb radiative capture reaction was measured in inverse kinematics using a stable beam of $^{84}$Kr at an incident energy of $2.7 A$ MeV. The effective relative kinetic energies of these measurements lie within the relevant energy window for the $γ$ process in supernovae. The central values of the measured partial cross sections of both reactions were found to be $0.17-0.42$ times the predictions of statistical model calculations. Assuming the predicted cross section at other energies is reduced by the same factor leads to a slightly higher calculated abundance of the $p$ nucleus $^{84}$Sr, caused by the reduced rate of the $^{84}$Sr($γ$,p)$^{83}$Rb reaction derived from the present measurement.
△ Less
Submitted 17 March, 2023;
originally announced March 2023.
-
Efficient algorithms to solve atom reconfiguration problems. II. The assignment-rerouting-ordering (aro) algorithm
Authors:
Remy El Sabeh,
Jessica Bohm,
Zhiqian Ding,
Stephanie Maaz,
Naomi Nishimura,
Izzat El Hajj,
Amer E. Mouawad,
Alexandre Cooper
Abstract:
Programmable arrays of optical traps enable the assembly of configurations of single atoms to perform controlled experiments on quantum many-body systems. Finding the sequence of control operations to transform an arbitrary configuration of atoms into a predetermined one requires solving an atom reconfiguration problem quickly and efficiently. A typical approach to solve atom reconfiguration probl…
▽ More
Programmable arrays of optical traps enable the assembly of configurations of single atoms to perform controlled experiments on quantum many-body systems. Finding the sequence of control operations to transform an arbitrary configuration of atoms into a predetermined one requires solving an atom reconfiguration problem quickly and efficiently. A typical approach to solve atom reconfiguration problems is to use an assignment algorithm to determine which atoms to move to which traps. This approach results in control protocols that exactly minimize the number of displacement operations; however, this approach does not optimize for the number of displaced atoms nor the number of times each atom is displaced, resulting in unnecessary control operations that increase the execution time and failure rate of the control protocol. In this work, we propose the assignment-rerouting-ordering (aro) algorithm to improve the performance of assignment-based algorithms in solving atom reconfiguration problems. The aro algorithm uses an assignment subroutine to minimize the total distance traveled by all atoms, a rerouting subroutine to reduce the number of displaced atoms, and an ordering subroutine to guarantee that each atom is displaced at most once. The ordering subroutine relies on the existence of a partial ordering of moves that can be obtained using a polynomial-time algorithm that we introduce within the formal framework of graph theory. We numerically quantify the performance of the aro algorithm in the presence and in the absence of loss, and show that it outperforms the exact, approximation, and heuristic algorithms that we use as benchmarks. Our results are useful for assembling large configurations of atoms with high success probability and fast preparation time, as well as for designing and benchmarking novel atom reconfiguration algorithms.
△ Less
Submitted 11 December, 2022;
originally announced December 2022.
-
Efficient algorithms to solve atom reconfiguration problems. I. The redistribution-reconfiguration (red-rec) algorithm
Authors:
Barry Cimring,
Remy El Sabeh,
Marc Bacvanski,
Stephanie Maaz,
Izzat El Hajj,
Naomi Nishimura,
Amer E. Mouawad,
Alexandre Cooper
Abstract:
We propose the redistribution-reconfiguration~(red-rec) algorithm to efficiently compute control protocols to assemble compact-centered configurations of atoms in two-dimensional arrays of optical traps with lattice geometries. The red-rec algorithm redistributes atoms among pairs of donor-receiver columns and reconfigures each column using an exact displacement-minimizing algorithm, harnessing pa…
▽ More
We propose the redistribution-reconfiguration~(red-rec) algorithm to efficiently compute control protocols to assemble compact-centered configurations of atoms in two-dimensional arrays of optical traps with lattice geometries. The red-rec algorithm redistributes atoms among pairs of donor-receiver columns and reconfigures each column using an exact displacement-minimizing algorithm, harnessing parallel control operations that simultaneously actuate multiple traps to reduce the execution time. We numerically quantify the performance of the red-rec algorithm, both in the absence and in the presence of loss, using realistic physical parameters and operational constraints. We show that the number of traps required to prepare a compact-centered configuration of atoms on a grid with a mean success probability of one half scales as the 3/2 power of the number of desired atoms, highlighting the challenges of assembling configurations of tens of thousands of atoms. We further demonstrate that faster preparation times can be achieved by rejecting configurations of atoms containing fewer atoms than a given threshold. The red-rec algorithm admits an efficient implementation that can readily be deployed on real-time control systems to assemble large configurations of atoms with high mean success probability and fast preparation times.
△ Less
Submitted 7 December, 2022;
originally announced December 2022.
-
Impacts of the direct URCA and Superfluidity inside a Neutron Star on Type-I X-Ray Bursts and X-Ray Superbursts
Authors:
Akira Dohi,
Nobuya Nishimura,
Hajime Sotani,
Tsuneo Noda,
Helei Liu,
Shigehiro Nagataki,
Masa-aki Hashimoto
Abstract:
We investigate the impacts of neutrino cooling mechanism inside the neutron star (NS) core on the light curves of type-I X-ray bursts and X-ray superbursts. From several observations of NS thermal evolution, physical processes of fast neutrino cooling, such as the direct Urca (DU) process, are indicated. They significantly decrease the surface temperature of NSs, though the cooling effect could be…
▽ More
We investigate the impacts of neutrino cooling mechanism inside the neutron star (NS) core on the light curves of type-I X-ray bursts and X-ray superbursts. From several observations of NS thermal evolution, physical processes of fast neutrino cooling, such as the direct Urca (DU) process, are indicated. They significantly decrease the surface temperature of NSs, though the cooling effect could be suppressed by nucleon superfluidity. In the present study, focusing on the DU process and nucleon superfluidity, we investigate the effects of NS cooling on the X-ray bursts using a general-relativistic stellar-evolution code. We find that the DU process leads find the longer recurrence time and the higher peak luminosity, which could be obstructed by the neutrons superfluidity. We also apply our burst models to the comparison with {\it Clocked burster} GS 1826$-$24, and to the recurrence time of superburst triggered by carbon ignition. These effects are significant within a certain range of binary parameters and uncertainty of the NS equation of state.
△ Less
Submitted 20 March, 2023; v1 submitted 31 August, 2022;
originally announced August 2022.
-
New self-force method via elliptic partial differential equations for Kerr inspiral models
Authors:
Thomas Osburn,
Nami Nishimura
Abstract:
We present a new method designed to avoid numerical challenges that have impeded calculation of the Lorenz gauge self-force acting on a compact object inspiraling into a Kerr black hole. This type of calculation is valuable in creating waveform templates for extreme mass-ratio inspirals, which are an important source of gravitational waves for the upcoming Laser Interferometer Space Antenna missio…
▽ More
We present a new method designed to avoid numerical challenges that have impeded calculation of the Lorenz gauge self-force acting on a compact object inspiraling into a Kerr black hole. This type of calculation is valuable in creating waveform templates for extreme mass-ratio inspirals, which are an important source of gravitational waves for the upcoming Laser Interferometer Space Antenna mission. Prior hyperbolic partial differential equation (PDE) formulations encountered numerical instabilities involving unchecked growth in time; our new method is based on elliptic PDEs, which do not exhibit instabilities of that kind. For proof of concept, we calculate the self-force acting on a scalar charge in a circular orbit around a Kerr black hole. We anticipate this method will subsequently facilitate calculation of first-order Lorenz gauge Kerr metric perturbations and self-force, which could serve as a foundation for second-order Kerr self-force investigations.
△ Less
Submitted 6 September, 2022; v1 submitted 14 June, 2022;
originally announced June 2022.
-
Horizons: Nuclear Astrophysics in the 2020s and Beyond
Authors:
H. Schatz,
A. D. Becerril Reyes,
A. Best,
E. F. Brown,
K. Chatziioannou,
K. A. Chipps,
C. M. Deibel,
R. Ezzeddine,
D. K. Galloway,
C. J. Hansen,
F. Herwig,
A. P. Ji,
M. Lugaro,
Z. Meisel,
D. Norman,
J. S. Read,
L. F. Roberts,
A. Spyrou,
I. Tews,
F. X. Timmes,
C. Travaglio,
N. Vassh,
C. Abia,
P. Adsley,
S. Agarwal
, et al. (140 additional authors not shown)
Abstract:
Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilit…
▽ More
Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities.
△ Less
Submitted 16 May, 2022;
originally announced May 2022.
-
A survey on the parameterized complexity of the independent set and (connected) dominating set reconfiguration problems
Authors:
Nicolas Bousquet,
Amer E. Mouawad,
Naomi Nishimura,
Sebastian Siebertz
Abstract:
A graph vertex-subset problem defines which subsets of the vertices of an input graph are feasible solutions. We view a feasible solution as a set of tokens placed on the vertices of the graph. A reconfiguration variant of a vertex-subset problem asks, given two feasible solutions of size $k$, whether it is possible to transform one into the other by a sequence of token slides (along edges of the…
▽ More
A graph vertex-subset problem defines which subsets of the vertices of an input graph are feasible solutions. We view a feasible solution as a set of tokens placed on the vertices of the graph. A reconfiguration variant of a vertex-subset problem asks, given two feasible solutions of size $k$, whether it is possible to transform one into the other by a sequence of token slides (along edges of the graph) or token jumps (between arbitrary vertices of the graph) such that each intermediate set remains a feasible solution of size $k$. Many algorithmic questions present themselves in the form of reconfiguration problems: Given the description of an initial system state and the description of a target state, is it possible to transform the system from its initial state into the target one while preserving certain properties of the system in the process? Such questions have received a substantial amount of attention under the so-called combinatorial reconfiguration framework. We consider reconfiguration variants of three fundamental underlying graph vertex-subset problems, namely Independent Set, Dominating Set, and Connected Dominating Set. We survey both older and more recent work on the parameterized complexity of all three problems when parameterized by the number of tokens $k$. The emphasis will be on positive results and the most common techniques for the design of fixed-parameter tractable algorithms.
△ Less
Submitted 22 April, 2022;
originally announced April 2022.
-
New constraints on the neutron-star mass and radius relation from terrestrial nuclear experiments
Authors:
Hajime Sotani,
Nobuya Nishimura,
Tomoya Naito
Abstract:
The determination of the equation of state (EOS) for nuclear matter has been one of the biggest problems in nuclear astrophysics, because the EOS is essential for determining the properties of neutron stars. To constrain the density-dependence of the nuclear symmetry energy, several nuclear experiments, e.g., reported by the S$π$RIT and PREX-II collaborations, have recently been performed. However…
▽ More
The determination of the equation of state (EOS) for nuclear matter has been one of the biggest problems in nuclear astrophysics, because the EOS is essential for determining the properties of neutron stars. To constrain the density-dependence of the nuclear symmetry energy, several nuclear experiments, e.g., reported by the S$π$RIT and PREX-II collaborations, have recently been performed. However, since their uncertainties are still large, additional constraints such as astronomical observations must be crucial. In addition, it is interesting to see the effect of their reported value on neutron star properties. In this study, focusing on the relatively lower density region, we investigate the allowed area of the neutron-star mass and radius relation by assuming the constraints from S$π$RIT and PREX-II. Each region predicted by these experiments is still consistent with the allowed area constrained by the various astronomical observations. Our results show that terrestrial nuclear experiments must provide further constraints on the EOS for neutron stars, complementing astronomical observations.
△ Less
Submitted 16 March, 2022; v1 submitted 20 February, 2022;
originally announced March 2022.
-
Characteristic Basis Function Method Combined with Calderón Multiplicative Preconditioner for PMCHWT Formulation
Authors:
Tai Tanaka,
Kazuki Niino,
Naoshi Nishimura
Abstract:
We propose a novel characteristic basis function method for analyzing the scattering by dielectric objects based on the Poggio-Miller-Chang-Harrington-Wu-Tsai formulation. In the proposed method, the electric and magnetic currents are orthogonalized with the help of the singular value decomposition, and are used as dual basis functions in a way similar to the RWG and BC basis functions. We show th…
▽ More
We propose a novel characteristic basis function method for analyzing the scattering by dielectric objects based on the Poggio-Miller-Chang-Harrington-Wu-Tsai formulation. In the proposed method, the electric and magnetic currents are orthogonalized with the help of the singular value decomposition, and are used as dual basis functions in a way similar to the RWG and BC basis functions. We show that the use of the Calderón multiplicative precondtioner together with the proposed method can prevent from the poor convergence of the solution of the matrix equation in problems involving dielectrics. We considered three different shapes of dielectric scatterers for the purpose of validation. The numerical results agreed well with those obtained by the conventional method of moments and the proposed method was faster than the conventional method. These results indicate that the proposed method is effective for scattering analysis of the dielectrics.
△ Less
Submitted 22 November, 2021;
originally announced November 2021.
-
First Direct Measurement of an Astrophysical p-Process Reaction Cross Section Using a Radioactive Ion Beam
Authors:
G. Lotay,
S. A. Gillespie,
M. Williams,
T. Rauscher,
M. Alcorta,
M. Amthor,
C. A. Andreoiu,
D. Baal,
G. C. Ball,
S. S. Bhattacharjee,
H. Behnamian,
V. Bildstein,
C. Burbadge,
W. N. Catford,
D. T. Doherty,
N. E. Esker,
F. H. Garcia,
A. B. Garnsworthy,
G. Hackman,
S. Hallam,
K. A. Hudson,
S. Jazrawi,
E. Kasanda,
A. R. L. Kennington,
Y. H. Kim
, et al. (13 additional authors not shown)
Abstract:
We have performed the first direct measurement of the 83Rb(p,g) radiative capture reaction cross section in inverse kinematics using a radioactive beam of 83Rb at incident energies of 2.4 and 2.7 A MeV. The measured cross section at an effective relative kinetic energy of Ecm = 2.393 MeV, which lies within the relevant energy window for core collapse supernovae, is smaller than the prediction of s…
▽ More
We have performed the first direct measurement of the 83Rb(p,g) radiative capture reaction cross section in inverse kinematics using a radioactive beam of 83Rb at incident energies of 2.4 and 2.7 A MeV. The measured cross section at an effective relative kinetic energy of Ecm = 2.393 MeV, which lies within the relevant energy window for core collapse supernovae, is smaller than the prediction of statistical model calculations. This leads to the abundance of 84Sr produced in the astrophysical p process being higher than previously calculated. Moreover, the discrepancy of the present data with theoretical predictions indicates that further experimental investigation of p-process reactions involving unstable projectiles is clearly warranted.
△ Less
Submitted 14 September, 2021;
originally announced September 2021.
-
Parameterized complexity of reconfiguration of atoms
Authors:
Alexandre Cooper,
Stephanie Maaz,
Amer E. Mouawad,
Naomi Nishimura
Abstract:
Our work is motivated by the challenges presented in preparing arrays of atoms for use in quantum simulation. The recently-developed process of loading atoms into traps results in approximately half of the traps being filled. To consolidate the atoms so that they form a dense and regular arrangement, such as all locations in a grid, atoms are rearranged using moving optical tweezers. Time is of th…
▽ More
Our work is motivated by the challenges presented in preparing arrays of atoms for use in quantum simulation. The recently-developed process of loading atoms into traps results in approximately half of the traps being filled. To consolidate the atoms so that they form a dense and regular arrangement, such as all locations in a grid, atoms are rearranged using moving optical tweezers. Time is of the essence, as the longer that the process takes and the more that atoms are moved, the higher the chance that atoms will be lost in the process.
Viewed as a problem on graphs, we wish to solve the problem of reconfiguring one arrangement of tokens (representing atoms) to another using as few moves as possible. Because the problem is NP-complete on general graphs as well as on grids, we focus on the parameterized complexity for various parameters, considering both undirected and directed graphs, and tokens with and without labels. For unlabelled tokens, the problem is in FPT when parameterizing by the number of tokens, the number of moves, or the number of moves plus the number of vertices without tokens in either the source or target configuration, but intractable when parameterizing by the difference between the number of moves and the number of differences in the placement of tokens in the source and target configurations. When labels are added to tokens, however, most of the tractability results are replaced by hardness results.
△ Less
Submitted 26 July, 2021;
originally announced July 2021.
-
Effects of Nuclear Equation of State on Type-I X-ray Bursts: Interpretation of the X-ray Bursts from GS 1826-24
Authors:
A. Dohi,
N. Nishimura,
M. Hashimoto,
Y. Matsuo,
T. Noda,
S. Nagataki
Abstract:
Type I X-ray bursts are thermonuclear explosions on the neutron star (NS) surface caused by mass accretion from a companion star. Observations of X-ray bursts provide valuable information on X-ray binary systems, e.g., binary parameters, the chemical composition of accreted matter, and the nuclear equation of state (EOS) of NSs. There have been several theoretical studies to constrain the physics…
▽ More
Type I X-ray bursts are thermonuclear explosions on the neutron star (NS) surface caused by mass accretion from a companion star. Observations of X-ray bursts provide valuable information on X-ray binary systems, e.g., binary parameters, the chemical composition of accreted matter, and the nuclear equation of state (EOS) of NSs. There have been several theoretical studies to constrain the physics of X-ray bursters. However, they have mainly focused on the burning layers above the solid crust of the NS, which brings up issues of the treatment of NS gravitation and internal energy. In this study, focusing on the microphysics inside NSs, we calculate a series of X-ray bursts using a general-relativistic stellar-evolution code with several NS EOSs. We compare the X-ray-burst models with the burst parameters of a clocked burster associated with GS 1826-24. We find a monotonic correlation between the NS radius and the light-curve profile. A larger radius shows a higher recurrence time and a large peak luminosity. In contrast, the dependence of light curves on the NS mass becomes more complicated, where the neutrino cooling suppresses the efficiency of nuclear ignition. We also constrain the EOS and mass of GS 1826-24, i.e., stiffer EOSs, corresponding to larger NS radii, are not preferred due to too-high peak luminosity. The EOS and the cooling and heating of NSs are important to discuss the theoretical and observational properties of X-ray bursts.
△ Less
Submitted 5 January, 2022; v1 submitted 28 May, 2021;
originally announced May 2021.
-
Fission fragment distributions of neutron-rich nuclei based on Langevin calculations: toward r-process simulations
Authors:
Mizuki Okubayashi,
Shoya Tanaka,
Yoshihiro Aritomo,
Shoma Ishizaki,
Shota Amano,
Nobuya Nishimura
Abstract:
The nuclear fission of very neuron-rich nuclei related to the r-process is essential for the termination of nucleosynthesis flows on the nuclear chart and the final abundances. Nevertheless, most of the available fission data for the r-process calculations are based on theory predictions, including phenomenological treatments. In this study, we calculated a series of nuclear fission distribution f…
▽ More
The nuclear fission of very neuron-rich nuclei related to the r-process is essential for the termination of nucleosynthesis flows on the nuclear chart and the final abundances. Nevertheless, most of the available fission data for the r-process calculations are based on theory predictions, including phenomenological treatments. In this study, we calculated a series of nuclear fission distribution for neutron-rich nuclei away from the beta-stability line. As most of these nuclei are experimentally unknown, we are based on theoretical calculations based on the dynamical fission model with the Langevin method. We performed fission distribution calculations for neutron-rich actinoid nuclei, applicable to the r-process nucleosynthesis simulations. In the present paper, we compared the obtained mass and charge distributions with experimental data. We also show the results of the systematic behaviour of mass distribution for neutron-rich U and Fm isotopes.
△ Less
Submitted 19 May, 2021;
originally announced May 2021.
-
Statistical Analysis of the Relation between Coronal Mass Ejections and Solar Energetic Particles
Authors:
Kosuke Kihara,
Yuwei Huang,
Nobuhiko Nishimura,
Nariaki V. Nitta,
Seiji Yashiro,
Kiyoshi Ichimoto,
Ayumi Asai
Abstract:
To improve the forecasting capability of impactful solar energetic particle (SEP) events, the relation between coronal mass ejections (CMEs) and SEP events needs to be better understood. Here we present a statistical study of SEP occurrences and timescales with respect to the CME source locations and speeds, considering all 257 fast ($v_{CME}$ $\ge$ 900 km/s) and wide (angular width $\ge$ 60…
▽ More
To improve the forecasting capability of impactful solar energetic particle (SEP) events, the relation between coronal mass ejections (CMEs) and SEP events needs to be better understood. Here we present a statistical study of SEP occurrences and timescales with respect to the CME source locations and speeds, considering all 257 fast ($v_{CME}$ $\ge$ 900 km/s) and wide (angular width $\ge$ 60$^{\circ}$) CMEs that occurred between December 2006 and October 2017. We associate them with SEP events at energies above 10 MeV. Examination of the source region of each CME reveals that CMEs more often accompany a SEP event if they originate from the longitude of E20-W100 relative to the observer. However, a SEP event could still be absent if the CME is $<$ 2000 km/s. For the associated CME-SEP pairs, we compute three timescales for each of the SEP events, following Kahler (2005, 2013); namely the timescale of the onset (TO), the rise time (TR), and the duration (TD). They are correlated with the longitude of the CME source region relative to the footpoint of the Parker spiral ($ΔΦ$) and $v_{CME}$. The TO tends to be short for $|ΔΦ|$ $<$ 60$^{\circ}$ . This trend is weaker for TR and TD. The SEP timescales are only weakly correlated with $v_{CME}$. Positive correlations of both TR and TD with $v_{CME}$ are seen in poorly connected (large $|ΔΦ|$) events. Additionally, TO appears to be negatively correlated with $v_{CME}$ for events with small $|ΔΦ|$.
△ Less
Submitted 15 July, 2020;
originally announced July 2020.
-
Predicting Online Item-choice Behavior: A Shape-restricted Regression Perspective
Authors:
Naoki Nishimura,
Noriyoshi Sukegawa,
Yuichi Takano,
Jiro Iwanaga
Abstract:
This paper examines the relationship between user pageview (PV) histories and their item-choice behavior on an e-commerce website. We focus on PV sequences, which represent time series of the number of PVs for each user--item pair. We propose a shape-restricted optimization model that accurately estimates item-choice probabilities for all possible PV sequences. This model imposes monotonicity cons…
▽ More
This paper examines the relationship between user pageview (PV) histories and their item-choice behavior on an e-commerce website. We focus on PV sequences, which represent time series of the number of PVs for each user--item pair. We propose a shape-restricted optimization model that accurately estimates item-choice probabilities for all possible PV sequences. This model imposes monotonicity constraints on item-choice probabilities by exploiting partial orders for PV sequences, according to the recency and frequency of a user's previous PVs. To improve the computational efficiency of our optimization model, we devise efficient algorithms for eliminating all redundant constraints according to the transitivity of the partial orders. Experimental results using real-world clickstream data demonstrate that our method achieves higher prediction performance than that of a state-of-the-art optimization model and common machine learning methods.
△ Less
Submitted 24 May, 2020; v1 submitted 18 April, 2020;
originally announced April 2020.
-
r-Process enrichment in the Galactic halo characterized by nucleosynthesis variation in the ejecta of coalescing neutron star binaries
Authors:
Takuji Tsujimoto,
Nobuya Nishimura,
Koutarou Kyutoku
Abstract:
A large star-to-star variation in the abundances of r-process elements, as seen in the [Eu/Fe] ratio for Galactic halo stars, is a prominent feature that is distinguishable from other heavy elements. It is, in part, caused by the presence of highly r-process enriched stars, classified as r-II stars ([Eu/Fe]>+1). In parallel, halo stars show that the ratio of a light r-process element (Y) to Eu is…
▽ More
A large star-to-star variation in the abundances of r-process elements, as seen in the [Eu/Fe] ratio for Galactic halo stars, is a prominent feature that is distinguishable from other heavy elements. It is, in part, caused by the presence of highly r-process enriched stars, classified as r-II stars ([Eu/Fe]>+1). In parallel, halo stars show that the ratio of a light r-process element (Y) to Eu is tightly correlated with [Eu/Fe], giving the lowest [Y/Eu] ratio that levels off at r-II stars. On the other hand, recent hydrodynamical simulations of coalescing double neutron stars (cNSNSs) have suggested that r-process sites may be separated into two classes providing different electron-fraction distributions: tidally-driven dynamical ejecta and (dynamical or postmerger) non-tidal ejecta. Here, we show that a widely spanning feature of [Eu/Fe] can be reproduced by models that consider the different masses of tidally-driven dynamical ejecta from both cNSNSs and coalescing black hole/neutron star binaries (cBHNSs). In addition, the observed [Y/Eu] trend is explained by the combined nucleosynthesis in two kinds of ejecta with varying mass asymmetry in double NS systems. Our scenario suggests that massive tidally-driven dynamical ejecta accompanied by massive non-tidal part from cNSNSs or cBHNSs could alone accommodate r-II abundances, including an actinide boost in some cases. The event rate for cNSNSs estimated from our study agrees with the latest result of ~1000 (90% confidence interval of 110-3840) Gpc$^{-3}$yr$^{-1}$ by gravitational-wave detection, and a few events per Gpc$^3$ per year of cBHNSs associated with r-process production are predicted to emerge.
△ Less
Submitted 23 December, 2019;
originally announced December 2019.
-
Impact of Uncertainties in Astrophysical Reaction Rates on Nucleosynthesis in the $νp$ Process
Authors:
T. Rauscher,
N. Nishimura,
G. Cescutti,
R. Hirschi,
A. St. J. Murphy,
C. Fröhlich
Abstract:
The $νp$ process appears in proton-rich, hot matter which is expanding in a neutrino wind and may be realised in explosive environments such as core-collapse supernovae or in outflows from accretion disks. The impact of uncertainties in nuclear reaction cross sections on the finally produced abundances has been studied by applying Monte Carlo variation of all astrophysical reaction rates in a larg…
▽ More
The $νp$ process appears in proton-rich, hot matter which is expanding in a neutrino wind and may be realised in explosive environments such as core-collapse supernovae or in outflows from accretion disks. The impact of uncertainties in nuclear reaction cross sections on the finally produced abundances has been studied by applying Monte Carlo variation of all astrophysical reaction rates in a large reaction network. As the detailed astrophysical conditions of the $νp$ process still are unknown, a parameter study was performed, with 23 trajectories covering a large range of entropies and $Y_\mathrm{e}$. The resulting abundance uncertainties are given for each trajectory. The $νp$ process has been speculated to contribute to the light $p$ nuclides but it was not possible so far to reproduce the solar isotope ratios. It is found that it is possible to reproduce the solar $^{92}$Mo/$^{94}$Mo abundance ratio within nuclear uncertainties, even within a single trajectory. The solar values of the abundances in the Kr-Sr region relative to the Mo region, however, cannot be achieved within a single trajectory. They may still be obtained from a weighted superposition of different trajectories, though, depending on the actual conditions in the production site. For a stronger constraint of the required conditions, it would be necessary to reduce the uncertainties in the 3$α$ and $^{56}$Ni(n,p)$^{56}$Co rates at temperatures $T>3$ GK.
△ Less
Submitted 12 September, 2019; v1 submitted 7 September, 2019;
originally announced September 2019.
-
Uncertainties in $ν$p-process nucleosynthesis from Monte Carlo variation of reaction rates
Authors:
N. Nishimura,
T. Rauscher,
R. Hirschi,
G. Cescutti,
A. St. J. Murphy,
C. Fröhlich
Abstract:
It has been suggested that a $ν$p process can occur when hot, dense, and proton-rich matter is expanding within a strong flux of anti-neutrinos. In such an environment, proton-rich nuclides can be produced in sequences of proton captures and (n,p) reactions, where the free neutrons are created in situ by $\overlineν_\mathrm{e}+\mathrm{p} \rightarrow \mathrm{n}+\mathrm{e}^+$ reactions. The detailed…
▽ More
It has been suggested that a $ν$p process can occur when hot, dense, and proton-rich matter is expanding within a strong flux of anti-neutrinos. In such an environment, proton-rich nuclides can be produced in sequences of proton captures and (n,p) reactions, where the free neutrons are created in situ by $\overlineν_\mathrm{e}+\mathrm{p} \rightarrow \mathrm{n}+\mathrm{e}^+$ reactions. The detailed hydrodynamic evolution determines where the nucleosynthesis path turns off from N = Z line and how far up the nuclear chart it runs. In this work, the uncertainties on the final isotopic abundances stemming from uncertainties in the nuclear reaction rates were investigated in a large-scale Monte Carlo approach, simultaneously varying ten thousand reactions. A large range of model conditions was investigated because a definitive astrophysical site for the $ν$p process has not yet been identified. The present parameter study provides, for each model, identification of the key nuclear reactions dominating the uncertainty for a given nuclide abundance. As all rates appearing in the $ν$p process involve unstable nuclei, and thus only theoretical rates are available, the final abundance uncertainties are larger than those for nucleosynthesis processes closer to stability. Nevertheless, most uncertainties remain below a factor of three in trajectories with robust nucleosynthesis. More extreme conditions allow production of heavier nuclides but show larger uncertainties because of the accumulation of the uncertainties in many rates and because the termination of nucleosynthesis is not at equilibrium conditions. It is also found that the solar ratio of the abundances of ${}^{92}$Mo and ${}^{94}$Mo could be reproduced within uncertainties.
△ Less
Submitted 30 July, 2019;
originally announced July 2019.
-
Impact of Uncertainties in Nuclear Reaction Cross Sections on p-Nucleosynthesis in Thermonuclear Supernovae
Authors:
T. Rauscher,
N. Nishimura,
G. Cescutti,
R. Hirschi,
A. St. J. Murphy,
C. Travaglio
Abstract:
The propagation of uncertainties in reaction cross sections and rates of neutron-, proton-, and alpha-induced reactions into the final isotopic abundances obtained in nucleosynthesis models is an important issue in studies of nucleosynthesis and Galactic Chemical Evolution. We developed a Monte Carlo method to allow large-scale postprocessing studies of the impact of nuclear uncertainties on nucle…
▽ More
The propagation of uncertainties in reaction cross sections and rates of neutron-, proton-, and alpha-induced reactions into the final isotopic abundances obtained in nucleosynthesis models is an important issue in studies of nucleosynthesis and Galactic Chemical Evolution. We developed a Monte Carlo method to allow large-scale postprocessing studies of the impact of nuclear uncertainties on nucleosynthesis. Temperature-dependent rate uncertainties combining realistic experimental and theoretical uncertainties are used. From detailed statistical analyses uncertainties in the final abundances are derived as probability density distributions. Furthermore, based on rate and abundance correlations an automated procedure identifies the most important reactions in complex flow patterns from superposition of many zones or tracers. The method so far was already applied to a number of nucleosynthesis processes. Here we focus on the production of p-nuclei in white dwarfs exploding as thermonuclear (type Ia) supernovae. We find generally small uncertainties in the final abundances despite of the dominance of theoretical nuclear uncertainties. A separate analysis of low- and high-density regions indicates that the total uncertainties are dominated by the high-density regions.
△ Less
Submitted 16 August, 2019; v1 submitted 22 July, 2019;
originally announced July 2019.
-
Item Listing Optimization for E-commerce Websites based on Diversity
Authors:
Naoki Nishimura,
Kotaro Tanahashi,
Koji Suganuma,
Masamichi J. Miyama,
Masayuki Ohzeki
Abstract:
For e-commerce websites, deciding the manner in which items are listed on webpages is an important issue because it can dramatically affect item sales. One of the simplest strategies of listing items to improve the overall sales is to do so in a descending order of sales or sales numbers. However, in lists generated using this strategy, items with high similarity are often placed consecutively. In…
▽ More
For e-commerce websites, deciding the manner in which items are listed on webpages is an important issue because it can dramatically affect item sales. One of the simplest strategies of listing items to improve the overall sales is to do so in a descending order of sales or sales numbers. However, in lists generated using this strategy, items with high similarity are often placed consecutively. In other words, the generated item list might be biased toward a specific preference. Therefore, this study employs penalties for items with high similarity being placed next to each other in the list and transforms the item listing problem to a quadratic assignment problem (QAP). The QAP is well-known as an NP-hard problem that cannot be solved in polynomial time. To solve the QAP, we employ quantum annealing (QA), which exploits the quantum tunneling effect to efficiently solve an optimization problem. In addition, we propose a problem decomposition method based on the structure of the item listing problem because the quantum annealer we use (i.e., D-Wave 2000Q) has a limited number of quantum bits. Our experimental results indicate that we can create an item list that considers both sales and diversity. In addition, we observe that using the problem decomposition method based on a problem structure can lead to a better solution with the quantum annealer in comparison with the existing problem decomposition method.
△ Less
Submitted 27 March, 2019;
originally announced March 2019.
-
Boundary integral equation methods for the two dimensional wave equation in time domain revisited
Authors:
Mio Fukuhara,
Ryota Misawa,
Kazuki Niino,
Naoshi Nishimura
Abstract:
This study considers the stability of time domain BIEMs for the wave equation in 2D. We show that the stability of time domain BIEMs is reduced to a nonlinear eigenvalue problem related to frequency domain integral equations. We propose to solve this non-linear eigenvalue problem numerically with the Sakurai-Sugiura method. After validating this approach numerically in the exterior Dirichlet probl…
▽ More
This study considers the stability of time domain BIEMs for the wave equation in 2D. We show that the stability of time domain BIEMs is reduced to a nonlinear eigenvalue problem related to frequency domain integral equations. We propose to solve this non-linear eigenvalue problem numerically with the Sakurai-Sugiura method. After validating this approach numerically in the exterior Dirichlet problem, we proceed to transmission problems in which we find that some time domain counterparts of "resonance-free" integral equations in frequency domain lead to instability. We finally show how to reformulate these equations to obtain stable numerical schemes.
△ Less
Submitted 18 February, 2019;
originally announced February 2019.
-
The s-process nucleosynthesis in low mass stars: impact of the uncertainties in the nuclear physics determined by Monte Carlo variations
Authors:
G. Cescutti,
R. Hirschi,
N. Nishimura,
T. Rauscher,
J. den Hartogh,
A. St. J. Murphy,
S. Cristallo
Abstract:
We investigated the impact of uncertainties in neutron-capture and weak reactions (on heavy elements) on the s-process nucleosynthesis in low-mass stars using a Monte-Carlo based approach. We performed extensive nuclear reaction network calculations that include newly evaluated temperature-dependent upper and lower limits for the individual reaction rates. Our sophisticated approach is able to eva…
▽ More
We investigated the impact of uncertainties in neutron-capture and weak reactions (on heavy elements) on the s-process nucleosynthesis in low-mass stars using a Monte-Carlo based approach. We performed extensive nuclear reaction network calculations that include newly evaluated temperature-dependent upper and lower limits for the individual reaction rates. Our sophisticated approach is able to evaluate the reactions that impact more significantly the final abundances. We found that beta-decay rate uncertainties affect typically nuclides near s-process branchings, whereas most of the uncertainty in the final abundances is caused by uncertainties in neutron capture rates, either directly producing or destroying the nuclide of interest. Combined total nuclear uncertainties due to reactions on heavy elements are approximately 50%.
△ Less
Submitted 10 October, 2018;
originally announced October 2018.
-
Early chemical evolution of Zn driven by magnetorotational supernovae and the pathway to the solar Zn composition
Authors:
Takuji Tsujimoto,
Nobuya Nishimura
Abstract:
The site of Zn production remains an elusive and challenging problem in astrophysics. A large enhancement of the [Zn/Fe] ratios of very metal-poor stars in the Galactic halo suggests the death of short-lived massive stars, i.e., core-collapse supernovae (CCSNe), as one major site for Zn production. Previous studies have claimed that some specific CCSNe can produce Zn in sufficient quantities. Howe…
▽ More
The site of Zn production remains an elusive and challenging problem in astrophysics. A large enhancement of the [Zn/Fe] ratios of very metal-poor stars in the Galactic halo suggests the death of short-lived massive stars, i.e., core-collapse supernovae (CCSNe), as one major site for Zn production. Previous studies have claimed that some specific CCSNe can produce Zn in sufficient quantities. However, it remains unclear which models can withstand the critical test of observations. Using a Zn abundance feature similar to that of r-process elements in faint satellite galaxies, we find evidence that Zn production took place through much rarer events than canonical CCSNe. This finding can be unified with the implied decrease in the rate of Zn production with an increasing metallicity for Galactic halo stars, which narrows down the major site of Zn production in the early galaxy to magneto-rotational SNe (MR-SNe). On the other hand, in the later phase of galactic evolution, we predict that the major Zn-production site switched from MR-SNe to thermonuclear SNe (SNe Ia). According to this scenario, an accumulation of the contributions from two types of SNe eventually led to the solar isotope composition of Zn which mainly owes 66,68Zn to MR-SNe and 64Zn to SNe Ia triggered by He-detonation. The requirement of Zn production in SNe Ia sheds a new light on the hot debate on the scenario for SN Ia progenitors, suggesting that a He-detonation model might be one major channel for SNe Ia.
△ Less
Submitted 7 August, 2018;
originally announced August 2018.
-
Uncertainties in the production of p nuclides in SN Ia determined by Monte Carlo variations
Authors:
T. Rauscher,
N. Nishimura,
R. Hirschi,
A. St. J. Murphy,
G. Cescutti,
C. Travaglio
Abstract:
Several thousand tracers from a 2D model of a thermonuclear supernova were used in a Monte Carlo post-processing approach to determine p-nuclide abundance uncertainties originating from nuclear physics uncertainties in the reaction rates.
Several thousand tracers from a 2D model of a thermonuclear supernova were used in a Monte Carlo post-processing approach to determine p-nuclide abundance uncertainties originating from nuclear physics uncertainties in the reaction rates.
△ Less
Submitted 28 July, 2018;
originally announced July 2018.
-
Uncertainties in s-process nucleosynthesis in low mass stars determined from Monte Carlo variations
Authors:
G. Cescutti,
R. Hirschi,
N. Nishimura,
J. W. den Hartogh,
T. Rauscher,
A. St. J. Murphy,
S. Cristallo
Abstract:
The main s-process taking place in low mass stars produces about half of the elements heavier than iron. It is therefore very important to determine the importance and impact of nuclear physics uncertainties on this process. We have performed extensive nuclear reaction network calculations using individual and temperature-dependent uncertainties for reactions involving elements heavier than iron,…
▽ More
The main s-process taking place in low mass stars produces about half of the elements heavier than iron. It is therefore very important to determine the importance and impact of nuclear physics uncertainties on this process. We have performed extensive nuclear reaction network calculations using individual and temperature-dependent uncertainties for reactions involving elements heavier than iron, within a Monte Carlo framework. Using this technique, we determined the uncertainty in the main s-process abundance predictions due to nuclear uncertainties link to weak interactions and neutron captures on elements heavier than iron. We also identified the key nuclear reactions dominating these uncertainties. We found that $β$-decay rate uncertainties affect only a few nuclides near s-process branchings, whereas most of the uncertainty in the final abundances is caused by uncertainties in neutron capture rates, either directly producing or destroying the nuclide of interest. Combined total nuclear uncertainties due to reactions on heavy elements are in general small (less than 50%). Three key reactions, nevertheless, stand out because they significantly affect the uncertainties of a large number of nuclides. These are $^{56}$Fe(n,$γ$), $^{64}$Ni(n,$γ$), and $^{138}$Ba(n,$γ$). We discuss the prospect of reducing uncertainties in the key reactions identified in this study with future experiments.
△ Less
Submitted 3 May, 2018;
originally announced May 2018.
-
Incremental Optimization of Independent Sets under Reachability Constraints
Authors:
Takehiro Ito,
Haruka Mizuta,
Naomi Nishimura,
Akira Suzuki
Abstract:
We introduce a new framework for reconfiguration problems, and apply it to independent sets as the first example. Suppose that we are given an independent set $I_0$ of a graph $G$, and an integer $l \ge 0$ which represents a lower bound on the size of any independent set of $G$. Then, we are asked to find an independent set of $G$ having the maximum size among independent sets that are reachable f…
▽ More
We introduce a new framework for reconfiguration problems, and apply it to independent sets as the first example. Suppose that we are given an independent set $I_0$ of a graph $G$, and an integer $l \ge 0$ which represents a lower bound on the size of any independent set of $G$. Then, we are asked to find an independent set of $G$ having the maximum size among independent sets that are reachable from $I_0$ by either adding or removing a single vertex at a time such that all intermediate independent sets are of size at least $l$. We show that this problem is PSPACE-hard even for bounded pathwidth graphs, and remains NP-hard for planar graphs. On the other hand, we give a linear-time algorithm to solve the problem for chordal graphs. We also study the fixed-parameter (in)tractability of the problem with respect to the following three parameters: the degeneracy $d$ of an input graph, a lower bound $l$ on the size of the independent sets, and a lower bound $s$ on the solution size. We show that the problem is fixed-parameter intractable when only one of $d$, $l$, and $s$ is taken as a parameter. On the other hand, we give a fixed-parameter algorithm when parameterized by $s+d$; this result implies that the problem parameterized only by $s$ is fixed-parameter tractable for planar graphs, and for bounded treewidth graphs.
△ Less
Submitted 25 April, 2018;
originally announced April 2018.
-
Reconfiguration of graph minors
Authors:
Benjamin Moore,
Naomi Nishimura,
Vijay Subramanya
Abstract:
Under the reconfiguration framework, we consider the various ways that a target graph $H$ is a {\em minor} of a host graph $G$, where a subgraph of $G$ can be transformed into $H$ by means of {\em edge contraction} (replacement of both endpoints of an edge by a new vertex adjacent to any vertex adjacent to either endpoint). Equivalently, an {\em $H$-model} of $G$ is a labeling of the vertices of…
▽ More
Under the reconfiguration framework, we consider the various ways that a target graph $H$ is a {\em minor} of a host graph $G$, where a subgraph of $G$ can be transformed into $H$ by means of {\em edge contraction} (replacement of both endpoints of an edge by a new vertex adjacent to any vertex adjacent to either endpoint). Equivalently, an {\em $H$-model} of $G$ is a labeling of the vertices of $G$ with the vertices of $H$, where the contraction of all edges between identically-labeled vertices results in a graph containing representations of all edges in $H$.
We explore the properties of $G$ and $H$ that result in a connected {\em reconfiguration graph}, in which nodes represent $H$-models and two nodes are adjacent if their corresponding $H$-models differ by the label of a single vertex of $G$. Various operations on $G$ or $H$ are shown to preserve connectivity. In addition, we demonstrate properties of graphs $G$ that result in connectivity for the target graphs $K_2$, $K_3$, and $K_4$, including a full characterization of graphs $G$ that result in connectivity for $K_2$-models, as well as the relationship between connectivity of $G$ and other $H$-models.
△ Less
Submitted 24 April, 2018;
originally announced April 2018.
-
Reconfiguring spanning and induced subgraphs
Authors:
Tesshu Hanaka,
Takehiro Ito,
Haruka Mizuta,
Benjamin Moore,
Naomi Nishimura,
Vijay Subramanya,
Akira Suzuki,
Krishna Vaidyanathan
Abstract:
Subgraph reconfiguration is a family of problems focusing on the reachability of the solution space in which feasible solutions are subgraphs, represented either as sets of vertices or sets of edges, satisfying a prescribed graph structure property. Although there has been previous work that can be categorized as subgraph reconfiguration, most of the related results appear under the name of the pr…
▽ More
Subgraph reconfiguration is a family of problems focusing on the reachability of the solution space in which feasible solutions are subgraphs, represented either as sets of vertices or sets of edges, satisfying a prescribed graph structure property. Although there has been previous work that can be categorized as subgraph reconfiguration, most of the related results appear under the name of the property under consideration; for example, independent set, clique, and matching. In this paper, we systematically clarify the complexity status of subgraph reconfiguration with respect to graph structure properties.
△ Less
Submitted 16 March, 2018;
originally announced March 2018.
-
Sensitivity to neutron captures and beta-decays of the enhanced s-process in rotating massive stars at low metallicities
Authors:
N. Nishimura,
R. Hirschi,
T. Rauscher
Abstract:
The s-process in massive stars, producing nuclei up to $A\approx 90$, has a different behaviour at low metallicity if stellar rotation is significant. This enhanced s-process is distinct from the s-process in massive stars around solar metallicity, and details of the nucleosynthesis are poorly known. We investigated nuclear physics uncertainties in the enhanced s-process in metal-poor stars within…
▽ More
The s-process in massive stars, producing nuclei up to $A\approx 90$, has a different behaviour at low metallicity if stellar rotation is significant. This enhanced s-process is distinct from the s-process in massive stars around solar metallicity, and details of the nucleosynthesis are poorly known. We investigated nuclear physics uncertainties in the enhanced s-process in metal-poor stars within a Monte-Carlo framework. We applied temperature-dependent uncertainties of reaction rates, distinguishing contributions from the ground state and from excited states. We found that the final abundance of several isotopes shows uncertainties larger than a factor of 2, mostly due to the neutron capture uncertainties. A few nuclei around branching points are affected by uncertainties in the $β$-decay.
△ Less
Submitted 15 February, 2018;
originally announced February 2018.
-
Impacts of nuclear-physics uncertainties in the s-process determined by Monte-Carlo variations
Authors:
N. Nishimura,
G. Cescutti,
R. Hirschi,
T. Rauscher,
J. den Hartogh,
A. St. J. Murphy
Abstract:
The s-process, a production mechanism based on slow-neutron capture during stellar evolution, is the origin of about half the elements heavier than iron. Abundance predictions for s-process nucleosynthesis depend strongly on the relevant neutron-capture and $β$-decay rates, as well as on the details of the stellar model being considered. Here, we have used a Monte-Carlo approach to evaluate the nu…
▽ More
The s-process, a production mechanism based on slow-neutron capture during stellar evolution, is the origin of about half the elements heavier than iron. Abundance predictions for s-process nucleosynthesis depend strongly on the relevant neutron-capture and $β$-decay rates, as well as on the details of the stellar model being considered. Here, we have used a Monte-Carlo approach to evaluate the nuclear uncertainty in s-process nucleosynthesis. We considered the helium burning of massive stars for the weak s-process and low-mass asymptotic-giant-branch stars for the main s-process. Our calculations include a realistic and general prescription for the temperature dependent uncertainty for the reaction cross sections. We find that the adopted uncertainty for (${\rm n},γ$) rates, tens of per cent on average, effects the production of s-process nuclei along the line of $β$-stability, and that the uncertainties in $β$-decay from excited state contributions, has the strongest impact on branching points.
△ Less
Submitted 15 February, 2018;
originally announced February 2018.
-
Deep convolutional neural networks for segmenting 3D in vivo multiphoton images of vasculature in Alzheimer disease mouse models
Authors:
Mohammad Haft-Javaherian,
Linjing Fang,
Victorine Muse,
Chris B. Schaffer,
Nozomi Nishimura,
Mert R. Sabuncu
Abstract:
The health and function of tissue rely on its vasculature network to provide reliable blood perfusion. Volumetric imaging approaches, such as multiphoton microscopy, are able to generate detailed 3D images of blood vessels that could contribute to our understanding of the role of vascular structure in normal physiology and in disease mechanisms. The segmentation of vessels, a core image analysis p…
▽ More
The health and function of tissue rely on its vasculature network to provide reliable blood perfusion. Volumetric imaging approaches, such as multiphoton microscopy, are able to generate detailed 3D images of blood vessels that could contribute to our understanding of the role of vascular structure in normal physiology and in disease mechanisms. The segmentation of vessels, a core image analysis problem, is a bottleneck that has prevented the systematic comparison of 3D vascular architecture across experimental populations. We explored the use of convolutional neural networks to segment 3D vessels within volumetric in vivo images acquired by multiphoton microscopy. We evaluated different network architectures and machine learning techniques in the context of this segmentation problem. We show that our optimized convolutional neural network architecture, which we call DeepVess, yielded a segmentation accuracy that was better than both the current state-of-the-art and a trained human annotator, while also being orders of magnitude faster. To explore the effects of aging and Alzheimer's disease on capillaries, we applied DeepVess to 3D images of cortical blood vessels in young and old mouse models of Alzheimer's disease and wild type littermates. We found little difference in the distribution of capillary diameter or tortuosity between these groups, but did note a decrease in the number of longer capillary segments ($>75μm$) in aged animals as compared to young, in both wild type and Alzheimer's disease mouse models.
△ Less
Submitted 25 February, 2019; v1 submitted 2 January, 2018;
originally announced January 2018.
-
Uncertainties in the production of p nuclides in thermonuclear supernovae determined by Monte Carlo variations
Authors:
N. Nishimura,
T. Rauscher,
R. Hirschi,
A. St. J. Murphy,
G. Cescutti,
C. Travaglio
Abstract:
Thermonuclear supernovae originating from the explosion of a white dwarf accreting mass from a companion star have been suggested as a site for the production of $p$ nuclides. Such nuclei are produced during the explosion, in layers enriched with seed nuclei coming from prior strong $s$ processing. These seeds are transformed to proton-richer isotopes mainly by photodisintegration reactions. Sever…
▽ More
Thermonuclear supernovae originating from the explosion of a white dwarf accreting mass from a companion star have been suggested as a site for the production of $p$ nuclides. Such nuclei are produced during the explosion, in layers enriched with seed nuclei coming from prior strong $s$ processing. These seeds are transformed to proton-richer isotopes mainly by photodisintegration reactions. Several thousand trajectories from a 2D explosion model were used in a Monte Carlo approach. Temperature-dependent uncertainties were assigned individually to thousands of rates varied simultaneously in post-processing in an extended nuclear reaction network. The uncertainties in the final nuclear abundances originating from uncertainties in the astrophysical reaction rates were determined. In addition to the 35 classical $p$ nuclides, abundance uncertainties were also determined for the radioactive nuclides $^{92}$Nb, $^{97,98}$Tc, $^{146}$Sm, and for the abundance ratios $Y$(${}^{92}$Mo)/$Y$(${}^{94}$Mo), $Y$(${}^{92}$Nb)/$Y$(${}^{92}$Mo), $Y$(${}^{97}$Tc)/$Y$(${}^{98}$Ru), $Y$(${}^{98}$Tc)/$Y$(${}^{98}$Ru), and $Y$(${}^{146}$Sm)/$Y$(${}^{144}$Sm), important for Galactic Chemical Evolution studies. Uncertainties found were generally lower than a factor of two, although most nucleosynthesis flows mainly involve predicted rates with larger uncertainties. The main contribution to the total uncertainties comes from a group of trajectories with high peak density originating from the interior of the exploding white dwarf. The distinction between low-density and high-density trajectories allows more general conclusions to be drawn, also applicable to other simulations of white dwarf explosions.
△ Less
Submitted 11 December, 2017; v1 submitted 24 November, 2017;
originally announced November 2017.
-
Mass Ejection from the Remnant of a Binary Neutron Star Merger: Viscous-Radiation Hydrodynamics Study
Authors:
Sho Fujibayashi,
Kenta Kiuchi,
Nobuya Nishimura,
Yuichiro Sekiguchi,
Masaru Shibata
Abstract:
We perform long-term general relativistic neutrino radiation hydrodynamics simulations (in axisymmetry) for a massive neutron star (MNS) surrounded by a torus, which is a canonical remnant formed after the binary neutron star merger. We take into account the effects of viscosity, which is likely to arise in the merger remnant due to magnetohydrodynamical turbulence. As the initial condition, we em…
▽ More
We perform long-term general relativistic neutrino radiation hydrodynamics simulations (in axisymmetry) for a massive neutron star (MNS) surrounded by a torus, which is a canonical remnant formed after the binary neutron star merger. We take into account the effects of viscosity, which is likely to arise in the merger remnant due to magnetohydrodynamical turbulence. As the initial condition, we employ the azimuthally averaged data of the MNS-torus system derived in a three-dimensional, numerical-relativity simulation for the binary neutron star merger. The viscous effect plays key roles for the remnant evolution and mass ejection from it in two phases of the evolution. In the first $t\lesssim10$ ms, a differential rotation state of the MNS is changed to a rigidly rotating state, and as a result, a sound wave, which subsequently becomes a shock wave, is formed in the vicinity of the MNS due to the variation of the quasi-equilibrium state of the MNS. The shock wave induces significant mass ejection of mass $\sim(0.5-2.0)\times 10^{-2}M_\odot$ for the alpha viscosity parameter of $0.01-0.04$. For the longer-term evolution with $\sim 0.1-10$ s, a significant fraction of the torus material is ejected. The ejecta mass is likely to be of order $10^{-2}M_\odot$, so that the total mass of the viscosity-driven ejecta could dominate that of the dynamical ejecta of mass $\lesssim 10^{-2}M_\odot$. The electron fraction, $Y_e$, of the ejecta is always high enough ($Y_e\gtrsim0.25$) that this post-merger ejecta is lanthanide-poor; hence, the opacity of the ejecta is likely to be $\sim 10-100$ times lower than that of the dynamical ejecta. This indicates that the electromagnetic signal from the ejecta would be rapidly evolving, bright, and blue if it is observed from a small viewing angle ($\lesssim 45^\circ$) for which the effect of the dynamical ejecta is minor.
△ Less
Submitted 5 June, 2018; v1 submitted 6 November, 2017;
originally announced November 2017.
-
Abundance Uncertainties Obtained With the PizBuin Framework For Monte Carlo Reaction Rate Variations
Authors:
T. Rauscher,
N. Nishimura,
G. Cescutti,
R. Hirschi,
A. St. J. Murphy
Abstract:
Uncertainties in nucleosynthesis models originating from uncertainties in astrophysical reaction rates were estimated in a Monte Carlo variation procedure. Thousands of rates were simultaneously varied within individual, temperature-dependent errors to calculate their combined effect on final abundances. After a presentation of the method, results from application to three different nucleosynthesi…
▽ More
Uncertainties in nucleosynthesis models originating from uncertainties in astrophysical reaction rates were estimated in a Monte Carlo variation procedure. Thousands of rates were simultaneously varied within individual, temperature-dependent errors to calculate their combined effect on final abundances. After a presentation of the method, results from application to three different nucleosynthesis processes are shown: the $γ$-process and the s-process in massive stars, and the main s-process in AGB stars (preliminary results). Thermal excitation of nuclei in the stellar plasma and the combined action of several reactions increase the final uncertainties above the level of the experimental errors. The total uncertainty, on the other hand, remains within a factor of two even in processes involving a large number of unmeasured rates, with some notable exceptions for nuclides whose production is spread over several stellar layers and for s-process branchings.
△ Less
Submitted 3 September, 2017;
originally announced September 2017.
-
The s-process nucleosynthesis: impact of the uncertainties in the nuclear physics determined by Monte Carlo variations
Authors:
G. Cescutti,
N. Nishimura,
R. Hirschi,
T. Rauscher,
J. W. den Hartogh,
A. St. J. Murphy
Abstract:
We investigated the impact of uncertainties in neutron-capture and weak reactions (on heavy elements) on the s-process nucleosynthesis in low-mass stars and massive stars using a Monte-Carlo based approach. We performed extensive nuclear reaction network calculations that include newly evaluated temperature-dependent upper and lower limits for the individual reaction rates. We found beta-decay rat…
▽ More
We investigated the impact of uncertainties in neutron-capture and weak reactions (on heavy elements) on the s-process nucleosynthesis in low-mass stars and massive stars using a Monte-Carlo based approach. We performed extensive nuclear reaction network calculations that include newly evaluated temperature-dependent upper and lower limits for the individual reaction rates. We found beta-decay rate uncertainties affect only a few nuclei near s-process branchings, whereas most of the uncertainty in the final abundances is caused by uncertainties in the neutron capture rates. We suggest a list of uncertain rates as candidates for improved measurement by future experiments.
△ Less
Submitted 2 August, 2017;
originally announced August 2017.
-
Impacts of nuclear-physics uncertainty in stellar temperatures on the s-process nucleosynthesis
Authors:
N. Nishimura,
G. Cescutti,
R. Hirschi,
T. Rauscher,
J. Den Hartogh,
A. St. J. Murphy
Abstract:
We evaluated the uncertainty relevant to s-process nucleosynthesis using a Monte-Carlo centred approach. We are based on a realistic and general prescription of temperature dependent uncertainty for the reactions. We considered massive stars for the weak s-process and AGB stars for the main s-process. We found that the adopted uncertainty for (n,$γ$) rates, tens of per cent on average, affect the…
▽ More
We evaluated the uncertainty relevant to s-process nucleosynthesis using a Monte-Carlo centred approach. We are based on a realistic and general prescription of temperature dependent uncertainty for the reactions. We considered massive stars for the weak s-process and AGB stars for the main s-process. We found that the adopted uncertainty for (n,$γ$) rates, tens of per cent on average, affect the production of s-process nuclei along the $β$-stability line, while for $β$-decay, for which contributions from excited states enhances the uncertainty, has the strongest impact on branching points.
△ Less
Submitted 24 January, 2017;
originally announced January 2017.