The energy levels of the first few low-lying states of helium and lithium atoms in intense magnetic fields up to $\approx 10^8-10^9$~T are calculated in this study. A pseudospectral method is employed for the computational procedure. The methodology involves computing the eigenvalues and eigenvectors of the generalized two-dimensional Hartree-Fock partial differential equations for these two- and three-electron systems in a self-consistent manner. The method exploits the natural symmetries of the problem without assumptions of any basis functions for expressing the wave functions of the electrons or the commonly employed adiabatic approximation. It is seen that the results obtained here for a few of the most tightly bound states of each of the atoms, helium and lithium, are in good agreement with findings elsewhere. In this regard, we report new data for two new states of lithium that have not been studied thus far in the literature. It is also seen that the pseudospectral method employed here is considerably more economical, from a computational point of view, than previously employed methods such as a finite-element based approach. The key enabling advantage of the method described here is the short computational times which are on the order of seconds for obtaining accurate results for heliumlike systems.

We investigate the possibility that a magnetic field may be present in the star $o-$Ceti (hereafter, Mira) and that the field plays a role in the star's mass loss. The model presented here is an application of an earlier derived theory that has been successfully employed for intermediate and high-mass evolved stars, and is now extended to the low-mass end. The modelling shows that it is possible to obtain a hybrid magnetohydrodynamic-dust-driven wind scenario for Mira, in which the role of a magnetic field in the equatorial plane of the star is dynamically important for producing a stellar wind. The wind velocity and the temperatures obtained from the model appear consistent with findings elsewhere.

The energy levels of hydrogen and helium atoms in strong magnetic fields are calculated in this study. The current work contains estimates of the binding energies of the first few low-lying states of these systems that are improvements upon previous estimates. The methodology involves computing the eigenvalues and eigenvectors of the generalized two-dimensional Hartree-Fock partial differential equations for these one- and two-electron systems in a self-consistent manner. The method described herein is applicable to calculations of atomic structure in magnetic fields of arbitrary strength as it exploits the natural symmetries of the problem without assumptions of any basis functions for expressing the wave functions of the electrons or the commonly employed adiabatic approximation. The method is found to be readily extendable to systems with more than two electrons.

We present calculations for a magnetised hybrid wind model for Betelegeuse ($\alpha - $Orionis). The model is a direct application of our previously derived theory, combining a canonical Weber-Davis (WD) stellar wind with dust grains in the envelope of an AGB star \citep[see][]{Thirumalai2010}. The resulting hybrid picture provides a mechanism for solving the problem of lifting stellar material up from the photosphere \citep[e.g.][] {Harper2009,Guandalini2006,Jura1984} and into the circumstellar envelope. It also predicts wind velocities in agreement with current estimates. Our approach reveals that magnetic fields in supergiant stars like Betelgeuse \citep[see][]{Auriere2010}, may play a vital role in determining the nature of the stellar outflow and consequently, opens a new avenue of investigation in the field of hybrid stellar winds.

We present calculations for a magnetised hybrid wind model for Asymptotic Giant Branch (AGB) stars. The model incorporates a canonical Weber-Davis (WD) stellar wind with dust grains in the envelope of an AGB star. The resulting hybrid picture preserves traits of both types of winds. It is seen that this combination requires that the dust-parameter ($\Gamma_{d}$) be less than unity in order to achieve an outflow. The emergence of critical points in the wind changes the nature of the dust-driven outflow, simultaneously, the presence of a dust condensation radius changes the morphology of the magnetohydrodynamic (MHD) solutions for the wind. In this context, we additionally investigate the effect of having magnetic-cold spots on the equator of an AGB star and its implications for dust formation; which are seen to be consistent with previous findings.

We present a new pseudospectral algorithm for the calculation of the structure of atoms in strong magnetic fields. We have verified this technique for one, two and three-electron atoms in zero magnetic fields against laboratory results and find typically better than one-percent accuracy. We further verify this technique against the state-of-the-art calculations of hydrogen, helium and lithium in strong magnetic fields (up to about $2\times 10^{6}$ T) and find a similar level of agreement. The key enabling advantages of the algorithm are its simplicity (about 130 lines of commented code) and its speed (about $10^2-10^5$ times faster than finite-element methods to achieve similar accuracy).