Source author record

Kotaro Fujisawa

Kotaro Fujisawa appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

15 published item(s)

preprint2022arXiv

A novel formulation for the evolution of relativistic rotating stars

We present a new formulation to construct numerically equilibrium configurations of rotating stars in general relativity. Having in mind the application to their quasi static evolutions, we adopt a Lagrangian formulation of our own devising, in which we solve force balance equations to seek for the positions of fluid elements assigned to the grid points, instead of the ordinary Eulerian formulation. Unlike previous works in the literature, we do not employ the first integral of the Euler equation, which is not obtained by an analytic integration in general. We assign a mass, specific angular momentum and entropy to each fluid element in contrast to the previous methods, in which the spatial distribution of the angular velocity or angular momentum is specified. Those distributions are determined after the positions of all fluid elements (or grid points) are derived in our formulation. We solve the large system of algebraic nonlinear equations that are obtained by discretizing the time-independent Euler and Einstein equations in the finite-elements method by using our new multi-dimensional root-finding scheme, named the W4 method. To demonstrate the capability of our new formulation, we construct some rotational configurations both barotropic and baroclinic. We also solve three evolutionary sequences that mimic the cooling, mass-loss, and mass-accretion as simple toy models.

preprint2022arXiv

Fast optical flares from M dwarfs detected by a one-second-cadence survey with Tomo-e Gozen

We report a one-second-cadence wide-field survey for M-dwarf flares using the Tomo-e Gozen camera mounted on the Kiso Schmidt telescope. We detect 22 flares from M3-M5 dwarfs with rise times and amplitudes ranging from $5\, \mathrm{sec} \lesssim t_\mathrm{rise} \lesssim 100\,\mathrm{sec}$ and $0.5 \lesssim ΔF/F_{\star} \lesssim 20$, respectively. The flare light curves mostly show steeper rises and shallower decays than those obtained from the Kepler one-minute cadence data and tend to have flat peak structures. Assuming a blackbody spectrum with temperatures of $9,000-15,000\,\mathrm{K}$, the peak luminosities and bolometric energies are estimated to be $10^{29}\,\mathrm{erg\,sec^{-1}} \lesssim L_\mathrm{peak} \lesssim 10^{31}\,\mathrm{erg\,sec^{-1}}$ and $10^{31}\,\mathrm{erg} \lesssim E_{\rm bol} \lesssim 10^{34}\,\mathrm{erg}$, which constitutes the bright end of fast optical flares for M dwarfs. We confirm that more than 90\% of the host stars of the detected flares are magnetically active based on their H$α$ emission line intensities obtained by LAMOST. The estimated occurrence rate of the detected flares is $\sim 0.7$ per day per an active star, indicating they are common in magnetically active M dwarfs. We argue that the flare light curves can be explained by the chromospheric compression model; the rise time is broadly consistent with the Alfvén transit time of a magnetic loop with a length scale of $l_\mathrm{loop} \sim 10^4\,\mathrm{km}$ and a field strength of $1,000\,\mathrm{G}$, while the decay time is likely determined by the radiative cooling of the compressed chromosphere down to near the photosphere with a temperature of $\gtrsim 10,000\,\mathrm{K}$. These flares from M dwarfs could be a major contamination source for a future search of fast optical transients of unknown types.

preprint2022arXiv

Magnetic field sustained by the elastic force in neutron star crusts

We investigate the magneto--elastic equilibrium of a neutron star crust and magnetic energy stored by the elastic force. The solenoidal motion driven by the Lorentz force can be controlled by the magnetic elastic force, so that conditions for the magnetic field strength and geometry are less restrictive. For equilibrium models, the minor solenoidal part of the magnetic force is balanced by a weak elastic force because the irrotational part is balanced by the dominant gravity and pressure forces. Therefore, a strong magnetic field may be confined in the interior, regardless of poloidal or toroidal components. We numerically calculated axially symmetric models with the maximum shear--strain, and found that a magnetic energy $> 10^{46}$ erg can be stored in the crust, even for a normal surface dipole-field-strength ($<10^{13}$ G). The magnetic energy much exceeds the elastic energy ($ 10^{44} -10^{45}$ erg). The shear--stress spatial distribution revealed that the elastic structure is likely to break down near the surface. In particular, the critical position is highly localized at a depth less than 100 m from the surface.

