This study examines the nonlinear dynamics in tapping-mode atomic force microscopy (AFM) with tip-surface interactions that include van der Waals and Derjaguin-Muller- € Toporov contact forces. We investigate the periodic solutions of the hybrid system by performing numerical pseudo-arclength continuation. Through the use of bifurcation locus maps in the set of parameters of the discontinuous model, the overall dynamical response scenario is assessed. We demonstrate the influence of various dissipation mechanisms that are related with the AFM touching or lacking contact with the sample. Local and global analyses are used to investigate the stability of the stable solution in the repulsive regime. The impacting nonsmooth dynamics framed within a higher-mode Galerkin discretization is able to capture windows of irregular and complex motion.

Quantifying the nanomechanical properties of soft-matter using multi-frequency atomic force microscopy (AFM) is crucial for studying the performance of polymers, ultra-thin coatings, and biological systems. Such characterization processes often make use of cantilever's spectral components to discern nanomechanical properties within a multi-parameter optimization problem. This could inadvertently lead to an over-determined parameter estimation with no clear relation between the identified parameters and their influence on the experimental data. In this work, we explore the sensitivity of viscoelastic characterization in polymeric samples to the experimental observables of multi-frequency intermodulation AFM. By performing simulations and experiments we show that surface viscoelasticity has negligible effect on the experimental data and can lead to inconsistent and often non-physical identified parameters. Our analysis reveals that this lack of influence of the surface parameters relates to a vanishing gradient and non-convexity while minimizing the objective function. By removing the surface dependency from the model, we show that the characterization of bulk properties can be achieved with ease and without any ambiguity. Our work sheds light on the sensitivity issues that can be faced when optimizing for a large number of parameters and observables in AFM operation, and calls for the development of new viscoelastic models at the nanoscale and improved computational methodologies for nanoscale mapping of viscoelasticity using AFM.

This paper investigates the nonlinear dynamics in tapping-mode atomic force microscopy (AFM) with tip-surface interactions that include Van der Waals and Derjaguin-Müller-Toporov contact forces. We study the periodic solutions of the hybrid system by performing numerical pseudo-arclength continuation. The overall dynamical response scenario is evaluated via bifurcation loci maps in the set of parameters of the discontinuous model. We showcase the influence of different dissipation mechanisms activated when the AFM is in contact or out-of contact with the sample. The robustness of the stable solution in the repulsive regime is studied via local and global analyses. The impacting non-smooth dynamics framed within a higher-mode Galerkin discretization is able to capture windows of irregular and complex motion.

Increasing the signal-to-noise ratio in dynamic atomic force microscopy plays a key role in nanomechanical mapping of materials with atomic resolution. In this work, we develop an experimental procedure for increasing the sensitivity of higher harmonics of an atomic-force-microscope cantilever without modifying the cantilever geometry but instead by utilizing dynamical mode coupling between its flexural modes of vibration. We perform experiments on different cantilevers and samples and observe that via nonlinear resonance frequency tuning we can obtain a frequency range where strong modal interactions lead to 7-fold and 16-fold increases in the sensitivity of the 6th and 17th harmonics while reducing sample indentation. We derive a numerical model that captures the observed physics and confirms that nonlinear mode coupling is the reason for the increase of the amplitude of higher harmonics during tip-sample interactions.

We report the first study on the thermal behavior of the stiffness of individual carbon nanotubes, which is achieved by measuring the resonance frequency of their fundamental mechanical bending modes. We observe a reduction of the Young's modulus over a large temperature range with a slope -(173±65) ppm/K in its relative shift. These findings are reproduced by two different theoretical models based on the thermal dynamics of the lattice. These results reveal how the measured fundamental bending modes depend on the phonons in the nanotube via the Young's modulus. An alternative description based on the coupling between the measured mechanical modes and the phonon thermal bath in the Akhiezer limit is discussed.

In this work, we perform a comprehensive analysis of the robustness of attractors in tapping mode atomic force microscopy. The numerical model is based on cantilever dynamics driven in the Lennard–Jones potential. Pseudo-arc-length continuation and basins of attraction are utilized to obtain the frequency response and dynamical integrity of the attractors. The global bifurcation and response scenario maps for the system are developed by incorporating several local bifurcation loci in the excitation parameter space. Moreover, the map delineates various escape thresholds for different attractors present in the system. Our work unveils the properties of the cantilever oscillation in proximity to the sample surface, which is governed by the so-called in-contact attractor. The robustness of this attractor against operating parameters is quantified by means of integrity profiles. Our work provides a unique view into global dynamics in tapping mode atomic force microscopy and helps establishing an extended topological view of the system.

