Microbeams are building blocks for many microstructures as well as microelectromechanical systems (MEMS) and cannot accurately be modelled by classical continuum theories due to size effects based on their micro-scale. These size effects can be taken into account by the so-called higher order continuum theories. Modified Strain Gradient Theory (MSGT) and Modified Couple Stress Theory (MCST) are two commonly used theories, which extend the classical local continuum theories of grade one with the introduction of additional length scale parameters. In this contribution, the variational problem governing the elasticity of higher order beam formulation and the finite element implementation based upon, are briefly introduced. To this end, well known Euler-Bernoulli beam formulation assumptions are used. The size effect for gold-micro beams is demonstrated and the length scale parameters of gold microbeams for MSGT and MCST are identified form the existing experimental data from literature for the first time. As a novel aspect, significant size effect is demonstrated for the length-scales associated with the state of the art gold microbeam structures developed for MEMS applications, which reveals the necessity of the use of higher order theories at these length scales. Advantages and drawbacks of these theories are also identified. where the strain metrics are the normal strain tensor ε, dilatation gradient vector ∇tr(ε), deviatoric stretch gradient vector η
S(1), and rotation gradient tensor χ
S. Their respective work conjugates are the stress metrics, which in turn are defined as Cauchy stress tensor σ, pressure gradient vector p, traceless part of the double stress tensor τ
S(1), and couple stress tensor m S . The higher order stress-strain relationships are established via elastic constants (λ, µ) and length scale parameters (l 0 , l 1 , l 2 ) asMCST is usually elaborated as the special case of MSGT in which l 0 = l 1 = 0, which further reduces the number of unknown length scale parameters in MSGT to one (l 2 ). The length scale parameters for both MSGT (l 0 , l 1 , l 2 ) and MCST (l 2 ) are hereinafter referred to as l. For further sections it will be assumed that l 0 = l 1 = l 2 , hence l 0 will be used to refer to the length scale parameter for MSGT, whereas l 2 will be used for MCST.
Quantification of Length Scale ParametersFinite element codes are developed using variational principles. Then, experiments of Espinosa et al. [4] are simulated with these codes in order to come up with unknown length scale parameters for gold. An error parameter as the L2-norm of the residual vector is defined and sequential runs are performed for various values elastic moduli E and corresponding values of l that minimize the error function are determined. The minimum and maximum values of E are chosen according to the upper and lower limits reported in literature, i.e. 140 GPa and 20 GPa respectively. The length scale parameters yielding the minimum error for macroscopic elastic modulus of gold i.e. E=80 GPa are found as l 0 =3.6...