Source author record

Changjian Xie

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

ResearcherUnclaimed source record

math.NA Numerical Analysis physics.comp-ph astro-ph.GA physics.atom-ph physics.chem-ph

Catalog footprint

What is connected

5works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2021arXiv

A second-order numerical method for Landau-Lifshitz-Gilbert equation with large damping parameters

A second order accurate numerical scheme is proposed and implemented for the Landau-Lifshitz-Gilbert equation, which models magnetization dynamics in ferromagnetic materials, with large damping parameters. The main advantages of this method are associated with the following features: (1) It only solves linear systems of equations with constant coefficients where fast solvers are available, so that the numerical efficiency has been greatly improved, in comparison with the existing Gauss-Seidel project method. (2) The second-order accuracy in time is achieved, and it is unconditionally stable for large damping parameters. Moreover, both the second-order accuracy and the great efficiency improvement will be verified by several numerical examples in the 1D and 3D simulations. In the presence of large damping parameters, it is observed that this method is unconditionally stable and finds physically reasonable structures while many existing methods have failed. For the domain wall dynamics, the linear dependence of wall velocity with respect to the damping parameter and the external magnetic field will be obtained through the reported simulations.

preprint2020arXiv

Stereodynamic control of overlapping resonances in cold molecular collisions

Stereodynamic control of resonant molecular collisions has emerged as a new frontier in cold molecule research. Recent experimental studies have focused on weakly interacting molecular systems such as HD collisions with H$_2$, D$_2$ and He. We report here the possibility of such control in strongly interacting systems taking rotational relaxation in cold collisions of HCl and H$_2$. Using explicit quantum scattering calculations in full six dimensions it is shown that robust control of the collision dynamics is possible even when multiple (overlapping) shape-resonances coexist in a narrow energy range, indicating that cold stereochemistry offers great promise for many molecules beyond simple systems. We demonstrate a striking case where two prominent peaks in overlapping resonances are switched-off simultaneously by suitable alignment of the HCl molecule.

preprint2019arXiv

Second-order semi-implicit projection methods for micromagnetics simulations

Micromagnetics simulations require accurate approximation of the magnetization dynamics described by the Landau-Lifshitz-Gilbert equation, which is nonlinear, nonlocal, and has a non-convex constraint, posing interesting challenges in developing numerical methods. In this paper, we propose two second-order semi-implicit projection methods based on the second-order backward differentiation formula and the second-order interpolation formula using the information at previous two temporal steps. Unconditional unique solvability of both methods is proved, with their second-order accuracy verified through numerical examples in both 1D and 3D. The efficiency of both methods is compared to that of another two popular methods. In addition, we test the robustness of both methods for the first benchmark problem with a ferromagnetic thin film material from National Institute of Standards and Technology.

preprint2019arXiv

Two improved Gauss-Seidel projection methods for Landau-Lifshitz-Gilbert equation

In this paper, we present two improved Gauss-Seidel projection methods with unconditional stability. The first method updates the gyromagnetic term and the damping term simultaneously and follows by a projection step. The second method introduces two sets of approximate solutions, where we update the gyromagnetic term and the damping term simultaneously for one set of approximate solutions and apply the projection step to the other set of approximate solutions in an alternating manner. Compared to the original Gauss-Seidel projection method which has to solve heat equations $7$ times at each time step, the improved methods solve heat equations $5$ times and $3$ times, respectively. First-order accuracy in time and second-order accuracy in space are verified by examples in both 1D and 3D. In addition, unconditional stability with respect to both the grid size and the damping parameter is confirmed numerically. Application of both methods to a realistic material is also presented with hysteresis loops and magnetization profiles. Compared with the original method, the recorded running times suggest that savings of both methods are about $2/7$ and $4/7$ for the same accuracy requirement, respectively.

preprint2016arXiv

The ortho-to-para ratio of interstellar NH$_2$: Quasi-classical trajectory calculations and new simulations

Based on recent $Herschel$ results, the ortho-to-para ratio (OPR) of NH$_2$ has been measured towards the following high-mass star-forming regions: W31C (G10.6-0.4), W49N (G43.2-0.1), W51 (G49.5-0.4), and G34.3+0.1. The OPR at thermal equilibrium ranges from the statistical limit of three at high temperatures to infinity as the temperature tends toward zero, unlike the case of H$_{2}$. Depending on the position observed along the lines-of-sight, the OPR was found to lie either slightly below the high temperature limit of three (in the range $2.2-2.9$) or above this limit ($\sim3.5$, $\gtrsim 4.2$, and $\gtrsim 5.0$). In low temperature interstellar gas, where the H$_{2}$ is para-enriched, our nearly pure gas-phase astrochemical models with nuclear-spin chemistry can account for anomalously low observed NH$_2$-OPR values. We have tentatively explained OPR values larger than three by assuming that spin thermalization of NH$_2$ can proceed at least partially by H-atom exchange collisions with atomic hydrogen, thus increasing the OPR with decreasing temperature. In this paper, we present quasi-classical trajectory calculations of the H-exchange reaction NH$_2$ + H, which show the reaction to proceed without a barrier, confirming that the H-exchange will be efficient in the temperature range of interest. With the inclusion of this process, our models suggest both that OPR values below three arise in regions with temperatures $\gtrsim20-25$~K, depending on time, and values above three but lower than the thermal limit arise at still lower temperatures.