Optical diffraction tomography (ODT) reconstructs a sample's volumetric refractive index (RI) to create high-contrast, quantitative 3D visualizations of biological samples. However, standard implementations of ODT use interferometric systems, and so are sensitive to phase instabilities, complex mechanical design, and coherent noise. Furthermore, their reconstruction framework is typically limited to weaklyscattering samples, and thus excludes a whole class of multiple-scattering samples. Here, we implement a new 3D RI microscopy technique that utilizes a computational multi-slice beam propagation method to invert the optical scattering process and reconstruct high-resolution (NA>1.0) 3D RI distributions of multiple-scattering samples. The method acquires intensity-only measurements from different illumination angles, and then solves a non-linear optimization problem to recover the sample's 3D RI distribution. We experimentally demonstrate reconstruction of samples with varying amounts of multiple scattering: a 3T3 fibroblast cell, a cluster of C. elegans embryos, and a whole C. elegans worm, with lateral and axial resolutions of ≤250 nm and ≤900 nm, respectively. for text or data mining, so long as such uses are for non-commercial purposes and appropriate attribution is maintained. All other rights are reserved. Fluorescent imaging has enabled stunning visualizations of biological processes at a variety of size scales and resolutions, for studies of gene expression, protein interactions, intracellular dynamics, etc [1][2][3][4]. However, the fluorescent techniques require exogenous biological labels, and so do not directly give endogenous information about a sample's biological structure.Optical diffraction tomography (ODT) also targets 3D biological imaging. In contrast to fluorescent methods, ODT avoids the use of exogenous biological labels, and instead utilizes the intrinsic optical variation within a sample to reconstruct its 3D refractive-index (RI) distribution [5][6][7][8][9][10][11]. Hence, ODT avoids some of fluorescent imaging's main drawbacks, such as photobleaching, slow acquisition speed, low signal-to-noise (SNR) ratio, and complex samplepreparation protocol. Furthermore, RI imaging enables examination of the structural, mechanical, and biochemical properties of a sample, which are important for studies in morphology, mass, shear stiffness, and spectroscopy [9,[12][13][14][15].Standard implementations of ODT use either a rotating sample or a scanning laser beam to capture the angle-specific scattering arising from the sample [5,7,[16][17][18]. Under the assumption of weak scattering (i.e., 1st Born or Rytov approximations), 2D electric-field measurements directly yield information about the sample's 3D scattering potential [19][20][21]. Standard ODT reconstruction algorithms utilize the Fourier diffraction theorem to project the information contained in each electric-field measurement onto spherical shells (i.e., Ewald surfaces) in the 3D Fourier space of the sample's scattering potential [22,23]. ...