Spherical Harmonics Research Articles

Analysis of the three possible event horizons of anti-de Sitter space–time in the Klein–Gordon equation

In this paper, we calculate exactly the solution of Klein–Gordon equation in anti-de Sitter space–time. Due to spherical symmetry, the angular wave functions are given by spherical harmonics. The radial wave functions are given by local Heun functions. We analyze the behavior of the radial wave function in regions near three possible event horizon positions [Formula: see text] and investigate the decay rate, Hawking radiation and Hawking temperature for these three cases. Comparing the results, we verify that if the decay rate and Hawking radiation are greater, then the closer the event horizon is to the singularity of the black hole. Finally, we analyze the special case of a small black hole limit.

A 3D DEM modeling of soil-rock mixture considering spatial distribution orientation of blocks

Spatial distribution orientations of blocks can cause significant errors in the discrete element model (DEM) calculation of soil-rock mixture (SRM). To avoid this error, spherical harmonic (SH) series whose harmonization degrees fixed at 15 were proposed for block reconstruction. This research refers to the case-history of a deep excavation rift valley spanning from the Mabian to Zhaojue section of the Leshan-Xichang Expressway, mainly containing moderately-weathered silty mudstone, in the Leshan City, Sichuan Province, China. The appropriate degree of finite-term SH series is selected by the volume, surface area. 100 blocks were scanned on site, and sphericity and angularity of the blocks were calculated. The sphericity and angularity of 50 reconstructed blocks were considered for the error analysis of SH method. Moreover, stochastic polyhedron method was considered for comparing different block reconstructions. The maximum block placement angle was defined to control the spatial distribution orientations of the blocks. Large scale direct tests were carried. Numerical simulations of large-scale direct shear tests were conducted to study the influence of the spatial distribution orientations of the blocks on the mechanical properties of the SRMs. The results revealed that the finite-term SH series fixed at 15 accurately reflected the shape characteristics and mechanical behaviors of actual blocks. The spatial distribution orientations of the blocks had a minimal impact on the friction angle and cohesion of SRM constructed through the SH method. The SRMs developed via the SH method exhibited marginal variations in contact force and anisotropy index of contact across diverse block placement strategies. The evolution of coordination number was closer when employing the SH method under varied block placement methods. Blocks reconstructed by the SH method, could mitigate errors in DEM calculation caused by the spatial distribution orientations of the blocks.

FSH3D: 3D Representation via Fibonacci Spherical Harmonics

AbstractSpherical harmonics are a favorable technique for 3D representation, employing a frequency‐based approach through the spherical harmonic transform (SHT). Typically, SHT is performed using equiangular sampling grids. However, these grids are non‐uniform on spherical surfaces and exhibit local anisotropy, a common limitation in existing spherical harmonic decomposition methods. This paper proposes a 3D representation method using Fibonacci Spherical Harmonics (FSH3D). We introduce a spherical Fibonacci grid (SFG), which is more uniform than equiangular grids for SHT in the frequency domain. Our method employs analytical weights for SHT on SFG, effectively assigning sampling errors to spherical harmonic degrees higher than the recovered band‐limited function. This provides a novel solution for spherical harmonic transformation on non‐equiangular grids. The key advantages of our FSH3D method include: 1) With the same number of sampling points, SFG captures more features without bias compared to equiangular grids; 2) The root mean square error of 32‐degree spherical harmonic coefficients is reduced by approximately 34.6% for SFG compared to equiangular grids; and 3) FSH3D offers more stable frequency domain representations, especially for rotating functions. FSH3D enhances the stability of frequency domain representations under rotational transformations. Its application in 3D shape reconstruction and 3D shape classification results in more accurate and robust representations. Our code is publicly available at https://github.com/Miraclelzk/Fibonacci-Spherical-Harmonics.

An optimal ansatz space for moving least squares approximation on spheres