preprint2022arXiv

Singularity-Avoiding Multi-Dimensional Root-Finder

We proposed in this paper a new method, which we named the W4 method, to solve nonlinear equation systems. It may be regarded as an extension of the Newton-Raphson~(NR) method to be used when the method fails. Indeed our method can be applied not only to ordinary problems with non-singular Jacobian matrices but also to problems with singular Jacobians, which essentially all previous methods that employ the inversion of the Jacobian matrix have failed to solve. In this article, we demonstrate that (i) our new scheme can define a non-singular iteration map even for those problems by utilizing the singular value decomposition, (ii) a series of vectors in the new iteration map converges to the right solution under a certain condition, (iii) the standard two-dimensional problems in the literature that no single method proposed so far has been able to solve completely are all solved by our new method.

preprint2022arXiv

Strong toroidal magnetic fields sustained by the elastic crust in a neutron star

We investigate new solutions for magnetized neutron stars with a barotropic core in magnetohydrodynamic (MHD) equilibrium and a magneto-elastic crust, which was neglected by previous studies concerning stars in MHD equilibrium. The Lorentz force of the barotropic star is purely irrotational and the structures of magnetic fields are constrained. By contrast, a solenoidal component of the Lorentz force exists in the elastic crust and the structures of the magnetic fields are less restricted. We find that the minor solenoidal component in the elastic crust is important for sustaining the strong magnetic field in the core. Unlike previous studies, the toroidal magnetic field exists in the entire region of the core, and we obtain equilibrium states with large toroidal magnetic fields, where the toroidal magnetic energy is larger than the poloidal magnetic energy. The elastic force of the crust sustains an order of $10^{15}~\mathrm{G}$ toroidal magnetic field in the core, and the maximum strength of the toroidal magnetic field is approximately proportional to the crust thickness.

preprint2019arXiv

The optically thick rotating magnetic wind from a massive white dwarf merger product

WD J005311 is a newly identified white dwarf (WD) in a mid-infrared nebula. The spectroscopic observation indicates the existence of a neon-enriched carbon/oxygen wind with a terminal velocity of $v_{\infty,\rm obs}\sim 16,000\,\rm km\,s^{-1}$ and a mass loss rate of $\dot M_{\rm obs}\sim 3.5\times 10^{-6}\,M_\odot$ yr$^{-1}$. Here we consistently explain the properties of WD J005311 using a newly constructed wind solution, where the optically thick outflow is launched from the carbon burning shell on an oxygen-neon core and accelerated by the rotating magnetic field to become supersonic and unbound well below the photosphere. Our model implies that WD J005311 has a mass of $M_* \sim 1.1\mbox{-}1.3\,M_\odot$, a magnetic field of $B_* \sim (2\mbox{-}5)\times 10^7\,\rm G$, and a spin angular frequency of $Ω\sim 0.2\mbox{-}0.5 \,\rm s^{-1}$. The large magnetic field and fast spin support the carbon-oxygen WD merger origin. WD J005311 will neither explode as a type Ia supernova nor collapse into a neutron star. If the wind continues to blow another few kyr, WD J005311 will spin down significantly and join to the known sequence of slowly-rotating magnetic WDs. Otherwise it may appear as a fast-spinning magnetic WD and could be a new high energy source.

preprint2016arXiv

Hyperbolic Self-Gravity Solver for Large Scale Hydrodynamical Simulations

A new computationally efficient method has been introduced to treat self-gravity in mesh based hydrodynamical simulations. It is applied simply by slightly modifying the Poisson equation into an inhomogeneous wave equation. This roughly corresponds to the weak field limit of the Einstein equations in general relativity, and as long as the gravitation propagation speed is taken to be larger than the hydrodynamical characteristic speed, the results agree with solutions for the Poisson equation. The solutions almost perfectly agree if the domain is taken large enough, or appropriate boundary conditions are given. Our new method can not only significantly reduce the computational time compared with existent methods, but is also fully compatible with massive parallel computation, nested grids and adaptive mesh refinement techniques, all of which can accelerate the progress in computational astrophysics and cosmology.

