mSOUND is an open-source toolbox written in MATLAB. This toolbox is intended for modeling linear/ nonlinear acoustic wave propagation in media (primarily biological tissues) with arbitrary heterogeneities, in which, the speed of sound, density, attenuation coefficient, powerlaw exponent, and nonlinear coefficient are all spatially varying functions. The computational model is an iterative one-way model based on a mixed domain method. In this article, a general guideline is given along with three representative examples to illustrate how to set up simulations using mSOUND. The first example uses the transient mixed-domain method (TMDM) forward projection to compute the transient acoustic field for a given source defined on a plane. The second example uses the frequency-specific mixed-domain method (FSMDM) forward projection to rapidly obtain the pressure distribution directly at the frequencies of interest, assuming linear or weakly nonlinear wave propagation. The third example demonstrates how to use TMDM backward projection to reconstruct the initial acoustic pressure field to facilitate photoacoustic tomography (PAT). mSOUND (ht. tps://m-sound.github.io/mSOUND/home) is designed to be complementary to existing ultrasound modeling toolboxes and is expected to be useful for a wide range of applications in medical ultrasound including treatment planning, PAT, transducer design, and characterization. Index Terms-Frequency-specific mixed-domain method (FSMDM), heterogeneity, high intensity focused ultrasound, mSOUND toolbox, transient mixed-domain method (TMDM). I. INTRODUCTION N UMERICAL modeling of wave propagation provides an important avenue to help understand ultrasoundtissue interaction [1], [2], innovate medical ultrasound technologies [3], and facilitate the characterization and design of ultrasound transducers in an efficient and cost-effective manner [4]-[6].To date, many simulation toolboxes based on various algorithms have been made available in the public domain.