We revisit the moving least squares (MLS) approximation scheme on the sphere Sd-1⊂Rd\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathbb S^{d-1} \\subset {\\mathbb R}^d$$\\end{document}, where d>1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$d>1$$\\end{document}. It is well known that using the spherical harmonics up to degree L∈N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L \\in {\\mathbb N}$$\\end{document} as ansatz space yields for functions in CL+1(Sd-1)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal {C}^{L+1}(\\mathbb S^{d-1})$$\\end{document} the approximation order OhL+1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal {O}\\left( h^{L+1} \\right) $$\\end{document}, where h denotes the fill distance of the sampling nodes. In this paper, we show that the dimension of the ansatz space can be almost halved, by including only spherical harmonics of even or odd degrees up to L, while preserving the same order of approximation. Numerical experiments indicate that using the reduced ansatz space is essential to ensure the numerical stability of the MLS approximation scheme as h→0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$h \\rightarrow 0$$\\end{document}. Finally, we compare our approach with an MLS approximation scheme that uses polynomials on the tangent space of the sphere as ansatz space.

Wigner state and process tomography on near-term quantum devices

We present an experimental scanning-based tomography approach for near-term quantum devices. The underlying method has previously been introduced in an ensemble-based NMR setting. Here we provide a tutorial-style explanation along with suitable software tools to guide experimentalists in its adaptation to near-term pure-state quantum devices. The approach is based on a Wigner-type representation of quantum states and operators. These representations provide a rich visualization of quantum operators using shapes assembled from a linear combination of spherical harmonics. These shapes (called droplets in the following) can be experimentally tomographed by measuring the expectation values of rotated axial tensor operators. We present an experimental framework for implementing the scanning-based tomography technique for circuit-based quantum computers and showcase results from IBM quantum experience. We also present a method for estimating the density and process matrices from experimentally tomographed Wigner functions (droplets). This tomography approach can be directly implemented using the Python-based software package DROPStomo.

Method for prediction magnetic silencing of uncertain energy-saturated extended technical objects in prolate spheroidal coordinate system

Aim. Development of method for prediction by energy-saturated extended technical objects magnetic silencing based on magnetostatics geometric inverse problems solution and magnetic field spatial spheroidal harmonics calculated in prolate spheroidal coordinate system taking into account of technical objects magnetic characteristics uncertainties. Methodology. Spatial prolate spheroidal harmonics of extended technical objects magnetic field model calculated as magnetostatics geometric inverse problems solution in the form of nonlinear minimax optimization problem based on near field measurements for prediction far extended technical objects magnetic field magnitude. Nonlinear objective function calculated as the weighted sum of squared residuals between the measured and predicted magnetic field COMSOL Multiphysics software package used. Nonlinear minimax optimization problems solutions calculated based on particle swarm nonlinear optimization algorithms. Results. Results of prediction extended technical objects far magnetic field magnitude based on extended technical objects spatial prolate spheroidal harmonics of the magnetic field model in the prolate spheroidal coordinate system using near field measurements with consideration of extended technical objects magnetic characteristics uncertainty. Originality. The method for prediction by extended technical objects magnetic cleanliness based on spatial prolate spheroidal harmonics of the magnetic field model in the prolate spheroidal coordinate system with consideration of magnetic characteristics uncertainty is developed. Practical value. The important practical problem of prediction extended technical objects magnetic silencing based on the spatial prolate spheroidal harmonics of the magnetic field model in the prolate spheroidal coordinate system with consideration of extended technical objects magnetic characteristics uncertainty solved. References 48, figures 2.

Evaluation of GOCE/GRACE and combined global geopotential models using GNSS/levelling data over Nigeria