preprint2016arXiv

Rotational equilibria by Lagrangian variational principle: toward multi-dimensional stellar evolutions

We have developed a new formulation to obtain self-gravitating, axisymmetric configurations in permanent rotation. The formulation is based on the Lagrangian variational principle with a triangulated mesh. It treats not only barotropic but also baroclinic equations of state. We compare the various stellar equilibria obtained by our new scheme with those by Hachisu's self-consistent field scheme for the barotropic case, and those by Fujisawa's self-consistent field scheme for the baroclinic case. Included in these rotational configurations are those with shellular-type rotations, which are commonly assumed in the evolution calculation of rotating stars. Although radiation processes, convections and meridional flows have not been taken into account in this study, we have in mind the application of this method to the two-dimensional evolution calculations of rotating stars, for which the Lagrangian formulation is best suited.

preprint2015arXiv

A versatile numerical method for obtaining structures of rapidly rotating baroclinic stars: self-consistent and systematic solutions with shellular-type rotation

This paper develops a novel numerical method for obtaining structures of rapidly rotating stars based on a self-consistent field scheme. The solution is obtained iteratively. Both rapidly rotating barotropic and baroclinic equilibrium states are calculated self-consistently using this method. Two types of rotating baroclinic stars are investigated by changing the isentropic surfaces inside the star. Solution sequences of these are calculated systematically and critical rotation models beyond which no rotating equilibrium state exists are also obtained. All of these rotating baroclinic stars satisfy necessarily the Bjerknes-Rosseland rules. Self-consistent solutions of baro-clinic stars with shellular-type rotation are successfully obtained where the isentropic surfaces are oblate and the surface temperature is hotter at the poles than at the equator if it is assumed that the star is an ideal gas star. These are the first self-consistent and systematic solutions of rapidly rotating baroclinic stars with shellular-type rotations. Since they satisfy the stability criterion due to their rapid rotation, these rotating baroclinic stars would be dynamically stable. This novel numerical method and the solutions of the rapidly rotating baroclinic stars will be useful for investigating stellar evolution with rapid rotations.

preprint2015arXiv

Appearance of the prolate and the toroidal magnetic field dominated stars - Analytic approach

We have analyzed magnetized equilibrium states and showed a condition for appearance of the prolate and the toroidal magnetic field dominated stars by analytic approaches. Both observations and numerical stability analysis support that the magnetized star would have the prolate and the large internal toroidal magnetic fields. In this context, many investigations concerning magnetized equilibrium states have tried to obtain the prolate and the toroidal dominant solutions, but many of them have failed to obtain such configurations. Since the Lorentz force is a cross product of current density and magnetic field, the prolate shaped configurations and the large toroidal magnetic fields in stars require a special relation between current density and the Lorentz force. We have analyzed simple analytical solutions and found that the prolate and the toroidal dominant configuration require non force-free toroidal current density that flows in the opposite direction with respect to the bulk current within the star. Such current density results in the Lorentz force which makes the stellar shape prolate. Satisfying this special relation between the current density and the Lorentz force is a key for appearance of the prolate and the toroidal magnetic field dominated magnetized star.

preprint2015arXiv

Magnetic field configurations of a magnetar throughout its interior and exterior -- core, crust and magnetosphere

We obtained the magnetic field configurations, including both poloidal and toroidal components, throughout the interior and exterior of magnetars using a realistic equation of state. We divided the magnetized star into the hydromagnetic equilibrium core, Hall equilibrium crust and twisted force-free magnetosphere. We systematically and simultaneously calculated these regions under various boundary conditions using the Green function relaxation method, and noted the following interesting characteristics of these numerical results. First, the strength and structure of core magnetic fields affect the crustal magnetic fields. Second, the current sheet on the core-crust interface affects both internal and external magnetic field configurations. Third, the twisted magnetosphere makes a cross-point of magnetic field lines, such as X-point geometry, in the magnetosphere. The X-point geometry appears and disappears according to the strength of the twisted field in the magnetosphere or the core-crust boundary conditions. Our results mean that both Hall magnetohydrodynamics secular evolution and magnetospheric dynamical evolution are deeply affected by conditions of another region and the core-crust stress of magnetars.

