We study the numerical solution of the non-relativistic Schr\"{o}dinger equation for two-electron atoms in ground and excited S-states using pseudospectral (PS) methods of calculation. The calculation achieves convergence rates for the energy, Cauchy error in the wavefunction, and variance in local energy that are exponentially fast for all practical purposes. The method requires three separate subdomains to handle the wavefunction's cusp-like behavior near the two-particle coalescences. The use of three subdomains is essential to maintaining exponential convergence. A comparison of several different treatments of the cusps and the semi-infinite domain suggest that the simplest prescription is sufficient. For many purposes it proves unnecessary to handle the logarithmic behavior near the three-particle coalescence in a special way. The PS method has many virtues: no explicit assumptions need be made about the asymptotic behavior of the wavefunction near cusps or at large distances, the local energy is exactly equal to the calculated global energy at all collocation points, local errors go down everywhere with increasing resolution, the effective basis using Chebyshev polynomials is complete and simple, and the method is easily extensible to other bound states. This study serves as a proof-of-principle of the method for more general two- and possibly three-electron applications.

In dense stellar clusters, binary-single and binary-binary encounters can ultimately lead to collisions involving two or more stars. A comprehensive survey of multi-star collisions would need to explore an enormous amount of parameter space, but here we focus on a number of representative cases involving low-mass main-sequence stars. Using both Smoothed Particle Hydrodynamics (SPH) calculations and a much faster fluid sorting software package (MMAS), we study scenarios in which a newly formed product from an initial collision collides with a third parent star. By varying the order in which the parent stars collide, as well as the orbital parameters of the collision trajectories, we investigate how factors such as shock heating affect the chemical composition and structure profiles of the collision product. Our simulations and models indicate that the distribution of most chemical elements within the final product is not significantly affected by the order in which the stars collide, the direction of approach of the third parent star, or the periastron separations of the collisions. We find that the sizes of the products, and hence their collisional cross sections for subsequent encounters, are sensitive to the order and geometry of the collisions. For the cases that we consider, the radius of the product formed in the first (single-single star) collision ranges anywhere from roughly 2 to 30 times the sum of the radii of its parent stars. The final product formed in our triple-star collisions can easily be as large or larger than a typical red giant. We therefore expect the collisional cross section of a newly formed product to be greatly enhanced over that of a thermally relaxed star of the same mass.