AbstractGlobal Geopotential Models (GGMs) provide valuable information about Earth’s gravity field functionals, such as geoid heights and gravity anomalies. However, ground-based datasets are required to validate these GGMs at the regional and local scales. In this study, we validated the accuracy of GGMs by comparing them with ground-based Global Navigation Satellite System (GNSS)/levelling data for the first time in Nigeria. We employed two validation scenarios: with and without considering spectral consistency using the spectral enhancement method (SEM) to incorporate high and very high frequencies of the gravity field spectrum from the combined global gravity field model (XGM2019e_2159) and the residual terrain model (RTM) derived from the Shuttle Radar Topography Mission (SRTM) data, respectively. The results of this evaluation confirmed that the application of SEM improved the assessment of the GGM solutions in an unbiased manner. Integrating XGM2019e_2159 and SRTM data to constrain the high-frequency component of geoid heights in Gravity Field and Steady-State Ocean Circulation Explorer (GOCE)-based GGMs led to an improvement of approximately 10% in reducing the standard deviation (STD) relative to when SEM was not applied. GO_CONS_GCF_2_TIM_R6 at spherical harmonics (SH) of up to degree and order 260 demonstrated the lowest STD when compared to GO_CONS_GCF_2_DIR_R6 and GO_CONS_GCF_2_SPW_R5, with a reduction from 0.380 m without SEM application to 0.342 m with SEM implementation. In addition, four transformation models, namely, linear, four-parameter, five-parameter, and seven-parameter models, were evaluated. The objective is to mitigate the reference system offsets between the GNSS/levelling data and the GGMs and to identify the particular parametric model with the smallest STD across all GGMs. This effort reduced the GGMs misfits to GNSS/levelling to 0.30 m, representing a 15.3% decrease in STD. Notably, the XGM2019e_2159 model provides this improvement.

Comparison of terrestrial exospheric hydrogen 3D distributions at solar minimum and maximum using TWINS Lyman-α observations

Remote sensing observations of far-ultraviolet (FUV) emissions have been used to estimate 3D neutral hydrogen (H) density models of the terrestrial exosphere under solar minimum (2008) and solar maximum (2013 and 2015) conditions. Specifically, we used Lyman-alpha (Lyman-α, FUV at 121.6 nm) radiance data acquired by the Lyman-alpha detectors (LADs) onboard NASA’s TWINS satellites, collected over several days for each of the 3 years to provide sufficient coverage of the exospheric region. The datasets included Lyman-α measurements taken only above Earth radii (Re) of 3.75, assumed to be the optically thin region of the exosphere, where the measured Lyman-α intensity along a line of sight (LOS) is linearly proportional to the atomic hydrogen column density. Based on the calibration using multiple UV-bright stars, a significant decrease in TWINS1 LAD1/2 sensitivities was found for 2013 (LAD2 ∼1/2, LAD ∼1/3) and 2015 (LAD2 ∼1/3, LAD ∼1/7) compared to 2008. The calibration uncertainty was derived to be, on average, 5%. We estimated the 3D global hydrogen density distributions from these radiance data using two different tomographic inversion methods: first, a parametric fitting method based on a spherical harmonic function of order 3 and second, a high degree-of-freedom retrieval approach to validate our retrievals of H density. Both inversion methods consistently incorporate the effects of Lyman-α absorption within the exosphere and Lyman-α re-emission from Earth’s albedo. Our results reveal that H densities during solar maximum conditions are, on average, 30%–40% higher than those during solar minimum. All models showed a high concentration of atomic H on the dayside and nightside near the Sun–Earth line, which determines a nose/geotail structure consistent with theoretical effects from the solar radiation pressure. Furthermore, we identified that H-density enhancements during solar maximum with respect to solar minimum conditions occur at mid-to-high latitudes, particularly on the dusk side, while no significant enhancement seems to occur near the dayside nose.

Low-degree gravity field coefficients based on inverse GNSS method: insights into hydrological and ice mass change studies

The relative displacements of stations from a global network of Global Navigation Satellite System (GNSS) sites provide information on global mass transport. In this study, we use 19 years of global GNSS station displacements from the 3rd International GNSS Service reprocessing campaign to estimate the coefficients of the spherical harmonics of the Earth’s gravity field up to degree and order 8 using the inverse GNSS method based on elastic loading theory. The results indicate that the C30 coefficient can be derived based on GNSS station displacements as an alternative to solutions provided by Satellite Laser Ranging (SLR) and Gravity Recovery and Climate Experiment (GRACE). GNSS may support GRACE solutions that face the problems of deriving C30, which has fundamental meaning in estimating ice mass changes in polar regions. The recovery of Antarctic ice sheet mass change from January 2007 to December 2020 based on coefficients replaced by GNSS estimates results in a linear trend of − 152 ± 4 Gt/year, in comparison to − 149 ± 2 Gt/year for the replacement based on SLR from GRACE Technical Note #14. The results indicate that the spatial and seasonal patterns of terrestrial water storage changes derived from GNSS are consistent with those estimated using GRACE/GRACE Follow-On and SLR at a few-millimeter level in the Amazon and Brahmaputra River basin regions.