preprint2015arXiv

Magnetized stars with differential rotation and a differential toroidal field

We have succeeded in obtaining magnetized equilibrium states with differential rotation and differential toroidal magnetic fields. If an internal toroidal field of a proto-neutron star is wound up from the initial poloidal magnetic field by differential rotation, the distribution of the toroidal magnetic field is determined by the profile of this differential rotation. However, the distributions of the toroidal fields in all previous magnetized equilibrium studies do not represent the magnetic winding by the differential rotation of the star. In this paper, we investigate a formulation of a differential toroidal magnetic field that represents the magnetic field wound up by differential rotation. We have developed two functional forms of differential toroidal fields which correspond to a v-constant and a j-constant field in analogy to differential rotations. As the degree of the differential becomes very high, the toroidal magnetic field becomes highly localized and concentrated near the rotational axis. Such a differential toroidal magnetic field would suppress the low-T/|W| instability more efficiently even if the total magnetic field energy is much smaller than that of a non-differential toroidal magnetic field.

preprint2014arXiv

Equilibrium solutions of relativistic rotating stars with mixed poloidal and toroidal magnetic fields

Stationary and axisymmetric solutions of relativistic rotating stars with strong mixed poloidal and toroidal magnetic fields are obtained numerically. Because of the mixed components of the magnetic field, the underlying stationary and axisymmetric spacetimes are no longer circular. These configurations are computed from the full set of the Einstein-Maxwell equations, Maxwell's equations and from first integrals and integrability conditions of the magnetohydrodynamic equilibrium equations. After a brief introduction of the formulation of the problem, we present the first results for highly deformed magnetized rotating compact stars.

preprint2013arXiv

Coexistence of oppositely flowing multi-$φ$-currents: Key to large toroidal magnetic fields within stars

We will show the importance of coexistence of oppositely flowing $φ$-currents for magnetized stars to sustain strong toroidal magnetic fields within the stars by an- alyzing stationary states of magnetized stars with surface currents which flow in the opposite direction with respect to the bulk currents within the stars. We have imposed boundary conditions for currents and toroidal magnetic fields to vanish outside the stars. It is important to note that these boundary conditions set an upper limit for the total current within the stars. This upper limit for the total current results in the presence of an upper limit for the magnitude of the energy for the toroidal magnetic fields of the stars. If the stars could have the toroidal surface currents which flow in the opposite directions to the internal toroidal currents, the positively flowing inter- nal toroidal currents can become stronger than the upper limit value of the current for configurations without surface toroidal currents. Thus the energies for the toroidal magnetic fields can become much larger than those for the magnetized stars with- out surface toroidal currents. We have also analyzed the same phenomena appearing in spherical incompressible stars for dipole-like magnetic fields with or without sur- face toroidal currents by employing the zero-flux-boundary method. We have applied those configurations with surface toroidal currents to magnetars and discussed their flares through which magnetic helicities could arise outside the stellar surfaces.

preprint2013arXiv

Counter effects of meridional flows and magnetic fields in stationary axisymmetric self-gravitating barotropes under the ideal MHD approximation: clear examples - toroidal configurations

We obtain the general forms for the current density and the vorticity from the integrability conditions of the basic equations which govern the stationary states of axisymmetric magnetized self-gravitating barotropic objects with meridional flows under the ideal magnetohydrodynamics (MHD) approximation. As seen from the stationary condition equations for such bodies, the presence of the meridional flows and that of the poloidal magnetic fields act oppositely on the internal structures. The different actions of these two physical quantities, the meridional flows and the poloidal magnetic fields, could be clearly seen through stationary structures of the toroidal gaseous configurations around central point masses in the framework of Newtonian gravity because the effects of the two physical quantities can be seen in an amplified way for toroidal systems compared to those for spheroidal stars. The meridional flows make the structures more compact, i.e. the widths of toroids thinner, while the poloidal magnetic fields are apt to elongate the density contours in a certain direction depending on the situation. Therefore, the simultaneous presence of the internal flows and the magnetic fields would work as if there were no such different actions within and around the stationary gaseous objects such as axisymmetric magnetized toroids with internal motions around central compact objects under the ideal MHD approximation, although these two quantities might exist in real systems.