Despite numerous theoretical investigations on the mechanical properties of graphene, an accurate identification of its material behavior is still unattained. One hypothesis for this uncertainty is that modeling graphene as a static membrane cannot describe the strong coupling between mechanics and thermodynamics of this structure. Therefore, characterization methods built upon static models could not capture these effects. In this paper, we propose a new method for building a reduced order model for the dynamics of thermalized graphene membranes. We apply the proper orthogonal decomposition algorithm on time responses obtained from molecular dynamics simulations. As a result, a set of orthogonal modes is obtained which are then employed to build a reduced order model. The proposed model can describe the motion of the suspended graphene membrane over the whole spatial domain accurately. Moreover, due to its computational efficiency, it is more versatile for exploring the nonlinear dynamics of the system. This model is then employed for studying the nonlinear dynamics of graphene membranes at large amplitudes to extract Young's modulus. The obtained Young's modulus incorporates the effects of nano-scaled thermally induced dynamic ripples and hence, is temperature and size dependent. Our proposed atomistic modal order reduction method provides a framework for studying the dynamics and extracting the mechanical properties of other nano-structures at the molecular level.

This paper investigates the complex bifurcation scenario of electrically-actuated circular micro-plates subjected to differential pressure. Our analysis deals with the primary, secondary and ultimate saddle-node bifurcation points responsible for the device snap-through and addresses the pressure range in which the robustness of the two main stable configurations is undermined by minor coexisting attractors. By making use of basins of attraction and integrity profiles, safe dynamical regions of motion are evaluated with respect to the applied DC voltage. It is found that period-doubling bifurcations are accountable for a sensible reduction in the dynamical integrity, small variations of the DC voltage largely modify the response of the pressure sensor.

Diamond is a highly desirable material for state-of-the-art micro-electromechanical (MEMS) devices, radio-frequency filters and mass sensors, due to its extreme properties and robustness. However, the fabrication/integration of diamond structures into Si-based components remain costly and complex. In this work, a lithography-free, low-cost method is introduced to fabricate diamond-based micro-resonators: a modified home/office desktop inkjet printer is used to locally deposit nanodiamond ink as ∅50–60 µm spots, which are grown into ≈1 µm thick nanocrystalline diamond film disks by chemical vapor deposition, and suspended by reactive ion etching. The frequency response of the fabricated structures is analyzed by laser interferometry, showing resonance frequencies in the range of ≈9–30 MHz, with Q-factors exceeding 10⁴, and (f₀ × Q) figure of merit up to ≈2.5 × 10¹¹ Hz in vacuum. Analysis in controlled atmospheres shows a clear dependence of the Q-factors on gas pressure up until 1 atm, with Q ∝ 1/P. When applied as mass sensors, the inkjet-printed diamond resonators yield mass responsivities up to 981 Hz fg⁻¹ after Au deposition, and ultrahigh mass resolution up to 278 ± 48 zg, thus outperforming many similar devices produced by traditional top-down, lithography-based techniques. In summary, this work demonstrates the fabrication of functional high-performance diamond-based micro-sensors by direct inkjet printing.

Global properties of Multi-Degrees-of-Freedom (M-DoF) systems, in particular phase space organization, are largely unexplored due to the computational challenge requested to build basins of attraction. To overcome this problem, various techniques have been developed, some trying to improve algorithms and to exploit high speed computing, others giving up to possibility of having the exact phase space organization and trying to extract major information on a probability base. Following the last approach, this work exploits the method of “basin stability” (Menck et al., 2013) in order to drastically reduce the numerical effort. The probability of reaching the attractors is evaluated using a reasonable number of trials with random initial conditions. Then we investigate how this probability depends on particular generalized coordinate or a pair of coordinates. The method allows to obtain information about the basins compactness and reveals particular features of the phase space topology. We focus the study on a 2-DoF multistable paradigmatic system represented by a parametric pendulum on a moving support and model of a Church Bell. The trustworthiness of the proposed approach is enhanced through the comparison with the classical computation of basins of attraction performed in the full range of initial conditions. The proposed approach can be effectively utilized to investigate the phase space in multidimensional nonlinear dynamical systems by providing additional insights over traditional methods.

