First, I would recommend to go through some literature on low noise signal processing, just get well ground in the ideology there:
The questions asked there are usually something like: can the energy below Boltzmann kT/2 be detected? - Yes, it can be.
And commercial charge sensitive amplifiers routinely demonstrate noise well bellow kT/2. Why? Usually, the signal energy is released in much shorter time than thermal (dissipation) time constant of the system.
Can you detect 2.7K of microwave background using antenna and cabling at 300K? - Yes, as far as they have very losses.
But if you have 10% power loss in the antenna fider cable that will give you extra ~30K of the noise temperature.
Can very low signal be detected, e.g. 50pV - yes, depending on the effective bandwidth or same to say - integration time.
So systems having "low dissipation" (keyword corresponding to some books) can be used to detect very low signals, well below thermal noise (the important parameter being duration-bandwidth and characteristics of thermal= disspative coupling).
Sorry for long wandering around.
Coming to you case - I guess you right regarding parallel mode resonance. I guess you should connect your device to low noise RF MOSFET of capacitance comparable to static capacitance of your crystal. I would consider BF1207 Dual N-channel dual gate MOSFET or similar as the preamplifier. You may even consider several of them in parallel. One may also think of GaAs FET but their low frequency performance is questionable and they much harder to work with.
One may also consider JFETs like Moxtek's or low cost BF862 but at 4MHz they have no advantage over MOSFET (at lower frequency they would be preferred because of lower 1/f noise).
So after some pre-amplifier (you can use e.g. low cost MiniCircuits ZFL1000LN as the second stage after the MOSFET), you can use some high frequency lock-in amplifier, e.g.
http://www.thinksrs.com/products/SR844.htm.
If it is too expensive you can down-convert your signal using some simple I/Q mixer to some very low frequency (KHz) and sample it into the PC by some low cost DAQ. The mixer can be get at low cost as eval board. You can use low cost DDS synthesizer eval board to generate the signals.
To put rough estimate, you are starting with the noise level of ~ 0.5...1nV /sqrt(hz) vs 50pV signal, that means that integrating signal for ~100 sec you can get signal/noise ~1. So after averaging for an hour or so you may see clear signal.Important to note - the signal have to stay coherent (no much phase drift).
It is very critical to minimize interconnection of the quartz to mosfet - any extra capacitance, especially lossy one (like common FR-4 PCB material) will be a killer.