Calculation of the Gaunt Coefficients over Real Spherical Harmonics

Calculation of the Gaunt Coefficients over Real Spherical Harmonics

Analysis on scattering and inner near-field characteristics of a uniaxial anisotropic sphere by an off-axis high-order Bessel (vortex) beam

Based on the generalized Lorenz–Mie theory (GLMT) and the Fourier transform method, a theoretical approach is introduced to study the scattering of a uniaxial anisotropic sphere illuminated by an off-axis high-order Bessel (vortex) beam (HOBVB). According to the orthogonality of the associated Legendre function and exponential function, a concise expression of the expansion coefficients of the off-axis HOBVB in terms of the spherical vector wave functions (SVWFs) is derived that can effectively reconstruct the HOBVB with all conical angles. The differences of scattering characteristics of a uniaxial anisotropic sphere illuminated by an on-axis and off-axis HOBVB and a plane wave are exhibited. Influences of the topological charge, conical angle, particle size, and off-axis distance on the angle distributions of the radar cross-section (RCS), scattering and extinction efficiencies, and asymmetric factor are analyzed in detail. The unique internal and near-field distributions of a uniaxial anisotropic spherical particle illuminated by an on-axis and off-axis HOBVB are demonstrated. The results provide insights into the scattering and Bessel beam–matter interactions and may find important applications in optical propagation and optical micromanipulation, microwave engineering, target shielding, and near-field measurement.

Methodology for 3D sound synthesis of directional acoustic sources by higher-order ambisonics

This paper presents the 3D sound field synthesis using both the multimodal method and higher-order ambisonics (HOA) of the pressure field radiated by directional acoustic sources. Ambisonics is a technique for encoding and reproducing measured or modeled (virtual) sound pressure field, based on a decomposition of the acoustic field over spherical harmonics. The directional source considered in this work is an acoustic horn excited by a flat piston. The free-field radiation from this horn is first modeled accurately over a wide frequency range using the multimodal method, which requires relatively low computational resources. This radiated pressure field, collected on a dual-layer sphere of virtual sensors distributed over a Lebedev geometry, allows its projection into the ambisonic domain. The pressure field is then synthetized in the laboratory's 3D 5th order HOA spatialization sphere, which consists of fifty-six loudspeakers. This offers the ability of listening to the radiated sound using a higher-order ambisonic synthesis of a 'virtual' source before it is manufactured. To qualitatively evaluate the performance of the proposed procedure, the transfer function of the synthesized horn is measured around the listening point within the spatialization sphere. Key words : 3D sound synthesis, high-order ambisonics, directive sources, waveguides, horn, multimodal method

HRTF individualization using ear images from 3D meshes

The head-related transfer function (HRTF) describes how a human receives sound from different directions in space. It is unique to each listener. The best way to obtain an individualized HRTF is through direct measurement, but it requires sophisticated laboratory equipment and long measurement times. This work proposes an approach for HRTF individualization employing ear images from 3D meshes using a deep neural network and spherical harmonics transform (SHT). The method relies on the HUTUBS dataset, including 3D meshes and HRTFs. The model uses ear images to predict a low-dimensional representation of the HRTF. Initially, ten images of the right ear of each 3D mesh are taken in different positions. The model consists of two main parts. The first is a convolutional neural network (CNN), which is used to extract features from the ear images. The second part learns from the feature map obtained and predict the spherical harmonic coefficients. Finally, the individualized HRTF is obtained through inverse SHT. The performance of the method is evaluated by computing the log-spectral distortion (LSD) between the measured HRTF and the predicted one. The results show favorable LSD values compared to other models addressing the same problem.

3D audio recording of complex underwater sound fields