The bending rigidity of two-dimensional (2D) materials is a key parameter for understanding the mechanics of 2D NEMS devices. The apparent bending rigidity of graphene membranes at macroscopic scale differs from theoretical predictions at micro-scale. This difference is believed to originate from thermally induced dynamic ripples in these atomically thin membranes. In this paper, we perform modal analysis to estimate the effective bending rigidity of graphene membranes from the frequency spectrum of their Brownian motion. Our method is based on fitting the resonance frequencies obtained from the Brownian motion in molecular dynamics simulations, to those obtained from a continuum mechanics model, with bending rigidity and pretension as the fit parameters. In this way, the effective bending rigidity of the membrane and its temperature and size dependence, are extracted, while including the effects of dynamic ripples and thermal fluctuations. The proposed method provides a framework for estimating the macroscopic mechanical properties in other 2D nanostructures at finite temperatures.

Physical characteristics such as mass and stiffness of biological objects are emerging as new markers for severe diseases. Micromechanical resonators can be used to quantify multiple of these characteristics simultaneously. In this paper, we propose a methodology that utilizes higher flexural modes of vibration to perform simultaneous characterization of the density and elastic modulus of adsorbates. To demonstrate this concept, a polymeric block with a known dimension and anisotropy is written directly on the cantilever surface using a two-photon polymerization technique and characterised by modal analysis. Our method captures the effective bending stress exerted by non-isotropic materials which is masked in the atomic force microscopy indentation technique.

The accuracy of measurements in Amplitude Modulation Atomic Force Microscopy (AFM) is directly related to the geometry of the tip. The AFM tip is characterized by its radius of curvature, which could suffer from alterations due to repetitive mechanical contact with the surface. An estimation of the tip change would allow the user to assess the quality during imaging. In this work, we introduce a method for tip radius evaluation based on the nonlinear dynamic response of the AFM cantilever. A nonlinear fitting procedure is used to match several curves with softening nonlinearity in the noncontact regime. By performing measurements in this regime, we are able to maximize the influence of the tip radius on the AFM probe response, and this can be exploited to estimate with good accuracy the AFM tip radius.

Numerical integrations represent a time-consuming element in the long-term dynamics analysis of mechanical systems. This limits the resolution of the computations and the dimension of the system to be investigated numerically. In fact, even pushing memory resources to their thresholds, only few tools can deal with higher-dimensional systems. This work illustrates, in a preliminary manner, the results that can be obtained reducing the aforementioned constraints thanks to the implementation of algorithms based on a parallel computing approach. In particular, by focusing on basins of attraction, four applications are discussed. i) The full domain of attraction for a four-dimensional (4D) system describing a linear oscillator coupled with a nonlinear absorber is calculated. ii) The variation of a safe basin with respect to the system dimension is then analyzed. It is highlighted how 4D and 3D analyses provide more confident results with respect to 2D analyses. iii) The parametric variation of a 2D system with a reduced step is performed by building a 3D representation which allows to highlight a smooth transition between the states. iv) A convergence study of a basin of attraction resolution is carried out. The integrity factor is used as a comparison measure.

Linear and non-linear vibrations of a U-shaped hollow microcantilever beam filled with fluid and interacting with a small particle are investigated. The microfluidic device is assumed to be subjected to internal flowing fluid carrying a buoyant mass. The equations of motion are derived via extended Hamilton's principle and by using Euler-Bernoulli beam theory retaining geometric and inertial non-linearities. A reduced-order model is obtained applying Galerkin's method and solved by using a pseudo arc-length continuation and collocation scheme to perform bifurcation analysis and obtain frequency response curves. Direct time integration of the equations of motion has also been performed by using Adams-Moulton method to obtain time histories and analyze transient cantilever-particle interactions in depth. It is shown that exploiting near resonant non-linear behavior of the microcantilever could potentially yield enhanced sensor metrics. This is found to be due to the transitions that occur as a matter of particle movement near the saddle-node bifurcation points of the coupled system that lead to jumps between coexisting stable attractors.