Data in Figure 1e and Figure 4a:
Linear simulations of scattering and absorption cross-section for Si(+) and Si(-) nanohelices, oriented in x-,y- or z-direction and excited with left or right circular polarizations (LCP or RCP, respectively) are saved in files:
-Abs_Si_3turns_mesh2_negative
-Abs_Si_3turns_mesh2_positive
-Sc_Si_3turns_mesh2_negative
-Sc_Si_3turns_mesh2_positive

Abs (Sc) denotes absorption (scattering) spectra.
Positive (negative) denotes Si(+) (Si(-)) nanohelices.

Each file is a matrix of 6x100 elements, where 100 points correspond to 100 wavelengths between 200 and 1000 nm.
Lines are organized as follows:
- Line 1: x-oriented nanohelix, LCP excitation
- Line 2: x-oriented nanohelix, RCP excitation
- Line 3: y-oriented nanohelix, LCP excitation
- Line 4: y-oriented nanohelix, RCP excitation
- Line 5: z-oriented nanohelix, LCP excitation
- Line 6: z-oriented nanohelix, RCP excitation
Additional file lambda_200_1000 contains the wavelength points.

Total extinction spectra for two polarization is obtained by averaging the data for a random orientation.
CD and g factors are extracted from these extinction spectra defined in the formulas in the manuscript. 

Data in Figure 4b:
Four-dimensional complex electromagnetic field extraction from Lumerical and imported into MENP solver in MATLAB.
There are six matrices that represent x, y, and z-oriented nanohelices for LCP and RCP excitation in a broadband wavelength range.
Titles of the files contain information on the nanohelix orientation, while p90 (m90) stands for RCP (LCP) excitation.

Data in Figure 4c: 151x151x151 matrices of the absorption density were recorded for Si(+) sample excited at 730 nm with LCP or RCP;
to be plotted in MATLAB by using sliceplot and hot colormap.
Linear factor g_a_b_s is then calculated as defined in the equation of the manuscript.

Data in Figure 4d: for second harmonic simulations, electric field components and the intensity are recorded in xy and xz cross-sections,
under 730 nm excitation of the z-oriented nanohelix, for LCP and RCP excitation.
Files are named as follows:
E2 xz 365nm zor p90 data contains electric field intensity (E2) in the xz cross section, at 365 nm, for z-oriented nanohelix, excited with RCP (p90) at 730nm.
Ez xy 365nm zor m90 data contains real component of Ez in the xy cross section, at 365 nm, for z-oriented nanohelix, excited with LCP (m90) at 730nm.
Nonlinear g factor g_NL is then calculated from the normalized difference of the electric field intensity distributions in the xy and xz planes.