In this study, we developed a method to capture underwater sound fields and reproduce them for human listeners, providing an immersive 3D auditory experience. We recorded the sound field with an array of four hydrophones arranged in a tetrahedral structure. The side length of the tetrahedron (around 70 cm) was a compromise between capturing accurately the low frequency pressure gradients and spatial aliasing effects at higher frequencies. Subsequently, we transformed the signals into the domain of spherical harmonics facilitating a representation as first-order Ambisonics that is independent of the recording. For playback on devices such as headphones or loudspeaker systems, we used established decoding routines that enhance the immersive audio experience with techniques such as parametric decoding and diffuse sound field equalization. During an expedition to northern Norway, we recorded the underwater soundscape in deep fjords with vocalizing killer whales and singing humpback whales. The resulting binaural renderings for headphones of these sound scenes are astonishing and sparked interest among scientist, artists and general public. This research also helps to enhance our understanding and perception of these intelligent marine mammals and their unique acoustic underwater world, ultimately aiding in better protecting them from noise pollution.

Reconstruction of directional sources using physics-informed neural networks

In the context of augmented/virtual reality applications, replicating the directivity of sound sources is a critical step to ensure a high-quality immersive experience. However, accurate directivity measurements require complex experimental procedures involving a high number of microphones. In these scenarios, spatial interpolation allows to reconstruct the sound field with fewer sensors, which can be sparsely distributed around the source. The literature offers several interpolation schemes for sound field reconstruction, although they exhibit limited accuracy due to the lack of knowledge on the underlying physics. Recently, Physics-Informed Neural Networks (PINN) have been used to directly solve partial differential equations and thus better predict the evolution of dynamical systems. In this manuscript, we propose a spatial interpolation scheme based on PINN for the reconstruction of sound source directivity. PINN are used to solve the Helmholtz equation starting from acquisitions of the directivity over a limited set of sparsely distributed points. Results are compared with the outcomes of well-established methods based on Spherical Harmonics Decomposition (SH), Thin Plate Pseudo-Splines Interpolation (SPLI) and Piece-wise Linear Spherical Triangular Interpolation (TRI). PINN prove to provide better reconstruction of the directivity with respect to the proposed comparison methods, even when dealing with highly sparse sampling grids.

Conjugate points along spherical harmonics

Utilizing structure constants, we present a version of the Misiolek criterion for identifying conjugate points. We propose an approach that enables us to locate these points along solutions of the quasi-geostrophic equations on the sphere S2. We demonstrate that for any spherical harmonics Ylm with 1≤|m|≤l, except for Y1±1 and Y2±1, conjugate points can be determined along the solution generated by the velocity field elm=∇⊥Ylm. Subsequently, we investigate the impact of the Coriolis force on the occurrence of conjugate points. Moreover, for any zonal flow generated by the velocity field ∇⊥Yl10, we demonstrate that proper rotation rate can lead to the appearance of conjugate points along the corresponding solution, where l1=2k+1.∈N Additionally, we prove the existence of conjugate points along (complex) Rossby-Haurwitz waves and explore the effect of the Coriolis force on their stability.

Analysis of a Combined Spherical Harmonics and Discontinuous Galerkin Discretization for the Boltzmann Transport Equation

Abstract In [L. Bourhrara, A new numerical method for solving the Boltzmann transport equation using the PN method and the discontinuous finite elements on unstructured and curved meshes, J. Comput. Phys. 397 2019, Article ID 108801], a numerical scheme based on a combined spherical harmonics and discontinuous Galerkin finite element method for the resolution of the Boltzmann transport equation is proposed. One of its features is that a streamline weight is added to the test function to obtain the variational formulation. In the present paper, restricting our attention to the advective part of the Boltzmann equation, we prove the convergence and provide error estimates of this numerical scheme. To this end, the original variational formulation is restated in a broken functional space. The use of broken functional spaces enables to build a conforming approximation, that is the finite element space is a subspace of the broken functional space. The setting of a conforming approximation simplifies the numerical analysis, in particular the error estimates, for which a Céa’s type lemma and standard interpolation estimates are sufficient for our analysis. For our numerical scheme, based on ℙ k {\mathbb{P}^{k}} discontinuous Galerkin finite elements (in space) on a mesh of size h and a spherical harmonics approximation of order N (in the angular variable), the convergence rate is of order 𝒪 ⁢ ( N - t + h k ) {\mathcal{O}(N^{-t}+h^{k})} for a smooth solution which admits partial derivatives of order k + 1 {k+1} and t with respect to the spatial and angular variables, respectively. For k = 0 {k=0} (piecewise constant finite elements) we also obtain a convergence result of order 𝒪 ⁢ ( N - t + h 1 2 ) {\mathcal{O}(N^{-t}+h^{\frac{1}{2}})} . Numerical experiments in one, two and three dimensions are provided, showing a better convergence behavior for the L 2 {L^{2}} -norm, typically of one more order, 𝒪 ⁢ ( N - t + h k + 1 ) {\mathcal{O}(N^{-t}+h^{k+1})} .

$\mathrm{BC}_{2}$-Type Multivariable Matrix Functions and Matrix Spherical Functions

Matrix spherical functions associated to the compact symmetric pair (\mathrm{SU}(m+2), \mathrm{S}(\mathrm{U}(2)\times \mathrm{U}(m))) , m\geq 2 , having a reduced root system of type \mathrm{BC}_{2} , are studied. We consider an irreducible K -representation (\pi,V) arising from the \mathrm{U}(2) -part of K , and the induced representation \mathrm{Ind}_{K}^{G} \pi splits multiplicity-free. The corresponding spherical functions, i.e. \Phi \colon G \to \mathrm{End}(V) satisfying \Phi(k_{1}gk_{2})=\pi(k_{1})\Phi(g)\pi(k_{2}) for all g\in G , k_{1},k_{2}\in K , are studied by examining certain leading terms which involve hypergeometric functions. This is done explicitly using the action of the radial part of the Casimir operator on these functions and their leading terms. To suitably grouped matrix spherical functions we associate two-variable matrix orthogonal polynomials giving a matrix analogue of Koornwinder’s 1970s two-variable orthogonal polynomials, which are Heckman–Opdam polynomials for \mathrm{BC}_{2} . In particular, we find explicit orthogonality relations with the matrix polynomials being eigenfunctions to an explicit second-order matrix partial differential operator. The scalar part of the matrix weight is less general than Koornwinder’s weight.

A single layer representation of the scattered field for multiple scattering problems

The scattering of scalar waves by a set of scatterers is considered. It is proven that the scattered field can be represented as an integral supported by any smooth surface enclosing the scatterers. This is a generalization of the series expansion over spherical harmonics and spherical Bessel functions for spherical geometries. More precisely, given a set of scatterers, the field scattered by any subset can be expressed as an integral over any smooth surface enclosing the given subset alone. It is then possible to solve the multiple scattering problem by using this integral representation instead of an expansion over spherical harmonics. This result is used to develop an extension of the Fast Multipole Method in order to deal with subsets that are not enclosed within non-intersecting balls.

Frictional Contacts Between Rough Grains With Fractal Morphology

AbstractSurface morphology plays a crucial role in friction between two contacting geomaterial surfaces, yet many questions remain unanswered regarding how detailed frictional responses deviate from analytical solutions for smooth surfaces due to the presence of roughness. In this study, we revisit the Cattaneo‐Mindlin problem for contacts between two fractally rough elastic or elasto‐plastic spheres generated based on ultra‐high degree (e.g., up to 2,000) spherical harmonics with the corresponding wavelength less than a thousandth of the mean grain diameter. Transverse contacts are simulated by finite element method, validated by the extended Cattaneo‐Mindlin solution to full slide regime for smooth sphere contacts. Extensive simulations are conducted to study contacts between two rough spheres with various surface geometries, micro friction coefficients, normal contact distances, relative roughness, fractal dimensions, and wavelength ranges. Out results indicate that: (a) the new analytical solution can approximately predict the macro contact response except for extremely high relative roughness and narrow wavelength range; (b) deviations induced by roughness from smooth sphere contacts can be neutralized by plasticity, high normal contact interference, and high micro friction coefficient; and (c) fractal dimension impacts frictional contacts less than relative roughness. The main cause of these phenomena can be credited to the underlying microscale contact information. Contact area and stress distributions and their evolutions provide concrete evidence of these observed behavior. This work provides a pathway for applying computational contact mechanics to many geophysical fields, such as the asperity model in earthquake science and the mechanics of granular materials.